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ABSTRACT 


The  Growth  Potential  of  Corona  Discharges  from  Aircraft  Flying  in  Precipitation 

The  behaviour  of  positive  streamers  in  fields  near  and  above  the  stability 
field  has  been  studied  in  quasi-uniform  fields  of  up  to  18cm  in  extent  with  the  aim 
of  producing  a  model  of  the  behaviour  so  that  predictions  can  be  made  for  much 
greater  distances.  In  particular  the  evolution  in  air  of  free  charge  produced  by 
the  propagation  has  been  monitored  and  the  effect  of  pressure  and  propagation 
distance  on  this  evolution  studied.  A  frequency  (P)  or  coefficient  (?)  has  been  de¬ 
fined  and  measured  which  characterises  the  growth  of  free  charge.  The 
Schockley-Ramo  theorem  has  been  used  to  analyse  the  data  which  not  only  pro¬ 
vides  values  for  other  relevant  parameters  such  as  attachment  coefficient  and 
streamer  velocity  but  also  provides  a  novel  way  of  observing  prebreakdown 
streamer  modes.  A  model  has  been  constructed  which  adequately  accounts  for  the 
observed  behaviour  in  small  gaps  but  which  predicts  unrealistically  large  amounts 
of  free  charge  in  large  gaps  (>lm).  This  suggests  that  some  other  mechanism 
might  predominate  in  this  case;  an  increase  in  electron  density  leading  to  channel 
thermalisation  is  suggested  as  a  possibility. 
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Section  1:  Introduction 


The  most  remarkable  single  event  in  atmospheric  physics  is  surely  the  lightning 
discharge  so  that  it  comes  as  no  surprise  that  the  study  of  the  event  has  an  impressive 
history  together  with  an  equally  impressive  list  of,  in  many  cases,  legendary  investigators 
(Franklin,  C  T  R  Wilson,  Schonland  and  others).  Over  the  years  an  enormous  quantity  of 
data  has  been  garnered  from  numerous  programmes  conducted  around  the  world,  most 
notably  in  those  regions  where  the  lightning  stroke  is  an  devastating  as  it  is  spectacular. 
In  view  of  this,  what  is  surprising  is  that  the  entire  mechanism  from  the  initial  charge 
separation  through  onset  and  leader  development  to  the  final  arc  is  not  satisfactorily 
understood.  For  example,  there  is  still  a  certain  controversy  surrounding  even  the 
physics  of  the  charge  separation  whilst  the  measured  electric  fields  in  the  vicinity  of 
thunderclouds  (l-2kV/cm)  are  substantially  lower  than  those  encountered  in  laboratory 
long  gap  breakdowns.  The  investigation  reported  here  was,  and  still  is,  concerned  with 
the  onset  of  the  discharge.  That  is  to  say,  given  the  thundercloud  environment  -  with 
fields  of  measured  magnitudes  and  precipitation  of  various  sorts  -  together  with  a  know¬ 
ledge  of  laboratory  breakdown  processes  -  could  the  onset  of  the  discharge  be  reasonably 
well  predicted?  There  is  an  additional  question  involving  aircraft  as  will  shortly  become 
clear. 

Most  breakdowns  originate  with  corona  and  its  subsequent  development  or  growth; 
the  exception  being  breakdown  between  parallel  electrodes  involving  a  strictly  uniform 
field.  In  order  to  obtain  corona,  which  is  a  property  of  high  fields,  in  the  low  field  of  a 
thundercloud,  demands  the  assumption  that  precipitation  be  present  -  a  not  too  severe 
limitation  for  a  thundercloud!  The  geometrically  enhanced  electric  fields  at  the  surface 
of  the  precipitation  would  then  reach  the  level  required  for  corona  onset.  However,  it 
turns  out  that  the  fields  required  to  generate  corona  from  ice  and  water  particles  of  the 
anticipated  size  are  still  significantly  larger  than  those  commonly  encountered  in  a 
thundercloud.  The  exception  here  is  the  long  water  filament  drawn  out  when  water  drop¬ 
lets  collide  -  an  event  which  combines  a  low  onset  field  with  relative  rarity;  desirable 
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features  when  seeking  the  cause  of  what  is  after  ail  a  triggered  discharge.  Of  course, 
aircraft  offer  exactly  the  same  characteristics  -  indeed  the  onset  field  is  likely  to  be 
even  lower  -  so  if  the  argument  that  corona  is  a  discharge  prerequisite  is  acceptable  it 
follows  that  aircraft  can  trigger  a  lightning  flash.  It  is  probably  worth  adding  at  this 
point  that  the  overwhelming  experimental  evidence  does  currently  point  to  a  positive 
confirmation  of  aircraft  triggered  discharges.  Mazur,  for  example,  has  identified,  using 
radar,  the  source  of  several  '  hes  as  an  aircraft  flying  through  a  thunderstorm;  rocket 
triggered  discharges  provide  additional  support.  When  the  investigation  reported  here 
was  undertaken  in  1984  this  was  not  the  common  view  (e.g.  Clifford). 

Of  course,  corona  is  not  the  sole  requirement.  Generally  speaking,  corona  pro¬ 
duction  will  tend  to  dissipate  the  field  which  produces  it.  To  understand  corona  driven 
breakdown  from  essentially  isolated  conductors  (no  discrete  power  source)  a  distinction 
must  be  made  between  the  two  forms  of  corona  -  negative  and  positive  -  both  of  which 
will  usually  occur  simultaneously  at  opposite  extremes  of  such  a  conductor  located  in  a 
suitable  field.  Note  this  is  not  the  same  as  monopolar  corona  from  a  charged  conductor 
with  or  without  an  ambient  field.  The  key  difference  lies  in  the  concept  of  a  stability 
field;  this  is  the  minimum  ambient  field  in  which  the  corona  will  propagate.  Negative 
corona  requires  a  much  higher  field  for  propagation  than  positive  corona.  The  prop¬ 
agating  corona  is  usually  termed  a  streamer  and  the  positive  streamer  stability  field  is 
some  4kV/cm  only  at  atmospheric  pressure.  Since  this  field  is  pressure  dependent  then 
the  lower  pressures  in  a  thundercloud  lower  t,>e  stability  field  to  values  not  untypical  of 
measured  thundercloud  fields.  It  is  for  this  reason  that  the  growth  of  the  positive 
streamers  for  F ^  ^stability  ‘s  believed  to  be  °f  key  importance. 

This  investigation  then  is  concerned  almost  entirely  with  the  behaviour  of  such 
streamers.  Establishing  suitable  conditions  is  not  as  easy  as  would  appear  because  of  the 
large  voltages  involved.  Even  1MV  would  only  allow  propagation  distances  of  2.5m  to  be 
studied  and  then  because  of  the  size  of  the  electrodes  required  to  guarantee  quasi- 
uniform  fields  (">-7.5m  in  diameter)  only  measurements  at  atmospheric  pressure  could  be 
made.  The  investigation  has  been  confined  for  these  reasons  to  gaps  of  only  20cm  but 
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with  the  electrode  structure  contained  in  a  vacuum  chamber  which  allows  the  pressure  to 
be  varied  from  200  torr  to  760  torr  -  below  200  torr  streamers  are  difficult  to  generate. 
An  outline  of  the  arrangement  is  contained  in  Section  2.  The  theoretical  treatment  re¬ 
quired  to  interpret  the  results  is  in  Section  3  with  some  typical  results  in  Section  4.  Most 
of  the  data  obtained  is  contained  in  Section  5  with  some  discussion  of  these  results  in 
Section  6  and  a  broad  conclusion  presented  in  Section  7. 
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Section  2.  Experimental  Arrangement 

2.1  Electrode  System 

The  arrangements  consists  of  two  plane  parallel  electrodes.  A  third  needle  elect¬ 
rode  is  inserted  in  a  hole  in  the  plane  anode,  the  needle  being  earthed.  The  redesigned 
upper  electrode  (Fig.2.1)  is  a  25cm  diameter  aluminiun  toroid,  with  a  flat  face.  The  7cm 
diameter  edge  profile  is  intended  to  prevent  natural  corona  inception  from  the  cathode 
by  keeping  the  field  divergence  at  the  edge  low  and  the  12cm  clearance  from  the  GRP 
vessel  walls  prevents  surface  discharges  along  these  to  the  lower  electrode. 

This  is  a  42cm  diameter  copper  plate,  5mm  thick,  and  stands  on  three  spacers, 
resting  on  adjustable  feet.  Changing  spacers  allows  electrode  separations  from  5  to 
20cm  to  be  used.  A  6mm  diameter  stainless  steel  needle  protrudes  through  a  2cm  hole  in 
the  centre  of  the  lower  electrode,  set  by  line-of-sight  against  a  block,  to  be  12mm  above 
the  lower  electrode  surface  (for  easily  repeatable  consistency)  whatever  the  electrode 
spacing.  Both  are  earthed  via  (nominally)  5QJlresistors  (see  schematic  Fig. 2. 2).  The 
earth  point  is  a  copper  strip,  connected  to  a  mains  earth.  Subsequent  measurement  of 
the  resistance  of  the  signal  resistors  used  showed  them  to  be  55  and  56A  for  the 
conduction  and  displacement  circuit  respectively.  The  cables  used  to  connect  with  the 
oscilloscope  were  found  by  a  resonant  method  to  have  a  dynamic  impedance  of  76.fi.  All 
these  values  are  built  into  the  computer  program  which  finds  the  currents  from  the  volt¬ 
age  waveforms. 

Up  to-lOOkV  is  applied  to  the  upper  electrode  either  directly  or  via  a  200  J1  high 
voltage  resistor  from  a  Brandenburg  T-series  supply.  The  high  voltage  charging  resistor 
was  incorporated  to  prevent  the  formation  of  an  arc. 

Assuming  the  physical  circuit  between  the  HT  supply  and  the  electrodes  to  be 
circular  and  of  radius  40cm,  its  inductance,  L  =  |i(Jxr/2'v0.8pH.  The  capacitance  of  the 
electrodes,  C  =CQA/d  ~10pF  and  hence  the  time-constant  for  recharging  the  cathode 
directly, 1=  LC~3ns. 

If  the  cathode  charge  O  =  CV  ~15yC,  is  signficantly  depleted  inside  ^Jns,  the 
voltage  falls,  choking  off  the  discharge.  The  risetime  of  a  corona  pulse  may  be^-100ns 


and  a  time -to -spar*  is  typically  ~lps;  hence  corona  and  arc  processes  can  develop  unin¬ 
hibited  if  the  charging  resistor  is  left  out  of  the  circuit.  If  the  cathode  must  recharge 
via  the  200MJ1  resistor,  then  1  =CR  ~48ms  and  a  high  current  arc  cannot  be  sustained. 

2.2  Digitizing  System 

The  voltages  across  the  signal  resistors  are  taken  down  coaxial  cable  and  digitized 
by  a  Philips  PM315  125MHz  digital  storage  oscilloscope  (256  sample  memory).  Typical 
voltage  pulses  obtained  with  a  Tektronix  300MHz  bandwidth  analogue  oscilloscope  show 
some  oscillation  taut  no  short  timescale  events.  Though  the  PM315  has  an  internal  60MHz 
filter,  no  data,  therefore,  is  lost  on  account  of  this.  It  became  apparent  that  the  two 
channels  were  imperfectly  synchronized.  The  rise  in  the  conduction  current  circuit 
occurred  one  sample  later  than  the  rise  in  the  diplacement  circuit.  This  has  quite  severe 
consequences  for  the  analysis  routines  used,  as  discussed  later. 

The  high  speed  digital  oscilloscope  can  be  pre-triggered.  The  data  is  continually 
read  in  and  once  a  trigger  is  experienced,  the  read-in  of  data  stops.  Hence  data  can  be 
gathered  from  before  the  triggering  point.  This  is  an  improvement  over  traditional  oscil¬ 
loscope  measurements  where  the  very  first  portion  of  the  event  is  always  absent  from 
the  records.  (The  streamer  inception  is  irregular  and  so  repetitive  sampling  techniques 
cannot  be  used.) 

The  data  recorded  by  the  digital  scope  were  recorded  using  a  Z80-based  micro¬ 
computer  via  a  GPIB-418  IEEE  standard  interface.  The  Pascal  source  code  of  the  algor¬ 
ithm  for  handling  this  is  reprinted  in  the  Appendix  but  the  implementation  is  heavily 
machine-specific  and  the  code  is  particularly  inefficient.  This  aspect  of  the  measuring 
system  leaves  considerable  scope  for  improvement. 

The  (integer)  data  from  the  scope  are  stored  in  a  disk  file  and  converted  to  an 
array  of  real  data  representing  the  currents  calculated  using  the  measured  values  of  the 
signal  resistances  and  dynamic  cable  impedances  in  parallel. 
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2.3  Experiments  Carried  Out 

The  needle  projected  by  12mm  above  the  plane  to  generate  natural  corona.  In  this 
polarity,  positive  streamers  were  produced  in  the  high  field  induced  at  the  needle  end  and 
propagate  into  the  gap.  The  needle  is  at  earth  potential  and  this  falls  quite  rapidly  to  the 
potential  which  would  have  been  experienced  at  this  level  in  the  plane-parallel  gap  were 
there  no  needle  present.  This  potential  drop  is  En,  though  it  must  be  stressed  that  the 
needle  itself  is  at  earth.  The  potential  drop  in  the  enhanced  field  at  the  needle  end  de¬ 
pends  on  the  field  and  needle  height,  n,  and  so  keeping  the  needle  height  constant  will  en¬ 
sure  that  the  potential  drop  in  the  gas  at  its  end  is  the  same  whatever  the  gap  length  or 
cathode  potential.  This  is  rather  different  from  the  situation  in  an  inverted  point-plane 
gap  (the  high  voltage  is  usually  applied  to  the  point),  where  the  potential  drop  at  the 
needle-end  does  depend  on  the  applied  potential. 

To  produce  streamers  is  no  easy  task,  for  the  field  produced  in  the  gas  at  the 
needle-end  must  be  high  enough  and  of  a  large  enough  extent.  The  geometry  of  the 
needle  is  crucial  and  some  considerable  time  was  spent  producing  a  sharp  needle  with  a 
snail,  uniformly  rounded  tip.  The  idea  being  to  increase  the  size  of  the  critical  volume 
by  having  a  smoothly  rounded  tip,  but  to  increase  the  field  enhancement  by  having  a 
narrow  cone  angle.  The  eventual  needle  profile,  arrived  at  by  trial  and  error  had  30° 
cone-angle,  with  a  tip  of  radius  «-20pm  (measured  with  a  travelling  microscope).  It  was 
necessary  to  electrolyse  the  tip  in  order  to  get  a  smooth  profile  of  this  small  radius. 

This  gave  satisfactory  inception  at  high  pressures  but,  as  described  in  section  4, 
not  so  good  an  inception  at  lower  pressures,  where  the  critical  avalanche  must  be  larger 
for  inception.  On  occasion,  the  inception  was  very  infrequent,  as  the  critical  volume  for 
those  particular  conditions  was  very  snail.  To  overcome  this  difficulty,  an  ionizing 
-source  (50pCurie  Americium)  was  mounted  on  the  shaft  of  the  needle  just  below  the 
hole  in  the  plane  anode  to  provide  excess  starter  electrons.  The  ionizing  source  did  not 
affect  the  clearing  time  of  the  ions  and  so  under  conditions  where  the  source  was  not 
needed,  it  could  be  left  in  the  electrode  chamber  and  have  no  effect  at  all  on  the 
discharge. 
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The  primary  aim  of  this  work  was  to  investigate  £--«  influence  of  propagation 
length  on  the  corona  characteristics.  Accordingly  strearr  ers  were  generated  and  digital 
records  of  the  induced  currents  made,  with  various  different  pressures  and  electrode 
separations.  The  data,  stored  on  disk,  was  then  analyzed  to  provide  several  fundamental 
parameters  of  the  streamer's  propagation,  described  in  section  3. 

It  should  be  noted  at  this  stage  that  as  the  data  was  digital,  with  a  full  scale  de¬ 
flection  of  8  bits,  the  display  (on  a  monitor  with  similar  resolution)  was  scaled  so  that  as 
much  information  as  possible  was  visible.  An  autoranging  algorithm  was  used  to  find  a 
scale  factor  such  that  as  much  of  the  screen  was  used  ip  and  then  rounded  up  30  that  in¬ 
convenient  numbers  were  not  involved.  The  zero  level  was  determined  by  averaging  the 
first  5  points  in  the  record,  and  the  data  from  the  oscilloscope  converted  into  a  real 
number  array  for  further  analysis.  When  the  scaling  for  the  screen  display  was  obtained 
a  zero  level  was  chosen  such  that  the  real  zero  level  corresponded  to  a  line  on  the 
graticule  drawn  on  the  screen. 

The  effect  of  this  variable  scaling  though  is  that  when  making  direct  comparison 
of  different  pulses  on  the  screen  (or  on  the  screen  dumped  copies  used  herein  for  illust¬ 
rations)  the  scaling  must  always  be  borne  in  mind. 

The  electrodes  are  contained  in  a  50cm  diameter  GRP  vessel,  which  was  pimped 
out  using  a  throttled-back  rotary  vacuum  pimp.  Gas  was  let  in  through  a  mixing  unit, 
where  the  flow  could  be  controlled.  The  pressure  established  in  the  chamber  by  this  flow 
(~10l/min),  was  monitored  by  a  simple  rotary  dial  gauge.  Room  air  was  drawn  though  a 
silica-gel  dryer,  which  was  regularly  re-filled  with  fresh  silica  gell.  Though  the  influence 
of  humidity  was  not  studied  because  of  the  time  spent  analysing  the  large  quantity  of  dry 
gas  data,  the  facilities  exist  for  doing  so. 
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Section  3:  Theoretical  Considerations 

Introduction 

An  unavoidable  paradox  arises  when  studying  streamer  growth  heK-,viour  in  n 
quasi  -unifnrn  field  in  that  the  generation  of  the  streamer  requires  highly  divergent  field 
conditions,  albeit  localised.  In  general  terms,  the  analysis  of  charge  evolution  in  a 
uniform  field  from  the  behaviour  of  externally  induced  currents  is  relatively  straight¬ 
forward;  in  two  electrode  divergent  field  geometry  progress  can  only  he  made  giver,  at 
least  a  detailed  knowledge  of  the  field  distribution,  but  even  then  a  derived  charged 
distribution  may  be  ambiguous.  The  present  geometry,  detailed  in  Section  2,  enrtnirs 
elements  of  both.  The  treatment  presented  below  indicates  flow  the  two  contributions 
(uniform  and  non-uniform)  may  be  separated  leading  to  an  estimate  of  the  free,  charge 
and  its  rate  of  generation.  The  approach  also  highlights  the  importance  of  obtaining  the 
data  in  digitised  form  even  though  some  loss  of  resolution  is  inevitable. 


7.1.2  Induced  Current 

Differentiating  Eq.  3.1  yields  the  induced  current"'^*  drawn  by  S^: 

ik*  =  -qVVk*(r).r  =  qrk*(r).  £ 
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where  f^*(r'  is  the  corresponding  per  unit  field  .'jt  £  'dimension  ),  the  charge  q 

Having  instantaneous  vrlocity  v  .  Since  potential  fields  are  superposable,  generality  is 
not  Inst  by  summing  over  all  free  charge  in  the  system  and  for  on  axial  conducting 
filament  of  linear  charge  densityA(r)  (Fig. 3. 2),  the  induced  current  is  the  line  integral: 

ik  -  J Vr)E*(r)v( r)d4  *  J K(r)jCr)^  aj) 

£  is  a  vector  along  the  filament.  The  current,  j,  in  the  filament  is  given  by  j(r)  = 
^(r’v  'the  current  is  overwhelmingly  electronic,  since  ve'S>  v^). 


3.1.3  Parallel  electrodes 

For  plane  parallel  electrodes,  Fig.  3.3  shows  V  *,  F^*,  q^"  and  i^*  for  both 
electrodes  (k:  1,2).  As  a  charge  qmoves  from  Sj  to  S^,  falls  off  linearly  and  hence 
q*  is  progessively  transferred  from  to  S^.  That  is  to  say,  induced  charge  is  conserved. 

k'l  (3.4) 


Since  the  per-unit  field,  E*  =  1/d  everywhere  and  in  the  uniform  externally  applied 
field  E  the  charge  q  has  a  constant  velocity,  the  induced  current  is, 

i*  =  qE*v  =  qv/ii 


(3.3) 


The  current  is  almost  entirely  electronic,  since  vg»  v^  and  only  free  electrons 
contribute  significantly  to  the  induced  current.  By  the  same  token,  q  is  almost  entirely 
composed  of  free  electrons;  ions  -  including  negative  ions  -  make  little  contribution  to 
the  induced  current. 

This  induced  current  is  now  defined  as  ie,  and  will  be  referred  to  as  the  electron 
current.  This  will  he  used  as  a  theoretical  concept  in  geometries  where  it  cannot  he 
defined  as  an  induced  current  in  a  particular  electrode.  In  such  geometries,  i  represents 
the  current  due  to  all  free  electrons  weighted  uniformly  (i.e.  as  if  the  per-unit  field  were 
uniform  throughout  the  electrode  volume).  So  in  the  three  electrode  system  used  in  this 
study,  the  electron  current  is  a  concept,  namely  the  current  which  would  have  heen 


incurnd  in  the  lowrr  plane  electrode  if  the  point  were  not  there. 
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3.1.4  Significance  of  i„ 

To  further  understand  i  let  a  point  charge  q  move  between  plane  parallel 
electrodes  and  consider  that:  H 

[ied t-  %  fSx/d 

t,J 

the  dipole  moment  caused  by  this  translation  per  unit  electrode  separation  between  the 
electron  and  the  positive  ion  left  behind.  In  this  special  geometry, Sx/d  is  just  the  change 
in  per-unit  potential  along  the  path£x.  In  a  general  electrode  system  then,  if  one 
electron  created  by  ionization  in  the  gap  drifts  by  $x,  and  is  neutralized  at  an  electrode 
anc!  if  the  difference  in  V*  along  this  path  is  £v*  then  the  charge  transported  to  that 
particular  electrode  is  just  5v*q,  which  will  he  less  than  q.  (This  is  not  affected  by  the 
physical  arrival  of  the  charge,  for  the  electron  will  be  neutralized  by  the  charge  built  up 
on  the  electrode,  and  no  further  current  will  flow  to  the  electrode  through  the  external 
circuit.) 

Furthermore,  the  net  charge  in  the  gap  is  now  q.  The  remainder  of  the  induced 
charge  resides  on  the  other  electrodes.  Pare  must  therefore  be  exercised  when 
considering  the  induced  currents,  particularly  where  distributed  charge  is  concerned,  as 
free  charge  is  not  conserved  (some  may  attach,  for  example);  the  integral  of  the  free 
electron  current  may  not  then  be  easy  to  interpret. 


3.1.3  External  Current 

In  general,  the  conservation  of  induced  charge  (Cq.3.4)  is  the  same  as  conservation 
of  per-unit  potential,  since  q^*  =-qV^*: 


tV: 


which  means,  for  two  electrode  systems,  that  Vj*  +  V^*  =  1  Ej*  =  '^2*  or  ‘j*  =  ~,2*’ 
which  is  hardly  surprising  as  there  is  continuity  of  current  in  the  external  circuit.  Jxiote 


that  as  i *  is  found  by  differentiating  q*,  the  sign  of  i*  is  not  related  to  the  coordinate 
system  used,  but  to  the  external  circuit.  In  this  case  the  current  into  Sj  equals  the 
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current  out  of  Sj.  Tins  highlights  a  fundamental  limitation  of  current  measure”  .Tit  in 
two-electrode  systems,  i.e.  Since  ij*  =  there  is  only  one  degree  of  freedom.  The 

electron  current  Is  defined  as  the  induced  current  in  plane  parallel  geometry.  In  point - 
plane  geometry,  the  electrons  moving  near  the  point  are  weighted  considerably  more 
heavily  than  those  further  away,  and  the  conceptually  defined  electron  current,  the 
induced  current  if  the  electrons  are  all  given  equal  weight  is  not  recoverable  frne,  the 
actual  induced  current.  It  v/ill  be  shown  that  this  is  at  least  approximately  possihle  ir 
three  electrode  geometry  where  there  is  an  extra  degree  of  freedom  and  a  portieularly 
discrete  boundary  condition. 

It  can  he  added  that  the  same  analysis  is  applicable  to  any  tv.n  groups  of  all  the 

electrodes  in  a  system.  In  the  case  of  those  used  throughout  this  study,  the  externa! 

circuit  is  formed  from  the  needle  and,  the  lower  plane  'i_  and  i  .)  and  the  cab  ode  i  . 

'  c  r:  coth 

Then  the  external  current  is  (i  ♦  i  .),  and  is  equal  to  -i  .  ,  This  is  a  very  interesting 

c  c  cn  i  ii  7  j 

observation,  for  (i c+ i is  the  current  which  would  he  n  easured  by  physically  combining 
the  plane  and  [mint.  In  this  case,  however,  the  electrical  interconnection  of  these 
electrodes  is  effected  via  the  signal  resistors  and  an  earth  point.  They  are  electrically 
decoupled;  the  coupling  occurs  via  the  physical  processes  in  the  gap. 

3.2  Three  Electrode  System 


3.2.1  Three  Electrode  System 

Tig  3. ft  is  an  exaggerated  sketch  of  the  electrodes  used  in  this  study.  To  maintain 
consistency  with  what  follows,  the  point  will  he  referred  tn  as  S  ,  the  plane  anode  ns  5^ 
and  the  cathode  as  It  v/ill  be  shown  that  the  needle,  5^,  introduces  an  effect  ive 

boundary  condition  and  to  emphasize  the  generality  of  this,  a  purely  geometric  function, 
f(r),  v/ill  be  used. 

A  one  dimensional  filament  is  considered.  Two  coordinate  fran  es  are  used,  one, 
x,  having  its  origin  at  the  anode,  5^,  the  other,  r,  having  its  origin  at  the  ent!  of  the 
needle,  Sc,  where  x=n,  the  height  of  the  needle  above  the  anode  (Fig3.ftl. 
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T  :  ~  ; r  ■  i ' i ,  i  1 1 1 1 r  l . r r i ■ '  •  f 

Tn  find  r)(_  is  at  unit ,  and  S(i,  r,  at  zero  potential  Tcctiun  3.1.1'. 

A ceept  mg  for  new  th.it  tlm  per-unit  pntcntial  V  «,  only  has  an  appreciable  value  near  the 
needle  end,  r  <•  art  where  ar^  is  a  multiple  of  the  tip  radius!  ]this  v.ould  seem  reasonable 
if  the  needle  tip  were  an  isolr  ted  sphere  uf  radius  r(,  then  if  the  surface  is  at  ur.it 
potential,  beyond  the  surface  the  potential  will  tie  r^/r.  The  earthed  anode  will 
exaggerate  this  fall-eft.0),  the  general  form  expected  of  V'c*  is  given  by  the 
dimensionless  function  ffr)  .Tig. 3. 4'  which  has  the  following  properties: 

f(r=0)  x  !;  ffrrd-n;  =  0;  f'rmi-n)  =  0  =  f(n'(r  =  ri-n): 
f(r)-v0  unleso  r<  ar^(3.h) 


From  Fq.  3.3,  with  V  «  =  f(r), 

r-d-n, 

\  3  J Ecj dr  -  -J 


(3.7; 


r.  o 

and  is  hence  the  average  current  flowing  in  the  portion  of  the  channel  where  f(r)  is 
appreciable,  so  given  that  the  spatial  extent  of  f(r)  is  small,  i  will  be  the  same  as  the 
conduction  current  in  the  channel  base. 

The  per-unit  field  is  confined  to  a  small  region  nrnr  the  tip  and  a  charge  travailin' 
across  this  region  will  pass  through  a  difference  Av*  in  V*  of'v).  So  during  this  transit, 
the  whole  charge  q  is  transported  to  the  point.  If  the  size  of  this  region  is  so  small  that 
the  electrons  cross  it  quickly  compared  with  their  arrival  rate,  or  compared  with  the 
temporal  resolution  of  the  oscilloscope,  then  the  current  in  this  electrode  approximates 
to  the  arrival  rate  of  the  charge. 


3.2.3  Displacement  Current 

* 

V  ^  for  off-axis  charge  trajectories,  is  exactly  the  same  as  for  plane  parallel 

electroncs,  and  i  ,  =  i  (section  3.1.3),  but  because  the  needle,  5  ,  is  earthed  when  finding 
'  (!  p  C 
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V  (rather  than  linearly  rising  as  a  charge  across  the  gap,  as  in  section  3.1.3,  r ig.3.3' 
along  the  gap  axis  q^must  fall  abruptly  to  zero  as  r  *♦()  (Fig. 3. 5'.  The  charge  induced  nn 
the  plane  will  rise  'inearly  as  q  moves  across  the  gap,  only  to  fall  abruptly  hack  to  zero 
as  the  needle-end  is  approached.  The  charge  built  up  on  the  plane  as  q  traversed  the  gap 
is  transferred  rapidly  to  the  point.  (This  is  a  slight  simplification:  the  charge  built-up  on 
the  plane  is  (d-n)q/d,  and  that  remaining  on  the  cathode  is  nq/d  both  of  these  are 
transferred  to  the  point  as  the  test  charge  traverses  the  high  field  region  at  the  needle- 
end.  Gf  course  the  charge  induced  on  the  plane  is  a  function  of  the  position  of  q,  not  its 
path.  The  discussion  is  couched  in  terms  of  the  trajectory  of  q  to  make  visualization  of 
the  spatial  variations  clearer.}  From  the  point  of  view  of  the  conservation  of  induced 
charge,  (Eq.3.9),  it  is  self-evident  that  the  only  possible  solution  for  V  ^  is: 


Since  it  is  contributed  to  by  the  displacement  of  charge  across  the  gap,  this  is 
called  the  displacement  current.  Rearranging  Eq.3.10  we  have  ie>  the  current  which 
would  have  been  recorded  in  plane  parallel  geometry: 


ic  *  an: 

This  addition  of  a  negative  contribution  to  from  i  formalizes  the  rapid  drop  in 
V*d  near  the  needle  end,  and  will  turn  out  to  be  a  useful  analytical  tool.  This  rapid 
drop-off  in  V*d  is  due  to  an  effective  boundary  condition,  since  the  defining  equations 
for  V  d  require  the  needle  to  be  at  earth  potential.  This  is  where  the  i  term  in  Cq.3.10 
comes  from  and  is  a  fundamental  result. 


Illustrating  the  Shockley-Ramo  theorem:  Point  current  ip 
and  center  probe  current  lc  induced  by  one  O2  ion  drifting 
from  the  hyperboloid  point  to  the  plane  through  the  Laplacian 
field  in  87.6  kPa  dry  air  (Sigmond,  1978). 


Idealised  ^  -  id  current  profiles 
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This  technique  for  recovering  i  is  not  dependent  on  the  extent  of  f (r).  f.cr  is 
Eg.?. 7  for  ic.  The  only  place  in  this  discussion  where  the  spatial  extent  of  f(r)  is 
important  is  in  ascribing  i  to  he  the  physical  current  in  the  channel  hase.  This  will  be 
discussed  later  in  connection  with  the  space-charge  field  at  the  channel  base. 

3.2.A  Free  Charge  Reconstruction 

As  already  mentioned,  the  electron  drift  speed  is  »  ion  drift  speeds  and  only  the 
electronic  component  of  the  current  is  observed.  Mote  that  this  implies  that  free 
(electron}  ~harge  is  NOT  conserved;  an  electron  attaching  to  an  oxygen  molecule,  for 
example,  contributes  little  to  the  induced  current. 

If  the  free  electron  charge  in  the  gap  totals  G^,  and  IF  these  electrons  have  the 
same  drift  speed  v^  (the  field  is  uniform,  and  vr  is  found  from  well-established 
expressions)  then  ig  =  Q^v^/d  or: 


(3.12) 

and  the  free  electron  charge  in  the  gap,  G  can  be  found  by  simply  scaling  i  .  This 

9  e 

reconstruction  of  G^  from  instantaneous  value  of  the  induced  currents  relies  on  some 
fairly  gross  assumptions,  and  the  validity  of  these  is  considered  below  (section  3.3). 

3.2.5  Three  Electrode  Point-Plane  System 
Sigmond  (197C)  has  calculated  the  per-unit  potentials  and  fields  for  an  ingenious 
three  electrode  point-plane  system  used  by  Marodc  (1975a, b)  (Fig. 3.6).  The  cathode  is 
here  divided  into  two,  with  a  small,  isolated  central  portion.  Note  how  the  drop-off  of 
V  *  is  so  steep  that  the  charge  is  transported  almost  instantaneously  to  the  point  as  q 
traverses  the  high  field  region.  To  a  lesser  extent  this  occurs  as  q  reaches  the  central 
portion  of  the  cathode  too,  and  hence  the  currents  n. easureri  will  be  nearly  equal  to  the 
charge  arrival  rates,  as  outlined  above.  Thus  for  this  arrangement,  as  expected  for  the 
system  used  here,  the  current  to  the  point  is  effectively  equal  to  the  arrival  rate  of 
charge. 
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3.2 Polarity  of  Induced  Currents 

The  applied  field  is  supplied  by  Scath,  and  charges  of  either  polarity  drift  in 
opposite  directions  (according  to  the  coordinate  system  of  Fig.3.4)  to  produce  a  *ve 
conventional  current  j  ( directed  upward  in  the  figure).  Frnm  Fg.3.7  and  Cq.3.5  is 
-vey,  it  can  be  seen  that  i ^  and  ip  will  always  be  *ve,  and  hence  will  always  cause  a  -ve 
voltage  across  the  signal  resistor  R  in  the  external  circuit. 

All  references  to  currents  induced  from  now  on  will  refer  to  the  voltages  across 
the  resistors  as  though  they  were  currents.  So  a  -ve  conduction  current  waveform  (for 
example)  is  produced  by  an  upward  directed  conventional  current  in  the  base  of  the 
channel. 

The  only  *ve  voltage  (-ve  current)  induced  in  the  external  circuit  to  the  lower 
electrodes  i£  the  (-(d-n)/d,i  contribution  to  i  .  If  there  were  o  positive  voltarp  in  i 
therefore,  it  can  be  inferred  that  there  is  current  flowing  in  the  region  where  the  per- 
unit-potentiol  of  the  needle  is  appreciable.  This  is  a  result  of  there  being  an  effective 
boundary  condition  at  the  needle-end  and  is  unambiguous.  The  only  ambiguity  in  the 
analysis  lies  in  the  spatial  extent  of  the  needle-end  region.  In  the  sense  indicated  above, 
positive  electron  currents  ie,  and  conduction  currents,  i^  are  physically  not  possible. 

3.3  Limitations/Problems 

The  above  treatment  of  the  induced  currents  relies  on  some  fairly  gross 
assumptions;  however  it  will  now  be  shown  that  these  can  largely  be  avoided  v/hen  the 
particular  magnitudes  of  the  parameters  involved  are  considered  more  carefully. 


3.3.1  Space-Charge  Field 

flection  3.2 A  ignored  the  space-charge  field.  Firstly,  electron  velocities  in  the 
active  region  arc  ‘'■'107  cm/s.  The  contribution  from  the  space-charge  field  of  the  tip, 
where  the  electron  speed  is  the  highest  is  considered  belniw  'section  3.3.2).  Tven 
excepting  the  tip,  for  G^  to  be  recoverable  from  i^,  the  electron  speed,  and  therefore  the 
local  field  everywhere  must  he  the  san  e.  But  is  this  a  reasonable  assumption-?  Can  the 
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space-charge  field  of  the  channel  at  least  be  ignored'* 

Phelps  (l^Vt'O  asserts  that  the  field  along  the  channel  is  just  the  goon  etric,  or 
applied  field;  whereas  N'arodc  (1901)  emphasizes  that  conductivity  is  the  more  important, 
T  being  determined  such  that  there  is  current  continuity  along  trie  channel. 

Phelps  would  have  the  space-charge  distributed  by  the  applied  field:  N'arode,  the 
conductivity  (charge  density)  determining  the  field. 

A  simple  calculation  shows  that  for  a  channel  of  net  charge  ^10nC  (typical 
value  of Ti^dt),  the  space-charge  field,  KV/cm,  compared  with  nn  ambient 

field  of  some  4KV/cm.  Hence  the  space  charge  field  con  safely  be  ignored  in  the 
channel. 

There  is  an  exception.  It  will  be  recalled  that  the  development  of  further 
streamers  is  supposed  to  be  suppressed  by  a  positive  space-charge  local  to  the  needle- 
tip.  The  suppression  of  further  streamers  is  indeed  good  evidence  that  the  field  at  t fie 
tip  is  reduced  by  a  positive  space-charge.  Certainly  the  geometric  field  here  is  high, 
causing  electrons  to  be  swept  out  of  the  gap  more  efficiently  than  they  ore  swept  into 
this  region  from  the  rest  of  the  ch  nncl.  It  will  tend  to  accumulate  a  positive  charge 
until  the  local  field  is  uniform  along  the  channel.  But  this  will  occur  only  whilst  the 
electrons  in  the  needle-region  remain  unattached.  This  is  essentially  either  the 
development  of  compensation  zone  process  of  Phelps  or  the  secondary  streamer  process 
of  Marode. 

In  this  case  the  space-charge  Held  works  in  favour  of  the  theory  developed  here. 
V/ith  the  reduced  local  field  due  to  this  space- charge,  the  electron  drift  speed  will  be  the 
same  in  the  base  of  the  channel  as  it  is  in  the  remainder,  and  the  recovery  of  the  free 
charge  (section  3.2.4)  thus  makes  a  valid  assumption.  The  electron  drift  speed  can  be 
obtained  accurately  enough  by  using  the  applied  field  and  pressure  for  a  given  gas  -  it 
does  not  change  significantly  because  of  the  space-charge. 

Furthermore,  the  superposition  principle  permits  the  effect  of  individual  charges 
to  be  considered  in  isolation  and  simply  summed  up  (section  3.1.2),  so  the  space-charge 
does  not  influence  the  per-unit-potential;  V*  is  a  property  of  the  electrodes. 
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It  tins  been  indicated  (section  3.2.2',  that  the  time  laker,  to  cross  the  needle-region 
is  important.  If  this  occurs  within  the  sample  tin  e  of  the  oscilloscope,  the  charge 
transfer  can  be  considered  to  have  occurred  instantaneously.  Tfie  slinrtest  sample 
interval  used  was  ]f>ns  and  for  an  electron  drift  speed  of  -viCkm/s  'in  the  whole  length  of 
the  channel',  this  corresponds  to  a  needle-end  region  64U|im  in  extent.  Since  the  tip 
radius  of  the  needle  used  was ~20pm,  the  conduction  current  can  safely  lie  identified 
with  the  arrival  rate  of  charge  at  the  tip.  The  transit  time  of  ttie  enhanced  per-unit- 
field  region  is  clearly  less  than  it  v/ould  have  been  without  the  positive  space-charge,  but 
it  has  been  shewn  that  the  overall  effect  of  this  does  riot  affect  the  identification  of  the 
conduction  current  with  the  arrival  rate  of  charge  at  the  tip. 

3.3.2  Steamer  Active  Region 

A  major  problem  in  the  analysis  of  the  discharge  currents  is  the  different  drift 
speeds  in  the  active  region  and  in  the  channel.  The  field  in  the  active  region  (Hartmann, 
1974}  may  be  1 2f)V/m>.torr  or  more.  The  corresponding  electron  speed  is  wv330'xm/s,  or 
20  times  greater  than  thc~30km/s  in  the  drift  region  (the  channel}.  However,  typical 
free  charge  values  before  cathode  impact  may  be  up  tn  lOOnC,  1CJ  x  the  charge  in  a 
streamer  tip  (10*^  ^C}  and  thus  for  a  single  tip  the  channel  current  would  be  ~]0^x 
greater  than  the  tip  current. 


3.4  Application 

Turning  now  from  generalities  to  some  specific  predictions  about  the  inducer.' 
currents,  etc.  expected  from  particular  discharges,  there  are  three  cases  to  consider. 

The  current  may  be  confined  entirely  to  the  tip  region  (i)  or  distant  from  the 
tip  region  (ii),  or  a  combination  of  both  (iii). 


3.4.1  Case'i} 


If  the  current  were  confined  to  the  tip  region,  the  electrons  would  drift  to  tbp 
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anode  surface  very  quickly  -me  be  neutralized  there.  Thus  only  a  brief  pulr.i-  of  current 
would  hr  observed.  As  discussed  above,  there  would  be  a  large  contribution  to  the 
conduction  current,  because  of  the  relatively  large  per-unit  field,  accompanied  by  a  large 
positive  contribution  to  the  displacement  current.  This  positive-going  displacement 
current  is  the  clearest  evidence  for  this  phenomenon,  as  it  can  arise  from  no  other 
situation  (negative  currents  are  always  expected).  This  kind  of  event  would  look  much  os 
sketched  in  Fig.  3.8(1). 


3,4.2  Ccse  (ii) 

If  the  charge  is  distant  from  the  point,  then  no  contribution  would  be  expected  to 
the  conduction  current;  the  per-unit  field  witn  respect  to  the  tip  is  negligible  at  any 
significant  distance  from  the  needle  end.  The  most  practical  visualization  of  this  kind  of 
behaviour  is  charge  injection  from  the  cathode  arising  from  the  impact  of  the  streamer. 
In  this  case  there  would  be  a  continuum  of  current  from  whatever  processes  are  occur¬ 
ring  in  the  remainder  of  the  discharge  immediately  before  the  impact,  and  the  burst  of 
charge  released  from  the  cathode  will  be  superin'. posed  upon  this.  The  electrons  liber¬ 
ated  from  the  cathode  v/ill  then  attach  over  a  period,  producing  the  waveforms  in  Fig. 
3.8(H). 

As  the  electrons  from  the  cathode  are  distant  from  the  point,  there  will  be  no 
induced  conduction  current.  The  electron  current  i  will  be  very  much  like  the  dis¬ 
placement  current  and  both  will  be  negative.  The  situation  is  similar  to  plane-plane 
geometry,  for  which  the  displocen.ent  current  would  be  identically  equal  to  the  electron 
current.  An  electron  generating  event  in  mid-gap  would  have  produced  the  same  effect. 

3.4,4  Case  (iii) 

V.'here  there  is  a  channel  growing  into  the  gap  from  the  point,  the  situation  is 
more  complex.  The  relative  growth  and  attachment  rates  need  to  be  considered.  There 
are  three  main  cases,  (a)  below,  (b)  at,  and  (c)  well  above  the  stability  field.  Based  on  a 
simple  model  of  the  discharge  'section  3.5.6  ff)  the  electron  currents  are  indicated  in 
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Fig. 3.9  for  each  case.  As  the  channel  is  built  up  the  electron  current  increases  initially, 
falling  again  for  (a),  flattening  out  for  (b)  and  subsegucntly  rising  in  case  (c).  Thr  in¬ 
fluence  of  this  on  the  displacement  current  can  be  calculated  with  the  aid  of  Fq.  3.!fl  and 
id  is  sketched  for  (a),  (b)  and  (c),  assuming  the  cose  of  a  standard  conduction  current 
pulse  in  Fig.  3.9. 

It  can  be  seen  that  the  dying  streamer  has  a  displacement  current  wholly  positive, 
only  marginally  different  from  the  conduction  current.  The  stability  streamer  dis¬ 
placement  current  does  dip  below  the  avis,  and  a  little  charge  injection  from  the  cathode 
has  been  marked  in  (according  to  section  3.4.2).  Above  the  stability  fiel-d  the 
displacement  current  dips  very  strongly  below  the  avis,  and  the  charge  injection  has  been 
drawn  in  n  ore  heavily  also.  In  this  example  the  second  derivative  of  i^  has  changed  from 
a  concave  slope  below  E  t  ^  to  a  convex  slope  above  ^  If  the  attachment  is  very 
strong  (Fig. 3.9(d)),  then  a  point  of  inflex  ion  is  expected  to  appear. 


3.5  Further  Anal  ysis 


3.5.)  Generated  Charge 


Consider  a  channel,  length  Lj  and  charge  Gj  at  time  tj,  developing  to  (l.^,  C^, 

/’t» 

t~).  Then  since  I  i  dt  is  removed,  and  the  remainder  gets  attached,  one  may  write: 

Jt  ix 

>  v  «  /  /  «  «  N 


Q,  •  (c? 


where  the  generated  charge  £C  is  introduced  to  satisfy  charge  conservation.  In  a 
grov/ing  streamer  system,  this  is  in  effect  the  total  streamer  tip  current  integrated  over 
a  sample  interval.  In  all  calculations,  it  will  be  divided  by  the  sample  interval,  to  give  a 


charge  injection,  or  charge  generation  rate.  It  will  be-^O  after  the  streamer  transit, 
unless  sane  other  electron-detaching  or  ionizing  process  is  occurring. 


3.5.2  Attachment  Frequency 

To  find  from  Eg. 3.14  requires  a  knowledge  ofv  and  this  is  not  well  defined 

£3 

experimentally.  However,  if  the  channel  is  decaying  without  any  new  electron 
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detachment  or  ionization,  =  0.  inmost  cases,  the  channel  lias  decayed  substantially 
by  the  time  the  streamers  cross  the  gap  anclji^dt  in  Fq.  3.  !4  can  be  neglected  b>  coni  - 
parison  with  the  magnitude  of  G^.  Typical  values  are  for  to  change  by  2-3nF  at  an 
attachment  frequency  of  -'-3MHz  compared  with -»-320pC  extracted  fron.n  the  channel 
base  in  the  same  time.  So  unless  the  attachment  frequency  is  very  low  and  there  is 
appreciable  conduction  current  at  this  point,  can  be  found  by  performing  a  least- 
squares  exponential  fit  on  the  data.  An  example  of  this  is  shown  in  Fig3.ID  where  this 
has  been  performed  between  the  two  cursor  markers.  In  this  way,  using  the  fitted  value 
of  Vg  from  a  portion  of  the  discharge  when  =  0,  allows  to  be  found  for  portions  of 
the  curve  when  it  is  non-zero  (Fig.3.11). 


3.3.3  Streamer  Velocity 

The  large  pulse  of  electrons  released  when  the  streamer  arrives  at  the  cathode,  is 
seen  as  a  sharp  increase  in <£o  (Pi g. 3- J J ).  From  the  transit  time  of  the  streamer,  t^,  the 
mean  streamer  speed,  vg  =  (d-n)/t^  can  be  found.  This  time-of-flight  method  is  rather 
crude  compared  v/ith  considerably  more  sophisticated  photomultiplier  studies,  but  the 
results  are  useful  as  confirmation  or  for  comparison. 


3.5.4.  Streamer  Function,  A(x) 

In  attempting  to  identify  a  formal  growth  parameter  for  a  streamer  consider  that 
a  streamer  system  "creates"  an  elementary  line  charge  dQ  =  A(x)dx,  at  x,  at  time  t=x/vs 
(Fig. 3. 12):  this  charge  then  drifts  at  v^  back  toward  the  point,  arriving  x/v  seconds  later 
at  time  T  =  x/vs  +  x/ve  -  x/V  where  it  is  convenient  to  define  an  effective  velocity  V. 
After  attachment,  the  charge  dG  laid  down  becomes  upon  arriving  at  the  point,  dG'  = 
riGe"7x  =  A(x)e'1xdx.  Hut  with  x=VT,  dx=VdT,  and  dG'  =  A(TV)e“ 7"^VdT,  then: 


i  -  «  AfTvW^V 

c  di  1  ; 

A60-  it^je^/v 


or 


(3.16) 
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dQ'  =A(*)e-?*d* 

=  A(TV)e'  VdT 


dQ  =  A(x)ds 

(created  at  time  t  =  «/»,  ) 


Fig  3.12 


ic(T)  =  dQ'/dt  =  A(TV)e-  ,?TVV 


charge 
dQ'  *oe- 

81  =  o  at  time  T 

’r  Vve  =  */V 


1  =  M  ♦  ‘2 

d/*s  =  X/*»  +  X/»e  -  */*e 


X  4 


streamer  reaches  d  *8t 

charge  deposited  at  X  is  A(X)dX 
at  time  tj  =  X/j;, 


=>  X  =  (d/»s  *  x/spJV 
>(*,0  =  A(X)e' 


drifts  to  x  in 
time  t2  =  (X  *)/*e 
and  is  A(X)e  7*X  **dX 


Fig  3.13 


31 


Calculation  of  A'x'  from  i  is  only  possible  while  significant  current  survives  to 
re.ii'b  the  point.  Noise  otherwise  dominates  the  calculation  of  A'x',  and  the  limiter' 
information  yielded  by  this  technique  limits  its  use.  This  demonstrates  ttip  rinmnanre  of 
attachment  in  the  channel.  Compared  with  the  attachment  effect,  A'x'  tins  negligible 
influence  upon  tt>e  form  of  it. 


3.5. f>  Channel  Function,  A^x' 

The  instantaneous  electron  distribution  along  the  channel,  ^(x'i,  may  be  more  useful 
t  Nor  A(x'.  Clearly  C^=Ja(xxi'x,  and  A'  x)  is  related  to  A.'x)  by  the  following  relation,  'see 

-l(X-x) 


fig.  3.  13’ 


A^x,t)  =  A(x)e 


X’  f^+ir<)V’  13.17: 


In  this  case,  t:(f,'ic,  Substituting  this  into  fg.3.17  and  integrating  the  electron 

.  -,A; 

i  i large  for  a  streamer  which  has  extended  a  distance  o',  and  putting  A  X,  -J\.o  for  tf  e 
st man  er  function,  one  has: 

..'if'  =  J^\(x)c/x 

a 

o 


■3. !  P.' 


This  has  been  evaluated  with  the  following  values  of  the  coefficients  involved: 

-  If •(»".*,  ve=<*niWs,  vs  =  lrsOUm/s  ....  V=3i.5kin/s).  Ttic  result  is  plotter'  us  fur 

various  values  of  d  and3"in  Fig.3.ia  and  shows  that  this  general  picture  of  the  discharge 
is  reasonable  insofar  as  it  produces  results  which  agree  with  the  phenomena  and  v. it.’'  the 
general  idea  of  a  growing  streamer. 

If  the  streamer  is  dying,  2f<D,  then  Op  decreases  with  d.  G^  is  stable,  reaching  a 
constant  level  whence,  and  rises  rapidly  when  there  is  strong  growth, JT>0.  Initially,  as 
the  channel  is  being  established,  all  values  ofij-  imply  growth  in  G  The  stability  field 
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Section  4:  Typical  Results 


Introduction 

The  input  data  consists  of  just  two  time  dependent  currents  (viz.  conduction  and 
displacement)  but,  as  indicated  in  Section},  a  number  of  derived  quantities  are  obtained 
for  each  corona  event  (e.g.  free  electron  magnitude  Qy  generated  charge  S  Q). 
Naturally,  the  temporal  development  of  these  quantities  is  field  and  pressure  depen¬ 
dent.  In  this  section  some  typical  results  are  presented  for  various  field/pressure 
regimes  together  with  an  explanation  of  these  profiles  in  terms  of  the  microphysics  of 
the  event. 

4.1  Regular  Corona 

The  bulk  of  the  work  concerns  conventional  streamers  at  sub-breakdown  fields  and 
at  different  pressures  and/or  propagation  lengths.  Subsequently,  as  the  field  is  increased, 
a  number  of  corona  modes  are  observed,  which  are  associated  with  pre-breakdown  pheno¬ 
mena.  These  corona  modes  are  interesting  as  their  elucidation  bears  upon  the  topic  of 
the  glow-to-arc  transition  in  uniform  fields  which  is  commonly  encountered  in  overvolted 
gaps,  at  a  reduced  field  nearly  always  greater  than  34V/cm.torr  (for  air);  by  comparison, 
the  breakdown  in  the  present  study  occurs  at  8  to  lOV/an.torr.  We  deal  here  with  the 
sub-breakdown  streamer  event. 


4.1.1  Inception  field 

At  reduced  pressure,  the  lateral  diffusion  of  electron  avalanches  increases.  To 
generate  a  critical  avalanche  then,  a  greater  charge  is  required  in  this  larger  volume  to 
raise  the  tip  field  to  the  level  where  the  replication  mechanism  can  operate.  So  a  higher 
reduced  field  is  needed  before  onset  occurs.  Replotting  Phelps'  stability  field  data  as 
reduced  stability  field,  Estab^’  s^ows  on  the  contrary  that  the  stability  field  is  lower  at 
lower  pressure  (Fig.4.1).  The  combined  effect  of  these  two  variations  is  that  at  a  certain 
pressure  (about  ~500  torr  in  the  present  study),  the  onset  and  stability  fields  coincide. 
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Below  this  value,  onset  streamers  always  grow  throughout  their  life,  and  regrettably  no 
sub-stability  field  measurements  can  be  made. 

4.1.2  Pre-stability  Field  Streamers 

As  mentioned  above,  observations  of  this  regime  are  restricted  to  high  pressures. 
The  displacement  current  at  low  reduced  field  is  entirely  positive  (Fig.  4.2,4408.2) 
indicating  that  there  is  a  lot  of  current  at  the  base,  and  very  little  in  the  body  of  the 
channel.  (It  will  be  recalled  that  a  positive  displacement  current  is  fundamentally 
related  to  the  current  at  the  needle-tip,  section  3.2.6).  The  conduction  current  shows  no 
fine  structure  (more  of  this  later)  and  the  enall  electron  current  (Fig.  4.3)  ensures  that  ic 
and  ij  are  almost  anti-symmetric  (Fig.  4.4, #408.3). 

The  electron  charge,  Q  ,  falls  off  throughout  the  discharge  (Q  and  i  have  the 
9  9  e 

same  form,  as  they  are  related  by  a  simple  scaling  factor  (section  3.2.4)  and  all  the 
figures  are  drawn  with  a  dual  scale  representing  this).  The  form  of  in  fact  follows 
that  presented  in  Fig.3.14  for  the  case  of  a  strongly  decaying  streamer  system.  Common 
to  all  records  is  the  initial  rise  in  Q^.  Following  this  the  particular  features  of  interest 
are  the  peak  -*0.23ps  after  the  event  began  and  following  that  the  subsequent  strong 
decay  of  Q^.  The  simple  calculation  of  the  charge  in  the  gap  (section  3.5.6)  shows  a 
similarly  broad  peak,  which  indeed  (for  example  with  the  growth  parameter  used  in  the 
modeiyr-10),  occurs  at  around  -»0.27|is.  All  this  shows  clearly  that  the  low  reduced  field 
in  the  example  studied  cannot  support  streamer  propagation,  even  though  a  critical 
avalanche  can  be  created  at  the  needle  Up.  Once  the  high  field  region  at  the  needle  end 
has  been  left  behind,  the  streamer  travels  a  cm  or  so  using  its  internal  energy  to 
compensate  for  the  inadequate  supply  from  the  external  field  and  finally  dies.  Visual 
observation  of  the  corona  confirmed  that  it  did  not  propagate  to  the  cathode. 

It  has  been  observed  though  that  streamers  below  the  stability  field  may  still 

» 

survive  to  the  cathode.  Figure  4.5  shows  a  sequence  of  records,  as  E/p  is  increased.  The 
event  just  considered,  #408.2,  was  obtained  at  the  very  low  ambient  reduced  field  of 
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4.24V/ cm. torr.  At  a  slightly  higher  reduced  field,  4.5V/cm.torr,  the  streamers  reach  the 
cathode,  and  a  small  current  pulse  is  seen  demonstrating  this  (Fig.4.5ii).  This  record 
though  (4409.3)  shows  that  Qg  was  decreasing  throughout  the  streamer  propagation. 

Thus  although  it  reached  the  cathode,  this  streamer  was  not  growing,  and  according  to 
section  3.4  was  below  the  stability  field.  The  decay  is  more  clearly  exhibited  in 
Fig.4.5iii,  showing  Qg  for  event  #97.0.  This  was  for  a  longer  propagation  distance 
(16.4cm)  than  for  the  previous  example  (10.3cm),  and  so  the  reduction  in  Qg  has  a  longer 
time  to  manifest  itself.  The  reduced  field  of  4.67V/cm.torr  for  event£97.0  is  actually 
higher  than  the  4.5V/ cm  .torr  used  in  the  previous  example  (#409.3)  though  the  decay 
seems  stronger:  this  is  of  course  partly  caused  by  the  relative  timescales  of  the  plots  and 
partly  by  the  longer  total  discharge  time  in  the  latter  event.  Fig.4.5iv  shows  an  event  (# 
99.1)  where  the  streamers  grew  slightly  throughout  the  discharge;  Q  increases  con¬ 
tinually  up  to  cathode  arrival.  This  trend  is  confirmed  by  Figs.4.5v  and  4.5vi  where  the 
growth  in  Q^  increases  with  the  reduced  field. 

4.1,3  Streamers  at  the  Stability  Field  or  Greater 
Streamers  can  cross  a  gap  at  fields  below  the  stability  field.  There  is  a  range  of 
fields  at  which  a  proportion  of  the  streamers  cross  the  gap  and  a  proportion  do  not.  This 
can  be  clearly  seen  from  the  displacement  current  records,  which  either  do  or  do  not 
contain  a  sharp  negative  feature  as  the  streamers  resch  the  cathode  (or  not).  This 
feature  is  considered  here. 


4.1.3a  Injection  from  Cathode 

Fig.4.6  shows  a  sequence  of  streamer  displacement  and  conduction  current  records 
for  an  increasing  reduced  field.  They  all  show  a  sharp  negative  feature  in  the  dis¬ 
placement  current,  which  is  accompanied,  if  at  all,  by  only  a  relatively  insigificant 
feature  in  the  conduction  current.  The  peak  in  the  displacement  current  is  evidence  that 
there  has  been  a  sharp  increase  in  the  free  charge  in  the  gap,  for  ie  must  have  in- 
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creased.  The  lack  of  a  correspondingly  large  signal  simultaneously  in  ic  indicates  that 
the  charge  generation  responsible  for  this  feature  occurs  well  away  from  the  point 
(section  3.4).  It  has  been  assumed,  since  visual  observation  confirms  that  the  streamers 
cross  the  gap  when  this  negative  feature  in  i^  is  present,  that  this  generation  of  charge 
occurs  at  the  cathode  when  the  streamers  arrive  and  the  cathode  region  is  established. 

This  rapid  injection  of  charge  is  siperimposed  upon  the  now  growing  for  the 
system  as  a  whole.  The  two  are  partially  separable  by  calculating  Sq  as  described  in 
section  3.5.1.  Fig.4.5  shows  the  8 Q  records  corresponding  to  the  records  already 
presented  for  various  reduced  fields.  The  impulsive  nature  of  the  charge  injection  can  be 
seen  more  clearly  when  plotted  as  SQ  demonstrating  that  this  feature  increases  in  size 
with  C/p;  this  can  be  linked  with  the  ultimate  size  of  the  streamer  system.  It  would  not 
seem  unreasonable  that  the  charge  injection  from  the  cathode  should  increase  as  the 
number  and  size  of  the  streamers  arriving  there  increases. 

4.1.3b  Conduction  Current  Perturbation 

Simultaneously  with  the  large  spike  in  the  displacement  current  record  a  small 
perturbation  (Fig.4.7)  is  sometimes  seen  in  the  conduction  current  record.  Quite  what 
this  represents  is  an  open  question,  but  it  is  certainly  much  smaller  than  it  would  be  if 
the  current  pulse  responsible  for  the  displacement  current  feature  were  located  at  the 
needle-end  and  must  have  some  other  explanation. 

One  possibility  is  that  the  theoretical  analysis  presented  in  section  3  is  not 
perfectly  followed  by  the  practical,  real,  electrode  arrangement.  Some  induced  current 
could  be  being  picked  up  by  the  needle  from  charges  located  at  the  cathode;  for  example, 
if  E  *  were  not  zero  everywhere  far  from  the  needle.  However,  even  if  there  were  any 
such  contribution,  it  would  at  least  be  continuous  with  the  position  of  the  space-charge, 
and  not  vary  rapidly  far  from  the  point.  Therefore  the  perturbation  would  have  the  same 
shape  as  the  displacement  current  record,  which  in  fact  tails  off  by  far  the  more  slowly 
of  the  two.  This  feature,  furthermore,  would  always  be  present  whenever  the  streamers 
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cross  the  gap,  although  in  fact  it  is  not  always  seen,  ruling  out  this  particular  possibility. 

The  alternative  is  that  an  enhanced  field,  associated  with  a  potential  redistri¬ 
bution  brought  about  by  a  return  wave,  causes  the  electrons  in  the  channel  base  to  speed 
up,  resulting  in  a  slightly  increased  conduction  current.  Though  since  the  perturbation 
occurs  simultaneously  with  the  peak  in  i^  (or  rather  within  the  16ns  sample  window),  this 
implies  a  propagation  velocity  in  excess  of  llNtn/s,  compared  with  a  value  of  —TONtn/s 
reported  by  Marode  (1975a)  for  the  speed  of  the  wave  during  the  compensation  phase. 
This  velocity  is  —100  times  that  observed  by  Suzuki  )1971)  so  it  seems  there  is  some  con¬ 
siderable  disagreement  regarding  the  speed  of  ionizing  waves.  The  wave  observed  by 
Marode  occurred  under  the  same  circumstances  as  the  perturbation  discussed  here,  and 
so  is  perhaps  the  better  value  for  comparison.  Thus  it  seems  reasonable  to  identify  this 
perturbation  with  the  arrival  of  a  return  wave  -  as  concluded  by  Marode. 

4.1,4  Attachment  frequency 

Once  an  entire  channel  has  been  created,  no  further  increase  in  electron  charge 
Qg  is  expected  (unless  of  course  some  other  electron  detaching  process  occurs).  Fitting 
an  exponential  curve  to  the  later  portion  of  the  Qg  record  gives  very  close  agreement, 
showing  that  the  charge  residing  on  the  channel  is  progressively  attached.  Attempts  to 
fit  exponentials  to  streamers  which  do  not  cross  the  gap  must  be  treated  with  caution, 
for  the  channel  may  still  be  extending  although  electron  density  is  declining  (see  e.g.  the 
model  in  section  3.5.6).  The  result  (e.g.  Fig.4.5i  shows  the  result  of  such  a  calculation) 
may  indeed  be  exponential,  but  does  it  represent  only  attachment?  The  decay  of  can 
only  be  considered  to  be  attachment  after  cathode  impact,  so  (e.g.)  in  Fig.4.5ii  the 
calculation  of  the  attachment  frequency  based  on  a  fit  between  the  two  cursor  positions 
does  satisfy  the  assumption  that  fc=0.  This  topic  will  be  discussed  in  the  next  section. 
The  attachment  algorithm  is  described  in  the  appendix. 
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4.1.5  Low  Pressure  Specifically 

The  results  at  low  pressure  are  very  interesting.  Clearly  they  are  of  some 
importance  for  the  questions  this  work  addresses,  for  a  thundercloud  may  exist  in  a  low 
pressure  regime. 

Fig.4.8  shows  a  low  pressure  streamer  record.  It  is  completely  different  fro  those 
at  high  pressure,  and  one  could  be  forgiven  for  thinking  it  to  be  some  new  kind  of  event 
and  not  a  streamer  at  all. 

■^/p  for  two-body  attachment  and  l/p*  for  three-body  attachment  depend  only  on 
E/p.  As  the  pressure  is  reduced,  the  attachment  coefficient  is  predicted  to  fall,  and  this 
is  indeed  observed  to  happen.  This  appears  to  be  the  principal  effect  on  the  discharge  of 
reducing  pressure  and  is  expressed  in  a  number  of  ways. 

4.1.5a  Induced  Current 

The  conduction  current  is  due  to  the  electrons  produced  in  the  streamer  active 
region  drifting  along  the  channel  and  arriving  at  the  point.  The  magnitude  of  the  con¬ 
duction  current  will  thus  depend  on  the  number  of  electrons  "added"  to  the  channel,  A(x) 
(section  3.5.5),  and  upon  the  attachment  coefficient.  In  this  case,  the  reduced 
attachment  allows  more  electrons  to  survive  to  the  point,  enhancing  the  conduction 
current.  Now  it  is  inevitably  true  that  the  electron  current  will  also  be  enhanced  by  re¬ 
duced  attachment,  but  for  a  charge  at  the  base  of  the  channel,  Ec*»  Ee*  and  the 
increase  in  ic  due  to  enhanced  current  in  the  channel  base  will  be  much  greater  than  the 
corresponding  increase  in  ie-  The  displacement  current,  essentially  a  combination  of  ic 
and  ie,  will  hence  be  more  positive  than  at  higher  pressure.  This  generally  explains  the 
novel  shape  of  the  i^  records  at  low  attachment  rate  (pressures). 

This  can  be  characterized  by  plotting  the  ratio  of  the  positive  peak  in  the  dis¬ 
placement  current  to  the  negative  peak  in  the  displacement  current  against  pressure 
(Fig.4.9).  The  correlation  is  complicated  by  the  difference  in  streamer  arrival  times 
since  the  velocity  is  also  a  function  of  the  pressure  but  the  general  trend  is  clear. 
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4.1.5b  Time-to-Peak 

The  balance  of  attachment  and  growth  influences  time-to-peak  of  the  conduction 
current  as  well.  In  fact,  Fig.4.10,  at  extremely  low  pressure  and  high  E/p,  the  peak  in 
the  conduction  current  actually  occurs  after  that  in  the  displacement  current!  This  can 
be  explained  at  follows:-  Increased  E/p  would  be  expected  to  cause  an  enhanced  stream 
growth  rate.  At  higher  pressures  this  would  be  swamped  by  the  high  attachment,  but  re¬ 
ducing  the  pressure  (attachment)  allows  more  of  this  increase  to  survive  as  far  as  the 
needle.  The  point  where  attachment  dominates  over  growth  occurs  later  and  later  in  the 
discharge,  ultimately  (as  observed)  producing  a  time-to-peak  longer  than  the  streamer 
transit  time.  This  process  is  represented  by  Eq.3.15;  the  relationship  between  the  time 
to  peak  and  the  growth  rate  is  discussed  later. 

4.1,5c  X-like  Event 

There  is  evidence  (Fig.4.10)  for  a  process  somewhat  similar  to  the  5-process  des¬ 
cribed  by  Marode  (1975a):  A  short  almost  vertical  rise  in  the  conduction  current  is 
observed  to  occur  simultaneously  with  the  burst  of  current  already  associated  with  the 
arrival  of  the  streamers  at  the  cathode  (section  3.4).  This  is  occasionally  present  in 
measurements  of  this  kind,  where  the  peak  in  ic  occurs  after  the  streamer  transit.  But 
the  similarity  is  not  perfect.  In  Marode's  experiment,  the  5  -  event  produces  the  highest 
conduction  current  in  the  discharge,  whereas  in  the  present  extreme  case  the  conduction 
current  continues  rising.  This  bears  on  the  debate  surrounding  the  establishment  of  the 
glow  in  the  channel. 

Marode  argues  that  the  8 -event  accompanies  a  potential  redistribution  along  the 
channel,  the  compensation  phase.  This  is  supposed  to  be  due  to  the  difference  between 
the  cathode  fall  voltage  and  the  streamer  tip  potential  travelling  down  the  channel  as  an 
ionising,  or  return  wave,  after  which  a  resistive  phase  is  entered,  where  there  is  current 
continuity  along  the  entire  channel.  Marode's  measurements»indeed  support  this  view. 
There  cannot  be  current  continuity  in  the  present  case,  as  (Fig.4.10)  the  current  rises 


following  the  &  - type  event  observed,  as  discussed  below. 


4.1.5d  Current  Continuity  in  the  Compensation  Phase? 

Applying  the  analysis  of  section  3,  the  currents  induced  to  the  various  electrodes 
can  be  found  if  there  is  current  continuity  along  an  axial  channel,  for  letting  the  current 
everywhere  be  j,  one  can  write  for  the  displacement  current: 

And  thus  it  seems,  from  a  fundamental  analysis  of  the  induced  currents,  that  for  a 
single  axial  filament  if  the  current  is  the  same  throughout  the  gap,  then  there  would  be 
identically  zero  displacement  current.  The  above  analysis  of  course  needs  extending  to  a 
branched  streamer  system.  Considering  a  segment  r  :  r ^<r<r2  containing  m  channels  of 
current  where  there  is  continuity  of  current  =  J),  if  one  notes  that  the  first  term 
in  the  integral  in  Eq.  4.1  does  not  depend  on  the  radial  position  in  the  gap  (r  is  an  axial 
coordinate!),  and  that  the  f(r)  term  will  only  have  an  appreciable  value  near  the  point 
where  there  will  be  a  single  discharge  channel  (and  hence  the  radial  variation  in  f(r)  can 
also  be  ignored),  the  displacement  curent  Siy*  for  the  segment  involved  can  be  written: 

k*i  it 

With  r^O,  this  evaluates  to^ ^Tand  for  ri=0.  So  it  can  be 

seen  that  for  the  whole  discharge,  the  incremental^*  add  up  to  (d-n)J/d,  which  nicely 
cancels  the  -(d-n)J/d  term  from  the  needle-end  (r^=0).  Thus  it  follows  that  if  there  is 
current  continuity,  there  will  be  no  displacement  current  whatsoever. 

In  all  the  events  recorded,  there  is  always  appreciable  displacement  current  well 
after  the  streamers  reach  the  cathode,  and  any  return  wave  h8s  long  since  decayed 


away.  Therefore  it  can  be  inferred  that  current  continuity  is  NOT  established  along  the 
channel,  conflicting  with  the  observations  by  Marode  in  small  discharges. 

So  this  S-like  event  does  not  establish  current  continuity  along  the  channel  as  in 
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Marode's  work,  but  it  does  coincide  with  the  cathode  impact.  A  return  wave  has  already 
been  associated  with  the  perturbation  in  the  conduction  current  (section  4.1.3b).  It 
seems  reasonable  to  argue  that  this  S-like  event  is  also  a  return  wave,  but  that  the  re¬ 
latively  steep  rise  at  these  low  pressures  and  high  reduced  fields  is  probably  due  to  the 
greater  conductivity  of  the  channel  in  this  case. 

4.1.5e  Attachment  Rate 

There  is  a  tendency  at  low  pressure  and  high  E/p  to  exhibit  a  slight  drop  in  attach¬ 
ment  frequency  from  the  initial  value,  evident  in  Fig.4.9,  which  is  exaggerated  at  higher 
reduced  field  and  lower  pressures  (Fig  4.11).  It  is  suggested  that  the  appropriate  value  of 
the  attachment  frequency  is  the  later  one,  (a)  because  it  is  constant  for  much  more  of 
the  record  than  in  the  very  first  portion  (except  in  very  extreme  conditions)  and  (b) 
because  the  charge  generation  is  more  likely  to  be  zero  in  the  later  portion  of  the  curve. 

A  possible  explanation  for  the  decay  in  attachment  frequency  in  the  early  stages 
of  the  discharge  is  that  it  has  been  found  from  the  reduction  of  Qg.  In  the  early  stages, 
particularly  at  low  pressure  where  the  conduction  current  is  appreciable  for  longer  on 
account  of  the  reduced  attachment  rate,  there  will  be  an  additional  reduction  in  Qg  Are 
to  the  physical  removal  of  charge  from  the  channel.  So  Qg  will  drop  a  little  more  rapidly 
at  the  start  of  the  decay  of  the  channel,  over  and  above  that  expected  from  the  attach¬ 
ment  observed  later  on.  An  extreme  example  of  this  process  is  shown  in  Fig.4.12. 

Throughout  this  study,  the  attachment  rate  has  been  found  by  using  a  simple  least - 
squares-fit  to  the  Qg  data.  In  all  but  the  most  extreme  conditions  this  is  entirely  reason¬ 
able,  for  the  maximun  charge  extracted  from  the  base  of  the  channel  beween  samples 
might  be  ~320pC  with  a  16ns  sample  time  and  2&n  A  peak  conduction  current.  The  cor¬ 
responding  value  of  Qg  can  be  as  high  as  800nC,  though  something  more  like  50nC  is 
more  typical,  and  might  decay  at  typically  3MHz  in  16n?->2.3nC.  So  except  at  parti¬ 
cularly  low  attachment  rates,  where  the  conduction  current  persists  for  a  very  long  time 
and  Qg  decays  by  less  in  a  sample  interval  due  to  this  attachment,  no  difficulty  is  ex- 
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pected  to  arise  from  the  loss  of  charge  from  the  channel  base  and  the  least -squares-fit 
method  is  satisfactory. 

Many  different  streamer  current  profiles  are  seen,  depending  on  the  field, 
pressure  and  propagation  distance.  But  it  can  be  seen  from  the  foregoing  that  the 
general  features  of  all  of  them  are  similar.  The  variation  in  the  shape  of  the  individual 
curves  is  largely  due  to  the  different  attachment  rates  and  transit  times. 

4.1.6  Growth  Frequency 

The  first  thing  one  notices  about  the  Q  records  is  the  amount  of  noise.  The 
reason  for  this  is  that  in  calculating  Q  comparison  is  made  between  adjacent  samples  of 
Qq,  itself  produced  from  adding  "simultaneous"  samples  of  ic  and  i^.  The  scatter  in  &Q 
therefore  represents  sampling  error. 

Fig.4.5  shows  a  sequence  of  Qg  and  So  records  as  the  streamers  first  of  all  just 
cross  the  gap  and  then  grow  successively  more  and  more  strongly.  The  growth  of  a  large 
and  extremely  narrow  spike  in  $Q  is  very  satisfying  evidence  that  the  charge  injection 
from  the  cathode  has  been  very  successfully  separated  from  the  background  of  the  free 
charge  remaining  on  the  channel. 

The  general  behaviour  of  Qg  from  the  simple  model  outlined  in  section  3.5.6  is 
clearly  shown  in  this  sequence,  and  shows  up  in  $Q  also.  The  sub-stability  field 
streamers  die  all  the  way  to  the  cathode  and  this  shows  up  as  a  reducing  &  Q. 

Likewise  iQ  rises  for  growing  systems.  These  curves  are  found  by  using  the  attachment 
frequency  calculated  in  the  decaying  channel.  In  this  case  the  loss  of  charge  from  the 
base  of  the  channel  has  been  explicitly  taken  into  account  as  is  plain  if  the  defining 
equation  for  is  recalled: 

Q,  •  («,-  i.<tt  e  +  Sq. 

'  Jt(  (3.14) 

Thus  the  difficulty  raised  above  with  the  attachment  rate  at  low  pressure  has  been 
overcome.  In  this  case  the  correction  matters,  for  in  the  early  portion  of  the  discharge 
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there  is  more  conduction  current  than  later  on  when  the  attachment  is  found. 

The  synchronizational  error  responsible  for  a  single  narrow  negative  spike  at  the 
outset  of  the  Qg  record  seems  to  the  $Q  algorithm  to  be  a  sudden  injection  of  charge, 
followed  by  a  sudden  removal  of  charge.  The  result  is  that  there  is  a  large  negative  and 
then  positive  spike  in  S(3.  A  positive  spike  in  Qg  will  produce  a  positive  and  then  a 
negative  spike  in  £q.  These  correspond  to  the  single  spike  in  Qg  which  was  an  artefact, 
and  (again)  are  ignored. 

Finally,  one  notes  the  slight  tendency  for  5(3  to  edge  above  the  axis  right  after 

the  streamer  transit  (Fig.4.13).  This  is  due  to  the  attachment  frequency  falling  slightly 

from  its  value  immediately  after  impact  with  the  cathode.  As  discussed  above,  this  drop 

is  due  to  residual  conduction  current  making  the  attachment  seem  greater  early  on.  The 

value  of  the  attachment  frequency  used  in  calculating  <5q  is  lower  than  the  actual  value 

just  after  cathode  impact.  There  will  be  more  charge  removed  from  the  channel  than 

predicted  using  the  lower  attachment  coefficient.  Since  Sq  is  thought  to  be  zero  in  this 

region,  using  Eq.  3.14  to  find  $Q  will  result  in  a  positive  value.  v 

\ 

The  amount  of  information  represented  by  the  5(3  curves,  information  further¬ 
more  directly  related  to  the  growth  of  the  streamer,  and  corrected  for  the  removal  of 
charge  from  the  base  of  the  channel;  such  as  the  initial  decay,  take-off  of  growth  or 
decay  depending  on  the  external  field,  and  charge  injection  from  the  cathode,  shows  that 
of  the  available  quantities,  this  best  represents  the  growth  of  the  streamer. 


4.2  Pre-breakdown  Corona  Modes 

As  the  breakdown  field  is  approached,  new  corona  modes  are  observed  normally 
involving  a  change  in  the  observed  current.  These  phenomena  have  the  appearance  of 
choked-off  breakdown  processes. 

There  are  a  number  of  different  prebreakdown  modes;  the  corona  can  exist  in 
distinct  states,  which  can  be  stable  over  periods  of  minutes,  or  can  switch  intermittently 
from  one  state  to  the  other  every  few  seconds  or  so.  In  an  attempt  to  elucidate  the  dif- 
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ferent  conditions  under  which  the  various  modes  were  exhibited,  a  case  study  was  made 
of  the  500torr  behaviour  and  this  is  presented  in  section  5. 

The  general  form  of  the  various  modes  is  discussed  here.  It  should  be  remembered 
that  the  discharge  current  growth  is  limited  by  the  series  resistor. 

4.2.1  Terminology 

Since  some  of  the  mode  parameters  are  peculiar  to  the  present  electrode  geo¬ 
metry  some  definitions  may  prove  useful.  The  classic  streamer  discussed  above  is  called 
primary  streamer  to  make  plain  that  this  is  the  usual  mode  observed  (at  low  reduced 
field).  Secondary  streamers  as  identified  by  Marode  (1975),  Hudson  &  Loeb  (1961),  etc. 
are  not  studied  in  this  section  -  they  properly  belong  to  a  discussion  of  classic  streamer 
development.  All  the  prebreakdown  modes  involve  an  increase  in  the  gap  current. 

Mode  ii  This  has  been  termed  "second  corona"  and  has  a  large  conduction  current  event 
growing  out  of  the  primary  streamer  (before  the  channel  has  quite  decayed). 

Mode  iii  Called  a  "delayed  second  corona"  as  it  has  some  features  in  common  with  Mode 
(i).  As  might  be  expected,  there  is  a  distinct  period  between  the  primary  and  the  delayed 
second  corona  where  no  current  growth  occurs. 

Mode  iii;  Called  'Vnultiple  second  corona"  or  the  'Ynultiple"  or  "high-frequency"  mode 
since  this  represents  the  main  characteristic  of  this  mode;  it  is  basically  several  events 
of  type  (ii)  in  rapid  succession.  Fig.  4.14  summarizes  the  differences  between  the  various 
types. 


4.2.2  Mode  (i):  Second  Corona 

As  the  field  increases,  the  first  new  phenomenon  observed  following  the  return 
wave  involves  a  discontinuous  jimp  in  the  conduction  current  (Fig.  4.15).  This  always 
occurs  something  like  50  to  100ns  after  the  streamer  transit,  and  cannot  therefore  be 
simply  a  large  return  wave  (see  for  example  Fig. 4.15)  where  a  perturbation  in  the  con- 
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duction  current,  coinciding  with  the  return  wave  proper,  can  clearly  be  seen  80ns  before 
the  regeneration  of  the  channel). 


4. 2. 2. a  Regeneration  of  the  Channel 

The  second  corona  current  is  large,  compared  with  that  of  the  primary  streamer. 
The  first  feature  of  interest  in  the  second  corona  is  that  the  current  undergoes  a  large 
and  fast  regeneration. 

Meanwhile  the  electron  current,  i  which  was  dropping  away  very  rapidly  before 
this  regeneration  of  the  channel  shows  a  tiny  increase  by  comparison.  Since  both  i  and 
i  .  are  changing  very  rapidly,  and  in  opposite  directions,  even  this  feature  is  probably  a 
synchronization  artefact. 

The  displacement  current  shows  a  positive  excursion  simultaneous  with  the  re¬ 
generation  of  the  conduction  current;  as  discussed  in  section  3,  this  shows  unambiguously 
that  this  event  is  confined  to  the  region  of  the  needle-end.  Furthermore  if  such  a  large 
i  record  can  be  obtained,  with  little  or  no  corresponding  signal  in  ie,  then  the  spatial 
extent  of  f(r)  must  be  quite  small  indeed,  confirming  the  position  taken  in  section  3. 

4.2.2.b  Appearance  of  the  second  corona 

Visual  observations  of  the  discharge,  when  the  current  records  exhibit  this 
feature,  show  it  to  contain  a  single  bright  channel,  in  the  lower  portion  of  the  gap,  ex¬ 
tending  from  the  point.  This  regeneration  of  the  channel  base  clearly  represents  the 
early  phase  of  a  regeneration  of  the  whole  of  a  single  filament  in  the  lower  channel  (see 
below  under  "growth  of  second  corona",  section  4.2. 2. e) 

Qualitative  observations  suggest  that,  like  return  waves,  this  phenomenon  also  is 
more  readily  expressed  at  lower  pressures.  At  higher  pressures  the  channel  can  never¬ 
theless  be  regenerated,  but  with  lower  probability. 
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A.2.2.C  Attachment  rate/heatinq  from  the  discharge? 

It  can  readily  be  seen  from  the  curves  of  such  an  event  that  the  attachment 
frequency  is  lower  than  before;  the  slope  of  the  curve  is  considerably  shallower  after  the 
second  event.  Owing  to  the  possibility  that  the  channel  has  been  heated,  and  therefore 
that  the  electron  drift  velocity  (because  of  the  change  in  pressure)  is  undefined,  the  new 
attachment  coefficient  cannot  be  determined  from  the  attachment  frequency.  Using  the 
old  attachment  frequency  to  calculate  Sq  leads  to  a  negative  displacement  of  the  second 
portion  of  the  curve.  The  old  attachment  frequency  implies  a  net  negative  in-jection  of 
charge  -  a  consequence  of  the  lower  real  attachment  frequency.  This  feature  is  shown 
unambiguously  by  the  low  pressure  event  in  Fig.  4.16. 

4.2.2.d  Detachment  or  ionization? 

A  question  to  be  addressed  is  whether  or  not  this  very  rapid  growth  in  electron 
density  is  associated  with  a  new  ionizing  regime,  or  with  some  process  whereby  electrons 
detach  from  pre-existing  negative  ions. 

Ionization  and  detachment  both  lead  to  an  increase  in  the  electron  charge.  The 
distinction  between  the  two  is  rather  more  subtle,  the  electrons  have  been  removed  from 
different  species  and  thus  the  energetics  and  dependence  on  gas  composition  or  ionic 
state  became  important:  detachment  from  negative  ions  may  proceed  by  a  larger 
nanber  of  routes  than  the  ionization  of  neutrals,  and  requires  less  energy  (Elliason  & 
Kogelshatz,  1986). 

Neverthless,  the  peak  in  the  conduction  current  associated  with  the  regeneration 
of  the  channel  base  is  larger  than  the  initial  conduction  current  peak.  This  difference  in 
size  can  at  times  be  exceptionally  large  (see  e.g.  Fig.4.18)  and  it  seems  unlikely  that 
detachment  from  the  negative  ions  created  during  the  transit  of  the  primary  streamer 
can  account  for  all  the  electrons  freed  in  the  second  corona.  Hence  some  fresh 
ionization  must  occur  to  provide  the  balance  of  the  electrons. 

Finally  the  comparison  with  secondary  streamers  must  be  ruled  out  categorically; 
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secondary  streamers  are  entirety  a  decaying  current  regime.  Spectroscopic  evidence 
(Hartmann,  1974)  shows  the  electron  energy  is  too  low  for  there  to  be  any  ionization  in 
the  secondary  streamer.  This  feature  is  associated  with  a  decaying  channel  from  the 
primary  streamer  and  would  always  follow  primary  streamers  which  crossed  the  gap.  The 
second  corona  is  not  always  produced  and  is  an  increasing  current  regime,  and  cannot, 
again,  be  a  secondary  streamer. 

4.2.2.e  Growth  of  Regenerated  Channel 

After  the  initial  rapid  regeneration,  confined  to  the  base  of  the  channel  ie  then 
begins  to  increase,  though  more  slowly,  producing  a  broad  peak  in  the  i  record 
(Fig.4.15).  There  is  often  a  time  delay  of  some  ns  be,»re  this  occurs  and  further  attach¬ 
ment  in  the  body  of  the  gap  reduces  i  in  the  interim.  The  attachment  frequency  during 
this  period  does  decrease  slightly  as  can  be  seen  in  Fig.4.15  by  comparing  the  decay 
curve  before  and  after  the  anall  perturbation  (which  in  this  case  is  probably  an  artefact 
as  discussed  in  section  4.2.2.a) 

This  broad  peak  in  ie>  unrepresented  in  ic,  shows  (section  3.4)  that  the  ionization/ 
detachment  processes  are  occurring  beyond  the  point  region  and  by  implication  suggests 
that  a  new  ionizing/ detaching  regime  is  developing  probably  along  the  single  bright 
channel  observed  visually.  Other  evidence  for  a  growing  channel  is  the  width  of  the 
second  corona  conduction  current  pulse.  If  the  current  were  all  confined  to  the  base  of 
the  channel  -  as  it  clearly  is  at  the  outset  of  this  event  since  thre  is  no  correspondingly 
large  electron  current  signal  -  then  ic  would  contain  only  an  impulsive  signal  -  the  elect¬ 
rons  would  drift  to  the  anode  very  quickly.  In  order  to  provide  the  electrons  to  sustain 
the  conduction  current  (observed  in  second  corona  records  to  last  quite  as  long  as 
primary  streamer  conduction  currents),  there  must  be  a  channel  growing  from  the  point. 
All  this  points  to  the  spatial  growth  of  the  second  corona  event. 

The  delay  between  the  inception  of  the  second  corona  and  the  rise  of  the  broad 


peak  in  the  electron  current,  already  touched  upon  above,  is  a  very  interesting 
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phenomenon.  After  the  perturbation  associated  with  the  synchronizational  error  when 
the  second  corona  occurs,  the  attachment  frequency  can  clearly  be  seen  to  reduce 
slightly.  Once  the  second  corona  has  begun  to  propagate  into  the  gap,  its  growth  picks  up 
very  strongly,  indicating  that  some  property  of  the  gap  is  enhancing  its  growth.  The  view 
advanced  here  is  that  the  remaining  electrons  in  the  "old1  channel  are  responsible.  The 
density  of  free  electrons  near  the  point  will  be  low,  as  they  are  most  effectively  swept 
out  of  the  gap  here.  The  electrons  at  the  point  will  have  been  generated  further  out  in 
the  gap  as  the  time  lag  between  first  and  second  corona  increases,  and  consequently  will 
be  more  heavily  attached.  The  density  immediately  near  the  point  will  thus  increase  as 
the  second  corona  propagates  out  into  the  gap.  This  will  explain  the  enhanced  growth  as 
the  second  corona  travels  away  from  the  point;  the  conductivity  increases  as  the  point  is 
left  behind.  The  larger  free  electron  density  will  mean  that  a  streamer  tip,  for  example, 
will  not  need  to  produce  secondary  electrons,  and  will  propagate  very  much  more  vigor¬ 
ously  than  in  a  lower  electron  density. 

No  information  on  the  speed  of  this  phenomenon  is  available,  as  the  distance  it  has 
propagated  is  not  known,  and  streak  or  photomultiplier  techniques  were  not  used.  But 
visual  observations  suggest  that  the  length  of  the  new  channel  could  be  about  half  the 
gap  length,  occasionally  more,  and  the  current  records  show  this  occurs  in  times  of  about 
$  to  J  of  a  streamer  transit  time  suggesting  a  speed  about  the  same  or  slightly  less  than 
a  classic  primary  streamer. 


4.2. 2. f  Explanation? 

This  raises  the  intriguing  possibility  that  this  is  in  fact  a  new  streamer,  forming 
after  a  delay  much  shorter  than  the  ~ lOOps  expected  from  observations  of  primary 
streamers,  and  guided  into  a  single  channel  by  some  mechanism.  Certainly  the  regener¬ 
ation  responsible  for  it  produces  a  conduction  current  pulse  larger  than  seen  from 
primary  streamers  and  once  these  electrons  have  been  neutralized  at  the  needle-end, 
enough  positive  charge  will  be  left  to  produce  a  streamer  tip.  Though  how  concentrated 
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this  is,  is  again  left  open. 

The  growth  of  the  regenerated  channel,  and  its  propagation  into  the  gap  is  so 
similar  to  the  features  of  primary  streamers  that  this  new  corona  mode  has  been  termed 
a  "second  corona",  to  distinguish  it  from  a  secondary  streamer  and  to  tentatively  asso¬ 
ciate  it  with  "subsequent  corona”  in  long  point-plane  gap  impulse  work.  Certainly  the  se¬ 
quence  of  events;  (i)  a  burst  of  conduction  current  identified  as  being  near  the  tip,  (ii) 
free  electrons  generated  in  the  body  of  the  gap  and  (Hi)  a  continuing  supply  of  electrons 
to  the  point  from  the  discharge;  is  strongly  suggestive  of  at  least  an  ionizing  or  detaching 
event  which  propagates  out  into  the  gap. 

As  the  field  and  gap  length  increases  (section  5.3.3),  successively  larger 
regenerations,  or  further  coronae  are  seen  (Fig.  4.17).  This  violent  growth  of 
conductivity  in  the  gap  is  strongly  indicative  of  a  stifled  breakdown.  This  topic  is 
considered  in  more  detail  in  section  6. 

4.2. 2.q  Second  Corona  Reach  Cathode? 

An  Intersting  feature  remains  to  be  discussed.  Occasionally  a  small  increase  in  ie 
(Qg)  is  seen  (Fig.4.18)  in  the  later  stages  of  the  propagation  of  the  second  corona.  This  is 
not  present  in  ic  and  so  by  analogy  with  primary  streamers,  can  the  second  coronae 
possibly  have  reached  the  cathode?  This  feature  would  then  represent  electron  injection 
from  the  cathode. 

This  is  discounted  for  two  reasons,  firstly  the  feature  is  too  small:  charge 
injection  from  the  cathode  due  to  the  primary  streamer  arrival  is  much  greater,  and 
secondly  there  is  no  sign  of  a  return  wave  signature  in  the  conduction  current  record,  as 
has  already  been  shown  to  be  present  when  the  primary  streamer  reaches  the  cathode. 

As  discussed  in  section  3  this  last  feature,  having  no  corresponding  increase  in  the 
conduction  current  record  cannot  be  occurring  anywhere  near  the  needle-tip,  nor  as  just 
indicated  is  it  identified  with  the  arrival  of  the  second  corona  streamers  at  the  cathode 
and  so  it  is  thought  that  it  must  be  a  further  ionizing  event  started  entirely  in  mid-gap. 
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4.2.3  Mode  (ii);  Delayed  Second  Corona 

With  a  low  trigger  level  only  the  first  streamer  mode  was  seen,  but  by  turning  up 
the  oscilloscope  trigger,  on  a  short  time-base,  higher  current  pulses  could  be  observed, 
suggesting  a  bi-modal  distribution  of  corona.  These  pulses  had  conduction  currents  which 
rose  on  timescales  comparable  with  primary  streamers,  but  were  of  larger  amplitude. 

The  displacement  current  as  a  result  was  more  positive-going  and  less  negative-going 
than  its  primary  streamer  counterpart. 

Observations  with  longer  timebase  showed  that  rather  than  there  being  a  bimodal 
streamer  distribution,  the  second  corona  (Fig-4.19)  was  occurring  up  to  5ps  after  the 
primary  streamer  had  completely  decayed.  In  cases  where  the  second  corona  seemed  to 
be  in  isolation,  it  is  inferred  that  the  first  corona  had  escaped  the  oscilloscope 
"window".  The  delay  variability  is  discussed  below  (section  5.3).  The  visual  appearance 
of  the  corona  in  this  mode  was  exactly  the  same  as  for  classic  streamers,  and  contained 
none  of  the  bright  channels  observed  in  the  more  rapidly  generated  second  coronae  dis¬ 
cussed  above.  The  likelihood  of  observing  this  delayed  second  corona  was  very  low  com¬ 
pared  with  that  of  observing  the  accompanying  classical  corona  in  isolation. 

4.2.4  Mode  (iii):  Multiple  second  corona 

The  corona  could  switch  into  a  third,  high  frequency  mode,  where  not  just  a  few, 
but  very  corona  pulse  displayed  multiple  second  corona-like  pulses.  This  mode  was 
intermittent  and  could  cut  in  and  out  irregularly.  A  typical  current  record  is  shown  in 
fig.  4.20.  The  intermittency  suggests  it  is  dependent  on  some  gross  parameter  of  the 
gap,  which  can  stabilize  the  mode  -  the  discharge  had  to  be  running  for  a  while  before 
this  mode  appeared.  This  irregularity  made  observation  difficult  though. 

The  visual  appearance  of  the  corona  in  this,  third  mode  was  different  from  that 
observed  when  an  isolated  second  corona  pulse  was  present.  It  appeared  exactly  the 
same  as  normal  corona,  except  now  the  whole  discharge  wsmuch  brighter  than  usual. 
Isolated  trains  of  subsequent  corona  were  observed  right  down  to  onset,  though  with  low 
probability  (see  case  study,  section  5.3).  The  high  frequency  mode  is  clearly  just  this 


phenomena  occurring  repetitively. 

The  difference  in  the  appearance  of  the  corona  could  possibly  be  attributed  to  the 
second  corona  in  this  case,  being  started  in  an  inactive  channel.  When  the  event  causing 
the  inception  of  the  second  corona  occurs  in  a  channel  with  significant  electron  density, 
perhaps  there  is  a  tendency  to  guide  the  corona  along  the  highest  current  channel,  where¬ 
as  if  the  channel  has  decayed  so  far  that  there  is  insignificant  electron  density,  then 
perhaps  the  guiding  effect  is  less  pronounced. 

4.3  Breakdown 

Experiments  to  determine  the  breakdown  field  were  carried  out,  with  the  instru¬ 
ments  disconnected,  and  the  charging  resistor  out  of  the  circuit.  Determining  the  break¬ 
down  field  is  not  an  easy  measurement  with  a  d.c.  corona,  for  the  high  repetition  rate 
means  that  very  low  breakdown  probabilities  can  be  detected.  Defining  the  lower  limit  is 
virtually  impossible.  One  breakdown  in  a  five  minute  period  was  observed  at  3.62kV/cm, 
(6V/cm.torr  at  600  torr)  well  below  the  usual  level  of  about  7kV/cm  (or~9.2V/cm.torr  at 
S.T.P). 

The  probability  of  second  corona  was  also  investigated  and  the  results  (Fig.4.21) 
plotted  along  with  the  breakdown  field.  The  suggestion  is  that  the  formation  of  second 
corona  would  ultimately  lead  to  breakdown,  but  is  prevented  from  so  doing  by  the  series 
resistor.  This  is  certainly  the  picture  at  high  pressure. 

At  lower  pressures,  a  further,  interesting  observation  is  that  the  gap  will  not 
break  down  even  when  second  corona  is  observed,  and  the  resistor  is  out  of  the  circuit. 

The  explanation  may  lie  in  the  availability  of  energy  in  the  gap.  To  obtain  a  part¬ 
icular  reduced  field  does  not  require  as  high  a  voltage  on  the  cathode  at  lower  pressures 
as  it  does  at  higher  pressures.  So  the  stored  energy  in  the  gap  is  less  at  lower  pressure. 
Since  a  breakdown  event  is  faster  than  the  10ns  or  so  required  for  the  recharging  of  the 
gap  capacitance  then  the  principal  source  of  energy  for  the  breakdown  is  the  stored 
energy  in  this  capacitance.  At  very  low  pressures  not  enough  energy  is  present  at  re- 
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duced  field  otherwise  capable  of  producing  breakdown  in  the  absence  of  an  energy  limit. 

4.3.1  Energy  Limit 

In  the  event  shown  in  Fig.  4.16  (#418.3),  the  total  charge  drawn  by  the  external 
circuit,  Jledt +  n/d/»cdt,  is~23.6nC,  compared  with~110nC  stored  on  the  gap  capaci¬ 
tance  at  this  field.  This  amounts  to  something  like  a  20%  drop  in  applied  field  and  is  sig¬ 
nificant  enough  to  explain  the  failure  of  the  gap  to  breakdown.  Note  though  that  this  is 
not  Qg,  which  may  be  a  good  deal  higher  than  this  figure  before  the  applied  field  begins 
to  drop,  depending  on  J"icdt. 

4.4  Summary 

The  first  section  was  devoted  to  classic  corona,  corresponding  in  general  with  the 
description  to  be  found  in  the  literature  of  "positive  corona  streamers".  The  general 
features  of  current  observations  of  these  are  discussed,  in  the  light  of  the  induced 
current  analysis  developed  in  section  3.  The  wide  variability  of  the  appearance  of  the 
streamer  currents  as  the  field,  pressure  and  propagation  length  are  varied,  is  shown  to  be 
mainly  tkie  to  the  same  phenomena,  differences  being  due  to  changes  in  the  attachment 
frequency  and  transit  time. 

Thereafter,  it  was  found  that  as  the  reduced  field  increases,  different  corona 
modes  are  seen.  Indicative  of  further  current  growth  occurring  after  the  primary 
streamer  channel  has  begun  to  decay  these  phenomena  are  associated  with  pre¬ 
breakdown.  Close  analysis  of  the  induced  currents  (again  following  section  3),  shows  the 
mode  occurring  at  the  highest  field  to  be  due  to  a  streamer-like  event,  born  at  the  needle 
and  subsequently  propagating  out  into  the  gap.  Visually  this  mode  (mode  (i))  is  associated 
with  the  presence  in  the  discharge  of  discrete  bright  channels,  which  are  not  present  at 
lower  reduced  fields.  This  mode  appears  to  depend  upon  the  electron  density  in  the 
channel.  The  major  difference  beween  it  and  the  lower  reduced  field  modes  (ii  and  iii) 


with  longer  inception  times,  is  that  in  the  latter  case  the  channel  has  almost  decayed 
before  inception. 
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Breakdown  measurements  with  the  high  voltage  charging  resistor  either  in  or  out 
of  the  circuit  (Bicknell  &  Shelton,  1985)  show  that  the  breakdown  is  being  choked-off  by 
the  shortage  of  energy  in  the  observations  just  described;  the  gap  would  usually 
breakdown  if  the  resistor  were  not  present. 

The  next  section  considers  in  detail  the  way  in  which  parameters  characterizing 
the  classic  discharge  are  affected  by  the  ambient  conditions  and  particularly  by  the  prop¬ 
agation  length.  A  case  study  has  been  made  of  the  many  different  corona  modes  in  the 
pre-breakdown  regime,  to  try  and  elucidate  their  dependence  on  these  same  parameters. 
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Section  3:  Experimental  Results 


Introduction 

In  this  section  the  hulk  of  the  data  obtained  is  presented  in  graphical  form.  The 
considered  parameters  include  peak  conduction  current,  streamer  velocity,  attachment 
frequency,  maximum  free  charge,  electron  injection  from  the  cathode  and  attachment 
frequency  (coefficient).  An  interactive  computer  program  was  used  to  analyse  the  data 
(Appenriix/S)  which  v/as  stored  in  a  data  base  reproduced  in  Appendix  B.  The  data  refers 
to  v.hat  are  termed  here  "classic  streamers”  or  those  streamers  generated  between  onset 
and  breakdown.  As  breakdown  is  appruached  the  streamer  behaviour  becomes  more  com¬ 
plex  and  this  behaviour  is  considered  as  a  separate  item  at  thr  end  of  the  section.  Where 
appropriate  sor,  e  discussion  is  presented  but  more  detailed  discussion  is  reserved  for  sub¬ 
sequent  sections  as  indicated. 

5.1  Treatment  of  Data 

Inmost  cases,  four  corona  events  v  ere  analysed  for  each  set  of  experimental  con¬ 
ditions;  gap  length,  pressure  and  field.  F or  consistency,  dry  air  v/as  used  throughout.  For 
some  combinations  of  parameters,  10  or  more  events  v.ere  analysed.  In  some  cases,  not 
all  events  recorded  could  be  analysed  fully  and  only  some  of  the  derived  quantities  are 
available.  For  example,  while  all  streamers  have  a  velocity,  this  cannot  be  measured  by 
time-of-flight  for  those  which  do  not  cross  the  gap.  As  a  result  there  is  no  velocity  data 
for  certain  ranges  of  parameters.  Decause  the  reduced  onset  field  increases  with  falling 
pressure,  for  certain  fields  at  low  pressure  there  v.as  no  corona  whatever,  and  again  no 
parameters  are  available. 

Where  a  discrete  value  is  obtained,  (i.e.  all  parameters  except  attachment  and 
growth  frequencies,  then  the  mean  and  standard  deviation  of  the  measurements  were 
plotted. 

In  the  case  of  the  attacfmient  and  growth  frequencies,  being  obtained  from  a 
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least-sguares  technique,  m .  error  S’  coi;lc!  lie  obtained  for  each  experimental  point, 
Plotting}  these  individual  measurements,  since  the  range  of  experin ,  nntai  errors  was  wide 
and  the  values  not  a  little  scattered,  gives  a  confusing  picture.  Hence  weighted  mean 
of  the  means  and  standard  deviations  of  the  raw  derived  results  are  plotted,  where  the 
weight  is  the  reciprocal  of  the  fractional  error  of  the  point,  v,  j  =  x./ot.  The  means  and 
standard  deviations  then  become: 


X  ±  C  s 


(5.1) 


The  effect  of  this  weighting  is  to  emphasize  the  measurements  having  the 
smallest  error,  without  actually  eliminating  any  one  point.  A,  J  since  the  rn.eon  of  the 
individual  errors  is  used,  9  represents  the  range  of  the  observations,  not  the  error  of  the 
mean,  as  for  the  other  (discrete)  data.  Where  the  points  are  so  close  together  that  the 
error  bars  niake  the  graphs  confusing,  only  the  points  are  plotted.  When  this  is 
necessary,  the  spread  of  the  data  is  generally  wider  than  the  accuracy  of  any  one  par¬ 
ticular  measurement  anil  hence  the  underlying  variability  is  not  affected;  it  is,  on  the 
contrary,  highlighted. 

5.2  Classic  Streamers 

5.2.1  Conduction  Current 

fig. 5.1  reports  the  variation  of  peak  conduction  currents  with  El/p,  pressure  and 
propagation  length.  Data  belov.  400  torr  ore  not  plotted,  as  the  S-event  disrupts  the 
measurement. 


The  data  tenil  to  a  limiting  value  after  the  stability  field  is  attained,  except  that 
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Fig.  5.1 
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obtained  at  5.3<m,  which  may  not  have  levelled  off  at  the  reduced  fields  studied.  Within 
it  must  be  stressed  an  extremely  large  experimental  spread  and  at  il)0  torr  and  above, 
the  limiting  value  appears  not  to  depend  on  pressure,  but  nn  the  propagation  length. 

Fig. 5. 2  shows  a  plot  of  the  limiting  value  of  the  conduction  current  peak  against 
propagation  length,  d.  It  shov/s  a  minimum  at  somewhere  around  15  to  16  cms.  Though  is 
must  be  observed  that  the  range  of  currents  of  between  B  and  2LhnA  is  not  wide,  digital 
sampling  error  could  contribute  significantly  to  the  uncertainty  in  the  maximum  value  of 
i  measured,  and  the  large  error  involved  may  mean  that  the  rise  in  i^/max)  at  higher  d 
may  not  be  significant.  The  increase  in  the  limit  at  low  d  though,  is  too  steep  to  he 
ignored. 

Marode  (1975a)  reports  a  typical  streamer  current  maxima  of  around  50mA  or  so, 
for  prebreakdown  point-plane  systems  1cm  or  so  in  extent.  Though  the  influence  of 
vibrational  states  may  be  invoked  to  explain  some  of  this  increase  (from  enhanced  ion¬ 
ization,  or  other  mechanisms),  ttie  magnitude  of  the  increase  does  in  general  confirm,  the 
trend  to  higher  maxima  at  lower  gap  length. 

The  difference  compared  with  the  point-pione  work  of  b  arode  is  particularly 
interesting,  since  in  point-plane  work  it  is  considered  that  the  corona  current  is  de¬ 
pendent  on  the  potential  of  the  tip,  rather  than  the  average  field.  So  it  might  be  argued 
that  the  three  electrode  system  used  is  sufficiently  similar  to  point-plane  geometry  that 
the  actual  voltage  on  the  needle  is  crucial.  The  gap  in  the  present  study  is  "inverted" 
with  the  voltage  applied  to  the  cathode,  but  nevertheless  the  voltage  difference  to  estab¬ 
lish  a  particular  field  v/ill  increase  with  the  separation  of  the  electrodes.  Analogy  with 
point-plane  geometry  would  imply  a  uniformly  increasing  streamer  current  at  the  same 
field  as  the  gap  length  goes  up.  The  distinct  fall  in  current  as  distance  increases  very 
clearly  shov/s  this  analogy  to  be  groundless.  (It  can  also  be  argued  that  at  a  given  E/p 
and  d  the  voltage  point  above  the  cathode  voltage  must  drop  ns  the  pressure  drops,  to 
maintain  C/p.  This  would  imply  a  lower  value  of  ic(max)  using  the  analogy  with  point- 
plane  findings  -  again  this  is  not  seen  in  the  Figure.) 
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The  topic  of  the  conduction  current  maximum  is  token  up  in  thr  discussion,  and 
while  the  low  gap  length  ami  of  the  curve  (Fit;. 5.2,'  might  he  explicable  ns  n  geometric 
factor,  the  possibility  that  it  can  subsequently  increase  with  distanre  cannot  tie  explained 
by  supposing  residual  current  to  reach  the  point  from  the  rest  of  the  enlarged  streamer 
system,  for  the  peak  occurs  in  some  IOC  to  2C0ns  at  most.  Electron  drift  speeds  of  even 
50km/s  mean  that  at  most  the  first  cm  of  channel  only  can  be  responsible  for  deter¬ 
mining  the  peak  in  ic-  Thus  the  suggestion  of  a  variability  with  length  over  and  above  the 
first  few  cm  is  quite  extraordinary. 

5.2.2  Stream, er  Velocities 

Fig. 5.3  shows  plots  of  streamer  velocities  against  E/p,  for  various  propagation 
lengths  and  pressures,  found  by  time-of-flight.  With  the  exception  of  the  5.3cm  data,  all 
the  600  torr  velocities  lie  on  the  same  curve.  At  500  torr,  the  stody  is  the  same,  except 
that  the  tendency  to  lov/cr  streamer  velocity  at  5.3cms  is  also  shown,  though  less 
strongly,  by  the  10.3cm  data. 

Fig. 5. 4  shows  a  family  of  curves,  each  at  a  different  value  of  F/p,  showing  v$  as  a 
function  of  pressure,  for  the  18.3cm  data.  There  is  a  very  interesting  minimum  at  around 
400  torr,  but  only  for  an  E/p  of  8.5V/ cm. torr;  at  lower  E/p  there  is  only  a  point  of  inflex¬ 
ion  here. 

Streamer  propagation  velocity  is  an  easily  measured  but  clearly  very  complex 
parameter,  and  no  satisfactory  theory  of  streamer  velocity  exists.  The  rise  in  streamer 
velocity  with  reduced  field  is  self-evident,  for  the  electron  drift  velocity  will  be  higher. 
Explaining  the  minimum  in  velocity  when  plotted  against  pressure  at  a  given  reduced 
field  is  not  easy,  but  it  may  be  possible  to  show  qualitatively  how  the  minimum  comes 
about  by  considering  the  time  to  create  a  secondary  avalanche. 

A  particular  property  of  this  minimum  is  that  it  occurs  only  at  high  reduced 
field.  The  curves  might  be  considered  to  be  a  general  reduction  in  streamer  velocity  as 
pressure  increases,  with  an  E/p-dependent  enhancement  which  increases  strongly  with 
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pressure  above  5GC  torr. 


5.2.3  Free  Charge 

To  provide  an  easily  defined  parameter  to  characterize  the  streamer  growth,  the 
maximum  free  electron  charge,  D^fmax)  was  measured.  This  parameter  does  not  appear 
to  have  been  measured  in  this  way  by  any  other  investigator  and  hence  the  way  it  varies 
is  of  interest. 

The  maximum  in  Gp,  however,  is  something  of  a  hybrid,  representing  as  it  does  the 
free  charge  residing  on  the  channel  and  from,  cathodic  injection.  It  is  no  easy  task  to 
separate  the  cathodic  injection  by  sight  from  the  background  of  the  growing  streamer 
system,  and  so  the  easier-to-ibentify-parameter,  the  maximum,  v/as  measured.  This  is 
something  of  a  compromise,  hut  the  unambiguous  definition  of  the  maximum  v/as  pre¬ 
ferred  to  the  more  subjective  measurement  which  would  have  arisen  from  trying  to  de¬ 
fine  the  separation  point  betv/een  the  growth  in  and  the  charge  injection  from  the 
cathode.  Anyway,  the  charge  injection  can,  in  principle,  be  more  satisfactorily  measured 
from  the  size  of  the  narrow  spike  in  the  $G  records,  (see  below  section  5.2.4). 

Fig. 5. 5  shows  the  variation  of  G^(max)  as  a  function  of  E/p  at  various  propagation 
lengths  and  pressures.  In  the  figure  it  can  be  seen  that  after  the  stability  field  is  passed 
the  gap  charge  rises  exponentially  through  a  decade  for  an  0. IV/ cm. torr  increase  in 
E/p.  Except  for  the  5.3cm  data  this  seems  to  be  independent  of  propagation  length:  all 
the  data  lie  roughly  on  the  same  curve.  This  shows  the  strong  increase  in  charge  once 
the  system  begins  growing. 

Beyond  5V/ cm.  .torr  the  data  separates  and  a  strong  propagation  distance  depen¬ 
dence  sets  in,  with  weaker  growth  as  a  function  of  E/p.  RepJotting  the  data  os  a  function 
of  propagation  length  results  in  an  (almost)  exponential  increase  in  with  distance 
(Fig. 5. 6).  The  strong  growth  of  G^(max)  with  distance  would  suggest  either  that  the 
charge  on  the  channel,  Q^,  is  considerably  larger  than  that  injected  from  the  cathode,  or 
else  that  the  charge  injected  from  the  cathode  is  proportional  to  Q^.  This  might  be 


fcry: 

tsQDtorr 

550torr 

SOOtorr 

450torr 

400torr 

360torr 

300tonr 

250torr 


18Jan 


A 

▲ 

a 

B 

o 

e 

v 

▼ 


?OOtorT 


$ 


S' 


80 


accounted  for  ns  follov.s. 


0.2.3. a  Tip/Channol  Charge  Ratio 

Redefining  tt to  coordinate  frame  such  that  the  end  cf  the  streamer  is  at  y=0 
'rig. 5. 7;,  and  considering  the  electrode  to  he  on  infinite  distance  away  (i.e.  that  all 
electrons  attach  before  reaching  the  anode),  the  channel  function, "X(y)  can  he  written: 

O  O 

which  can  easily  he  evaluated  for  a  particular  form  of  A(y).  This  will  hold  fur  a  branched 
streamer  system,  as  well  as  for  a  single  channel;  for  considering  a  horizontal  slat  through 
the  discharge,  and  taking  the  idealized  view  of  the  pcr-unit-ficlc!  distribution  outlined  in 
section  3  the  induced  current  is  just  that  from  adding  the  currents  induced  from  t tie  in¬ 
dividual  filaments.  is  just  the  total  charge  in  all  those  filaments  and  the  attachment 
will  always  he  proportional  to  it.  The  growth  will  depend  on  the  size  of  the  individual 
streamer  tips,  hut  this  will  nffert  the  situation  by  changing  the  growth  function  A(y), 
which  is  what  one  hopes  to  measure  -  hence  there  is  no  need  to  build  into  A(y)  the  in¬ 
fluence  of  the  tip  potential  distribution. 

The  principle  question  is  whether  this  integral  will  give  a  result  proportional  to 
A(y=0),  that  is  to  the  size  of  the  tip.  Assuming  the  charge  injection  from  the  cathode  to 
depend  at  least  on  the  size  of  the  tip  also,  the  maximum  in  Gf,  evaluated  at  some  part¬ 
icular  reduced  field  would  then  be  proportional  to  the  total  streamer  tip  charge  even 
taking  the  channel  in  to  account.  With  a  constant  growth  rate,  A(y):^,  this  evaluates 
loA'p  (writing  Hv/V  s).  With  an  exponential  growth  rate,  A(y)=^exp  Jfljd-y)  (y  is  defined 
backwards  along  the  channel:  A(x)  ^AexpJTx;  \-d->  where  d  is  the  length  of  the  channel). 
This  evaluates  tn^  „ 

o  •  C.3) 

-JL**j>V4/(y(i-X)+  f). 


xc  < 


(i)  y.y':  charge  deposited 
is  A(d-y')dy 


(ii)  y*0  streamer  arrives 
at  time 
t*y'/vs 
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Equating  times 

y'/vs •  (y-y *)/ve 

or 

y’-(i-v/vs)y’ 

Fig. 5.7 


(iii) 


y:  electrons  drift 
to  y  in  time 
(y-y')/ve=  t 


(i) 


charge 

A(d-y ' )dy  e  1  ^  y  ^ 


Cii) 


> 


Fig. 5. 8 


Variation  of  positive  ion 
number  Ns  into  field  E 
along  the  streamer  path 
(after  Gallimherti) 


Fig.  5. 16 

Showing  a  small  dark 
current  between  the 
two  corona  events. 
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Using  n  powrr  function  in  not  us  troctnulc  u  pmhlrn.,  substituting 
into  Erj.5.2  'x  =  (.'->',  us  oScvc)  and  integrating  by  parts  yi e ! c.s: 

flD  oo 

A  -  -M (d- **{>(- 

v.ritmy  pL(d -0-\)'f)n*y>(-fy)dy 

the  above  gives  the  rerurrrr.ee  relation: 

V  -  M*- 
*  £ 

If  U  ot A(y'sC;,  the  size  of  the  streamer  tip  charge,  then  this  implies  Un«Cdn.  Clearly 
Cq.  5.4  is  inconsistent,  with  this  oven  for  n=l,  as  then  l  J ,  =  ~-jr~  In  SI  units 

,  fi  t1 

typical  values  are:  d-vO.lm  and  jl*. 2Qn  .(l-V/vs)vG.9,  and  so  these  terms  are  of  nearly 
equal  size.  Uj  is  not  even  approximately  «e  o’,  as  would  be  required  if  the  tip  charge 
were  to  be  proportional  to  the  channel  charge. 

So  the  maximum  nap  charge  and  the  system  size  prior  to  impact  will  be  propor¬ 
tional  if  a  constant,  or  an  exponential  growth  is  assumed,  but  not  if  a  power-law  is  used 
for  the  streamer  function,  A(y).  The  almost  exponential  growth  in  Q^(n.ax)  with  distance 
suggests  an  exponential  growth  function,  A(y)  as  the  two  would  then  be  self-consistent. 
This  analysis  depends  on  the  channel  having  decayed  to  zero  free  charge  at  its  base  by 
the  time  the  cathode  is  reached.  This  is  clearly  not  the  ease  for  the  5. Jan  systems 
studied  (see  e.g.  Fig. 5. 5',  and  so  this  analysis  breaks  down  at  5. 3cm.  This  will  explain  the 
fact  that  the  5.3  data  do  not  coincide  with  the  rest  at  fields  slightly  above  stability  field 
and  will  also  explain  the  apparently  steeper  growth  of  the  spatial  O  (max)  curves 
(Fig. 5. 6)  in  the  early  stages. 

The  almost  exponential  spatial  growth  (Fig. 5.0)  underscores  the  difference  be¬ 
tween  point -plane  and  uniform  field  corona.  In  this  study,  the  streamer  growth  seems  to 


n(/-V/vQUv-' 
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be  exponential  v/ith  distance,  implying  for  constant  velocity  a  growth  rote  proportional 
to  the  system  size.  This  can  bo  compared  for  example  with  the  calculated  t ip  size  ob¬ 
tained  hy  C.allimberti  (1972c)  (Fig. 5. 3)  which  rises  and  then  falls  in  the  highly  non- 
uniform  rod-plane  field. 


5.2.4  Charge  Injection 

The  charge  injecLed  by  secondary  processes  when  the  streamer  reaches  the  cath¬ 
ode  results,  as  discussed  above,  in  a  large,  narrow  spike  iri  the  %G  curves.  This  appears  to 
be  a  well-defined  parameter  and  its  size  has  been  rr, ensured.  Put  a  slight  v.arning  ought 
to  be  attached  to  this  data,  the  records  are  digital  and  the  peak  in  Sg  is  very  narrow  -  a 
few  samples.  So  whilst  it  is  fair  to  say  that  only  one  experimental  parameter  (i^)  is 
varying  at  all  rapidly  when  Sg  is  found,  the  error  in  the  maximum  may  still  be  fairly 
large. 


Fig. 5. 10  shows  the  variation  of  this  parameter  with  pressure,  E/p  and  propagation 


shows  how  strongly  both  depend  on  the  system  growth.  The  magnitude  of  the  charge  in¬ 
jection  rate  is  very  high,  •%.  amps,  hut  of  a  short  duration,  ~10's  of  ns.  A  ]  amp  pulse  of 
current  lasting  (typically)  32  to  40nS  will  inject  32  to  4GnC  of  negative  charge  into  the 
gap. 

Even  ICO  streamer  tips  will  have  a  total  charge  of  perhaps  InC  (10  '  ^  coul.  each). 
Thus  it  can  be  seen  that  the  “lec'ron  injection  of  somc*x.50nC  from  the  cathode  is  very 
efficient  by  comparison  v/ith  the  streamer  active  region.  This  is  not  very  surprising  if  it 


is  considered  that  the  ionization  potential  of  oxygen  is—J6eV,  compared  with  a  work 
function  of  ~4.2eV  for  aluminium,  and  that  the  collision  frequency  for  production  of  a 
photo-electron  from  the  cathode  is  many  orders  of  magnitude  higher  than  that  for  ion¬ 
ization  of  a  gas  molecule. 

In  fact,  the  removal  of  electrons  from  the  base  of  the  channel  will  result  in  a  net 
positive  charge  on  the  channel,  and  integrating  the  conduction  current  yields  a  value 
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of  ~  2G  tn  3GnC  at  most  for  1 1  ie  conduction  charge.  The  injection  of  electrons  fro r  Iff 
cathode  then  results  in  a  change  froni  a  net  positive  to  a  nei  negative  charge  on  tl no 
ctianr.ol. 

In  summary,  the  charge  in  the  streamer  active  reijion  n  ay  be  a  little  below  ~lnC, 

liut  results  in-~50nC  heing  injected  from  the  cathode  due  to  the  greater  photo- 

efficiencv.  The  free  charge  in  the  nap,  G,  ,  is  —500  to  lOOOnC  in  extreme  cases  but  the 

net  positive  charge  on  the  channel  is  only~20  tn  30nC  at  most,  and  the  injection  from 

tire  cathode  reverses  the  polarity  of  this  net  charge.  The  charge  injection  from  the 

cathode  is  so  small  by  comparison  with  G  that  the  value  of  C  'max)  will  be  but  little 

g  g 

altered  by  this.  This  explains  the  difficulty  observing  the  injection  of  charge  from  G^ 
records  except  at  low  E/p. 


5.2.5.  Attachment  Frequency 

An  exponential  fit  routine  was  used  to  analyse  the  data.  The  end-points  could 
be  varied  in  real-time,  and  a  fit  obtained  over  any  desired  portion  of  tbe  record. 

Whether  the  data  actually  follows  an  exponential  form  is  immaterial  to  the  fit  routine, 
which  happily  calculates  a  decay  or  growth  frequency  for  the  data  selected.  The  routine, 
though,  does  find  an  accuracy  figure,  following  the  method  outlined  by  Darford  (1967). 
The  appropriateness  of  the  calculated  values,  therefore,  depends  to  no  small  amount  on 
the  portion  of  the  curve  selected,  but  is  represented  by  the  calculated  error. 

The  selection  of  data  is  important  because  the  derivation  of  the  attachment  fre¬ 
quency  requires  there  to  be  no  generation  of  free  charge  in  the  gap  and  that  the  electron 
drift  speed  be  known.  The  data,  plotted  in  Fig.  5.1 1,  were  not  calculated  after  any 
secondary  phenomena  (like  regenerations),  for  the  possible  channel  heating  might  affect 
the  local  pressure,  and  hence  the  electron  drift  speed.  At  low  pressure,  where  there  is 
apparently  a  marked  drop  in  attachment  early  in  the  mature  stage  of  the  discharge  re¬ 
sulting  from  significant  removal  of  charge  from  the  gap  at  the  needle-end,  the  end  por¬ 
tion  of  the  Gp  record,  where  the  conduction  current  hasn  ost  fully  decayed,  is  preferred 
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for  the  calculation  of  attachment  frequency.  In  cases  w here  the  channel  was  regener¬ 
ated  this  is  not  possible,  for  a  secondary  process  has  intervened,  hut  these  data  might  he 
intercom parable  and  so  have  been  included.  The  attaclmierit  frequency  sought  is  a  lower 
bound  on  the  low  pressure  values. 

In  general,  the  calculation  was  extended  to  as  much  of  the  discharge  as  possible 
whilst  maintaining  a  satisfactory  fit.  The  selection  of  the  portion  of  the  decay  for  the 
fit  was  not  normally  a  problem  since  usually  the  decay  in  Gn  was  satisfactorily  exponen¬ 
tial  over  the  whole  of  the  record. 

The  low  E/p  data  at  6CC  torr  warrant  discussion.  Approaching  the  stability  field, 
it  would  seem  that  the  attachment  rate  falls  steeply  as  the  reduced  field  approaches  this 
value,  at  v.’hich  the  attachment  rises  abruptly  to  a  very  high  level.  The  key  is  that  Ihis 
occurs  at  the  stability  field  -  at  which  the  channel  is  only  just  sustaining  itself  and  so 
below  this  decay  of  G^  sets  in  before  the  channel  is  dead.  The  generation  of  charge  is 
not  zero  and  a  true  attachment  rate  is  not  found.  This  data  has  been  included  in  the 
figure  to  illustrate  this  process. 

Ttie  data  have  been  converted  from  attachment  frequencies  into  attachment  co¬ 
efficients  as  this  is  still  the  more  widely  used  parameter.  F urther,  the  reduced  attach¬ 
ment  coefficient  is  plotted,  for  this  removes  the  inverse  pressure  dependence  resulting 
from  the  greater  mean  free  path.  If  a  purely  two-body  attachment  process  operates, 
then  all  the  curves  should  coincide  if  *^/p  is  plotted  against  E/p. 

This  is  almost  the  case,  as  it  con  be  seen  from  Fig. 5. 11  that  *^/p  is  only  lov.'er  at 
greatly  reduced  pressure.  It  would  seem  that  no  clear  distinction  can  be  marie  from  this 
data  as  to  the  nature  of  the  attachment  process;  it  clearly  involves  mainly  two  body  pro¬ 
cesses  but  is  not  exclusively  so. 

Tire  range  of  attachment  coefficients  measured  corresponds  to  those  found  by 
other  workers,  as  does  the  scatter  of  experimental  points,  Figure  5.12(a)  showing  a  com¬ 
pilation  of  reduced  attachment  coefficient  measurements  is  reproduced  here,  with  the 
range  of  observations  covered  by  the  scales  in  Fig.5.12  indicated  by  the  small  rectangle. 


Reduced  attachment  coefficient  ^/p  (m.torr)'^ 


87 


89 


The  published  data  are  evidently  scarce,  but  the  range  of  n  casuren.ent  is  within  the 
range  nf  measurement  reported  by  other  authors. 

8.2.6  Charge  Injection  Rate 

As  discussed  above  (section  3),  applying  charge  conscrvtion  to  the  gap  yields  the 
parameter  So,  the  charge  generation,  or  injection  rate.  Since  the  data  v/as  very  noisy 
and  almost  any  rapidly  growing  function  could  be  used  to  characterize  the  data  v/ith 
equally  good  (had)  accuracy,  for  simplicity  and  computational  convenience  the  expo¬ 
nential  fit  routine  used  to  obtain  the  attachment  frequency  v.as  also  used  to  find  a 
growth  frequency, P,  for  the  charge  generation  rate,  defined  from  the  equation: 

St 1--  sJ1*  5.; 

This  may  prove  to  have  been  on  unfortunate  choice,  considering  the  remarks  in 
section  6.1  concerning  the  unsuitability  of  an  exponential  for  the  streamer  growth.  The 
alternative:  a  polynomial,  contains  a  great  many  degrees  of  freedom  and  considerable 
thought  must  go  into  deciding  the  form  of  equation  to  use  in  the  fitting  routine,  if  mean¬ 
ingful  results  are  required.  As  a  first  attempt  at  characterizing  the  behaviour  of  this  en 
tirely  new  parameter,  Sg,  an  exponential  fit  is  on  balance  as  adequate  as  any  other, 
furthermore,  it  has  the  advantage  of  providing  the  instantaneous  growth  rate  directly. 

Since  the  charge  generation  rate  depends  on  the  size  of  the  active  region,  the 
growth  frequency  of  this  parameter  also  represents  the  temporal  growth  rate  of  the 
streamer  tip  charges. 

Plotting  this  very  secondarily  deriveo'  parameter  (Fig.5.13)  shov:s  no  clear  trend 
with  pressure,  or  system  size,  but  a  broadly  linear  growth  against  C/p. 

Considering  the  strong  influence  of  pressure  on  streamer  velocity  and  streamer 
diffusion  radius,  it  is  by  no  means  clear  that  the  streamer  growth  will  be  independent  of 


pressure. 
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5.3  Pre-Breakdov. n;  A  case  r.twr,'y 


Turning  from  straightforward,  classic  streamers  to  the  various  pre-breakdown 
modes  discussed  in  section  A,  the  data  is  fragmentary,  and  difficult  to  synthesize. 
Accordingly,  to  try  and  clarify  the  relativity  of  the  processes,  a  particular  set  of  experi¬ 


mental  conditions  was  selected  for  closer  study.  All  the  500  torr  data  we  re  selected  for 
the  low  onset  field  at  this  pressure  means  that  low  reduced  field  measurements  are  avail¬ 


able. 


Fig.  5.1£  shows  a  multiple  streamer  event,  obtained  at  just  above  onset.  The 
current  pulses  arc  distinct;  they  do  not  run  into  each  other.  The  repetition  frequency  is 
much  higher  than  is  normal  for  streamers  (usually  kHz).  Records  of  this  character  were 
difficult  to  obtain  as  they  are  seldom  observed  at  low  field. 

The  displacement  current  record  displays  the  tell-tale  character  of  this  type  of 
corona  (Mode  00):  as  the  discharge  progresses  the  pulses  have  successively  larger 
positive-going  peaks  in  the  i^  record,  the  negative-going  ones  being  successively 
smaller.  This,  as  discussed  in  section  3,  indicates  that  the  discharge  is  located  progess- 
ively  closer  to  the  needle-end.  Visually  though,  the  corona  is  similar  to  classic  corona  os 
might  be  expected  since  the  multiple  mode  occurs  with  such  a  low  probability.  The 
current  records,  furthermore,  ore  indisinguishable  from  those  of  the  high  frequency 
modes  (Mode  (iii)),  where  the  streamers  are  again  visually  similar  but  on  the  whole 
brighter  than  classic  streamers.  Again,  this  is  probably  due  to  the  probability  of  the  mul¬ 
tiple  streamer  event  being  greater  in  this  mode. 

To  try  and  characterize  the  behaviour  of  this  mode,  a  "dark-time"  t^  was  defined 
(Fig. 5. 15)  in  analogy  with  long  gap  impulse  studies  (e.g.  C.allimberti,  1979).  Either  the 
transit  of  the  first  streamer,  or  the  negative  peak  in  a  subsequent  corona  is  taken  os  the 
start  of  the  dark  period,  the  end  as  the  beginning  of  fr°sh  current  grow  th.  In  Table  5.1 
these  dark  times,  the  positive  and  negative  i^  peaks,  the  maximum  value  uf  G^  and  the 
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ver/hor  scale 
y  J  A. 09  nC 
mi  2. 00  uS 


integral  of  the  corresponding  conduction  current  pulse  are  given,  for  all  the  corona  of 
this  type  observed  at  5GC  torr. 

The  total  nun, her  of  pulses  observed  is  small,  but  it  becomes  clear  that  for  most 
records,  the  second  corona  has  a  larger  conduction  charge  (G^  =  Ti ^dt)  than  the  first.  The 
dark  time  is  less  between  the  second  and  third,  than  between  the  first  and  second 
coronae.  Event  number  87.9,  however,  has  a  smaller  second  corona  conduction  charge, 
Qc,  and  has  a  longer  dark  time  after  it.  An  inverse  relationship  then  appears  to  hold  in 
nearly  all  cases:  other  factors  being  constant,  the  dark  time  is  reduced  and  the  con¬ 
duction  charge  increases  as  the  discharge  progresses. 

In  general, also,  the  electron  charge,  G^,  decreases  with  each  subseguent  corona 
event  (the  exceptions  being  events  85.10  and  438.0,  which  have  the  longest  dark  times, 
during  which  it  is  possible  that  the  cathode  has  recharged  a  little  more,  permitting  a 
stronger  growth  in  the  subsequent  corona}. 

A  further,  interesting  observation  is  that,  there  is  a  "dark  current"  in  the  con¬ 
duction  current  records,  of  between  180  and  280jiA  (Fig. 5. 16),  though  this  feature  is  not 
always  present.  Not  enough  data  exists  to  make  out  any  clear  trend. 

5.3.2  High  reduced  field/Mode  (i) 

At  higl  sr  reduced  field,  the  prebreakdown  current  records  (Fig. 5. 17)  look  com¬ 
pletely  different.  The  second  corona  occurs  much  sooner  than  at  lower  E/p  and  may  re¬ 
peat  on  very  short  timescales.  The  events  run  into  each  other,  rather  than  being  distinct, 
as  observed  at  low  Field.  And  it  will  be  recalled,  the  discharge  is  concentrated  into  a 
single  channel. 

The  same  parameters  as  for  the  multiple  mode  are  shown  in  Table  5.2.  Defining  a 
second  dark  time,  t^o,  is  difficult  for  this  data,  as  the  end  of  the  second  corona  is  dif¬ 
ficult  to  make  out.  Looking  at  Table  5.2  one  first  notes  that  the  dark  time,  t^j  is  on 
the  whole  shorter  than  at  the  lower  reduced  field  in  Table  5.1.  Between  the  lines  of 
Table  5.2  (excepting  764.3),  the  trend  is  the  other  way,  the  dark  period  increasing 
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slightly  with  increasing  reduced  field.  !’ut  this  may  la  ji.it  -j  reflection  uf  the  following: 

It  was  suggested  in  section  5.3.!  above  that  an  inverse  relationship  exists  between 
Gc  and  t^.  In  ttiis  case  the  conduction  charges  measurer'  increase  very  roughly  v.ith.  C/p 
and  this  may  contribute  to  reversing  the  gross  trend  tov.ard  lower  t  at  higher  r/p 
which,  it  must  again  be  stressed,  is  large  and  unambiguous.  The  influence  of  Gc  is  very 
much  a  secondary  effect  by  comparison. 

No  clear  influence  of  electron  charge  on  dark-time  is  evident.  One  notes,  how¬ 
ever,  that  C!^  is  generally  reduced  for  each  subsequent  corona,  while  CC^  increases.  List¬ 
ing  the  reference  numbers  of  the  records  in  order  of  increasing  G  and  then  in  order  of 

increasing  t  „  one  gets  (  764,  4 ! 9,  41G,  366!  for  C  and  (  -IG,  764,  419,  366!  for 
ok  g 

t(l,  .  The  electron  charge  is  a  m. ensure  cf  the  size  of  the  streamer  system,  and  this  lack  of 
any  correlation  with  T  indicates  that  the  second  corona  inception  is  independent  of  t’r 
scale  of  the  streamer  system. 

In  fact,  the  suggestion  from  all  this  data  is  that  the  gross  dependence  of  the  car- 
time  on  the  applied  field  arises  because  the  mechanism  responsible  for  the  format  ,u:  1 

the  second  corona  is  the  drift  of  space-charge  accun  ulnted  near  the  needle  tip.  ,v.  I- 
will  clear  faster  in  higher  applied  reduced  Lnpiacinn  fields.  This  else  explains  cv 
sccondary  trend  towar d  longer  dark  time  at  higher  conduction  charge  in  t*'c  1  , ‘ 
data  set;  there  is  more  charge  to  clear  from  the  point. 

flut  at  low  field,  successive  eorena  pulses  occur  ■■  d*  lr.v«  t  ..  . 

size  of  the  conduction  charge.  This  is  suggestive  of  the  _,v  en 
states.  Being  current  dependent,  the  greater  density  ,.f  ■.  br.  t  • 

conduction  charges  v/ill  increase  tt  e  e  •  fe-t  i  v>-  mr  . 

corona  possible  at  ari  earlier  stage  m  the  on  ft  ■  f  t*  .  . 

The  difference  between,  fhi  visual  . 

being  iriucti  like  a  higfi  frequency  .  lass.'  r: 
bright  cftannel,  extending  see  e  ■■.  ,  -  *  . 

channel  when  the  regeneration  .  /  .  ' 
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that  the  channel  is  still  active.  The  effect  of  this  may  be  to  guide  the  second  ccronn 
axially  along  the  highest  conductivity  path;  along  the  main  branch  of  the  streamer 
system,  that  is.  That  this  channel  is  so  very  intense  makes  sense  considering  the  exis¬ 
tence  of  free  electrons  along  it:  when  the  streamer  arrives  at  a  point  in  the  channel,  it 
does  not  have  to  generate  secondary  electrons,  there  are  plenty  there  already.  The 
streamer  growth  would  he  that  much  more  vigorous  as  a  consequence. 

A  question  remains:  if  the  space  charge  can  clear  from  the  point  in  time  to  pro¬ 
duce  corona  only  a  few  ps  after  the  primary  streamer,  why  does  this  happen  only  occas¬ 
ionally,  the  normal  interstreamer  time  being  of  the  order  of  ms?  The  suggestion  is  that 
the  decay  by  lateral  diffusion,  or  whatever  of  the  vibrational  states  decreases  the  effec¬ 
tive  ionization  coefficient,  to  such  a  level  that  the  space-charge  has  much  further  to 
travel  away  from  the  point  before  ionization  is  again  possible  there. 

It  may  also  he  argued  that  the  streamer  channel,  if  still  active,  may  generate  a 
more  efficient  source  of  starter  electrons  for  the  discharge  than  would  normally  he 
available  to  an  avalanche  developing  in  virgin  air  (by  photodetaclmient,  or  by  detach¬ 
ment  from  negative  ions  which  would  otherwise  diffuse  av/ay).  Hence  some  critical  time 
exists,  beyond  which  second  corona  cannot  form,  and  the  usual  ion  clearing  time  («.ms) 
determines  the  corona  frequency. 

5.3.3  Growth  and  N’agnjtude  of  Second  Corona 
As  already  shown,  the  electron  charge,  C^,  is  very  strongly  related  to  the  scale  of 
the  discharge.  A  question  being  addressed  is  how  the  discharge  scale  influences  break¬ 
down  and  streamer  growth.  Though  no  two  measurements  in  this  data  set  has  been  made 
at  comparable  reduced  field,  looking  at  the  higher  reduced  field  data  reveals  a  very  clear 


trend. 
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Taking  the  second  corona  conduction  charge  ms  a  ii  ecsure  of  the  size  of  the 
second  corona,  and  considering  the  influence  of  electron  charge  from  the  first  corona,  it 
seems  the  larger  the  electron  charge,  ,,  the  larger  tt.e  serond  corona  conduction 
charge,  G  o>  whatever  tfie  reduced  fields.  Looking  at  the  current  records  visually  con¬ 
firms  this  general  behaviour  (Fig. 5.  IE).  The  records  in  longer  gaps  are,  on  the  whole, 
more  energetic,  the  second  corona  growing  more  vigorously  than  in  shorter  gaps.  This  is 
a  very  distinct  trend,  and  the  influence  of  G)  is  very  large  indeed,  the  second  enrnna 
conduction  charge  can  range  between  and  2DnC  as  a  result  of  the  enhanced  C  ,  in 
fact  ranging  from  being  larger  or  smaller  than  the  first  conduction  charge  as  n  result. 

Sun  mi  ary 

The  trends  in  the  data  studied  suggest  that  'a}  the  scale  cf  the  discharge  has  little 
influence  on  the  second  corona  inception;  (b)  the  dark  time  before  the  second  corona  in¬ 
ception  is  reduced  as  E/p  increases,  indicating  that  the  clearing  of  the  ions  from  the 
needle-tip  region  controfs  the  second  corona  inception;  (c)  since  a  large  conduction 
charge  can  reduce  the  dark  time,  it  is  suggested  that  vibrational  states  createu  in  the 
primary  discharge  will  increase  the  effective  ionization  coefficient,  leading  to  quicker 
second  corona  inception.  This  is  advanced  as  the  controlling  process  -  if  no  second 
corona  has  been  produced  before  the  vibrational  states  diffuse  out  of  the  channel  (after 
30ps),  the  ions  have  to  clear  in  the  usual  ms  time  interval.  Finally  the  difference  be¬ 
tween  Mode  (i)  and  Modes  (ii)  and  (iii)  is  suggested  to  be  due  10  the  channel  still  being 
active  for  mode  (i).  With  a  non-zero  residual  free  electron  density  in  the  channel,  a 
second  streamer  would  have  a  reservoir  of  starter  electrons,  would  grow  more  strongly 
than  usual  and,  furthermore,  would  lie  guided  down  the  old  channel,  producing  n  single 
bright  unbranchcd  discharge. 
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Section  6.  Discussion 


6.1  Propagation  in  uniform  field 

6.1.1  Conduction  Current  Characteristics 
To  explain  the  time-to-peak  and  peak  value  of  the  streamer  conduction  current  is 
a  more  interesting  problem  than  it  would  at  first  seem.  Differentiating  Eq.  3.15  for  the 
conduction  current  yields: 


dij*)  .  vV7tVf  4  A(tV)-  »7^V)l 

TT  *  Wtv)  (  J 


Eq.6.1 


6.1.1, a  Problems  with  an  Exponential  Model 
2Tx 

If  A(x)  is  an  exponential,  A(x)  =  Be  ,  dA(tV)/e(  (tV)  =  3TA(tV),  and  solution  of 
Eq.6.1  to  find  the  maximum  of  i^  requires  that 3=^ .  In  other  words,  an  exponentially 
growing  generation  function  A(x)  would  lead  to  the  current  surviving  to  reach  the  point 
always  either  outstripping  the  attachment,  or  be  overtaken  by  it,  according  to  whether? 
is  >  or  <  respectively.  This  is  apparently  in  conflict  with  the  choice  made  earlier  to 
fit  an  exponential  to  the  charge  generation  rate  curves,  $Q  and  suggests  that  the  expo¬ 
nential  growth  rate  assumed  throughout  the  discussion  of  $Q  cannot  be  correct. 

In  fact  the  measured  (assumed  exponential)  growth  frequency  of  £Q  is  —5  to 
10MHz,  greater  than  the  attachment  frequency,  on  the  whole.  Consequently,  run-away 
growth  in  ic  might  be  indicated.  But  note  that  (Eq.3.16),  iQ(t)=A(v#t)vs.  The  measured 
values,  P,  of  the  growth  frequency  of  SQ  are  obtained  by  fitting  an  exponential  relation¬ 
ship,  iQ(t)=BeTt  (Eq.  5.5)  to  the  data.  Substituting  in  Eq.  3.16,  and  putting  x=v#t,  yields: 


AW  -  (‘VirJeT'/'i 


Eq.6.2 


which  corresponds  with  the  assuned  exponential  growth  in  A(x):  A(xKAe  (section  5.5.6) 
given  that.A.=B/v#  and  Y=Vh9.  Hence  J",  the  (spatial)  growth  coefficient  for  A(x)  (~ 
25m ’1)  is  below  the  attachment  coefficients  (»j-*10Chn’*)  measured  at  the  same  time,  and 
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the  tail  of  the  conduction  current  will  be  predicted  to  be  a  strongly  decaying  current 
regime,  as  is  indeed  observed. 

6.1. l.b  Timc-to-Peak/Power-Law  Model 
The  behaviour  of  i^  as  defined  by  Eq.3.15,  can  easily  be  made  to  peak  by  using 
any  function  growing  more  slowly  than  exponentially.  For  example,  as  discussed  by 
Bicknell  <fc  Shelton  (1986),  a  simple  power  law  of  the  form: 

A(x)  «■  A*-*'  Eq.6.3 

yields  a  time-to-peak,  t  ^  =  n/tjV  and  the  value  of  this, 

i  („) .  Av(v)V’VV' 

*•  ^  Eq.6.4 

Though  this  still  cannot  explain  the  clear  levelling-out  of  ic(max)  as  E/p  increases 
(Fig.5.2,  section  5.2.1),  for  it  is  clear  that  if  A(x)  grows  more  strongly,  that  is  if  Aor  n 
increase,  then  ic(max)  will  increase  also,  a  prediction  of  this  simple  model  completely  at 
odds  with  observations.  (Though  at  least  a  peak  of  some  kind  is  indicated;  an  improve¬ 
ment  over  the  exponential  model.) 


6.1. l.c  Timescales 

The  behaviour  of  the  peak  in  the  conduction  current  is  consdierably  complicated 
by  the  timescales  involved.  The  peak  occurs  within  —200ns  at  most  and  ignoring  the 
more  rapid  drift  in  the  tip  region,  this  corresponds  to  a  distance  of  -8mm  at  electron 
speeds  of  40km/s.  Hence  the  portion  of  channel  which  contributes  to  the  peak  is  close 
enough  to  the  needle-end  for  the  conditions  there  to  determine  the  magnitude  of  the 
peak. 

The  proposition  that  it  is  the  larger  scale  of  the  streamer  systems  which  is  res¬ 
ponsible  for  the  dependence  of  the  conduction  current  maximum  on  the  electrode  spacing 
is  completely  ruled  out  by  the  foregoing,  as  by  the  time  the  conduction  current  has 
peaked,  the  streamer  will,  on  the  whole,  have  advanced  just  4cm  (at  200len/s)  and  the  re¬ 
maining  portion  of  the  channel,  not  having  been  formed  cannot  possibly  contribute  to  the 
size  of  the  peak  in  the  conduction  current.  This  ignores,  of  course,  low  pressure  stream- 
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art,  where  the  low  attachment  rate  causes  the  peak  to  occur  after  the  streamer  transit. 
Even  than,  for  electrons  to  bevel  the  shortest  channel  length  of  5.3cm  requires  ~1.3ps. 

Whatever  the  pressure,  from  the  above  discussion  it  can  be  categorically  stated 
that  a  mature  channel  does  not  bear  on  the  conduction  current  peak,  just  its  lower  por¬ 
tion,  and  generally  only  the  first  cm  or  so  is  involved.  In  thundercloud  streamer  systems, 
this  argument  is  more  forceful:  the  part  of  the  channel  deciding  the  peak  in  ic  is  de¬ 
coupled  from  the  remainder  of  the  streamer  system. 

6.1.1.d  Laplacian  Field 

How  then  can  there  possibly  be  the  obseved  dependence  of  i^tmax)  on  propagation 
distance  reported  in  section  5.2.1?  Or  even  a  peak  of  any  kind,  since  the  streamer 
growth  seems  to  be  at  least  roughly  exponential  (section  5.2.3)?  Worse  still  the  simple 
power-law  model,  Eq.6.3  for  A(x)  implies  that  ic(max),  whilst  at  least  defined,  will  not 
depend  on  the  scale  of  the  discharge  since  this  is  not  taken  into  account,  nor  can  it  ex¬ 
plain  the  levelling-out  as  E/p  increases. 

Returning  to  the  variation  of  the  maximum  in  ic  with  electrode  spacing,  this  has 
already  been  shown  not  to  be  a  result  of  the  longer  propagation  in  the  longer  inter-elec¬ 
trode  gaps,  and  it  is  inferred  that  the  Laplacian  field  distribution  is  somehow  different 
near  the  point  in  such  situations. 

The  les  Renardieres  Groupe  (Resuitats  de  1973)  have  calculated  point-plane  per- 
unit  fields  for  a  5m  and  10m  electrode  spacing  with  the  same  conical  point  (Fig.6.1).  The 
general  conclusion  is  that  the  enhanced  field  region  extends  further  out  from  the  point  in 
the  small-gap  case.  Hence  it  can  at  least  be  inferred  that  the  small  gap  stramers  get 
off  to  a  better  start  and  on  account  of  the  extended  tip-region  will  travel  into  the  re¬ 
mainder  of  the  gap  with  a  larger  initial  tip  potential.  The  anallness  of  the  region  of 
channel  contributing  to  the  peak  in  ic  means  that  effects  such  as  this,  operative  near  the 
needle-end,  can  have  a  large  influence  on  the  value  of  the  ptlak.  This,  at  least  in  prin¬ 
ciple,  explains  the  low  gap-length  arm  of  the  minimum  in  Fig .5,1. 

In  fact,  the  measurements  of  Sadik(1980)  show  that  streamers  launched  from  a 
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needle  at  a  few  kV  above  the  anode  potential  grow  initially  as  they  propagate  into  the 
gap.  In  this  case,  the  needle  is  at  earth  potential,  but  the  ambient  potential  in  the  gap  at 
the  location  of  the  needle-end  were  it  not  there  would  have  been~En,  where  n  is  the 
height  of  the  needle-end.  This  will  drop  away  in  a  distance  (see  the  le  Renardieres 
Groupe  results)  which  increases  as  the  electrode  spacing  falls.  So  the  potential  gradient 
at  the  end  of  the  needle  is  broadly  independent  of  the  gap  spacing,  but  falls  off  more 
slowly  the  closer  the  electrodes  are.  This  is  the  effect  sought.  The  enhanced  growth 
initially  will  lead  to  an  increasing  conduction  current,  falling  thereafter  as  attachment 
takes  over  from  the  (comparatively)  weaker  growth  further  into  the  gap  once  the  needle 
has  been  left  behind.  By  then  the  maximum  has  been  determined,  as  the  first  cm  or  so 
only  of  the  channel  contributed  to  the  maximcm  and  the  maximum  will  not  depend  on  the 
applied  reduced  field. 

The  growth  observed  in  conduction  current  peak  as  gap  length  increases  cannot 
continue  indefinitely  and  so,  given  the  targe  spread  in  the  data,  this  is  ignored  as  not 
being  significant.  Aftidealized  distribution  would  decrease  from  some  high  value 
reaching  a  constant  level  as  the  gap  passed  some  critical  length. 

6.1.2  Streamer  Growth 

The  investigation  of  streamer  growth,  which  was  the  main  object  of  this  study,  has 
most  often  been  undertaken  in  non-uniform  point-plane  geometry.  So  the  slant  of  this 
discussion  must  necessarily  be  rather  different  from  that  normally  encountered  in 
streamer  physics  literature. 

Studies  in  point-plane  geometry  have  tended  to  concentrate  on  the  range  of  the 
streamers,  for  the  monotonic,  steep  fail  in  the  Lsplacian  field  leads  to  streamers  dying 
throughout  their  life  under  these  conditions.  By  contrast  the  streamers  can  be  expected 
to  grow  at  a  uniform  rate  in  the  (quasi -)uni form  applied  fields  used  in  this  investigation. 
This  is  dramatically  exemplified  by  the  experimental  electron  charge  curves  (Fig.5.5), 
which  show  that  increases  by  an  order  of  magnitude  as  a  function  of  the  discharge 
scale,  indeed  increasing  exponentially  with  streamer  propagation  length. 
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6.1. 2.a  Growth  Frequency, 

In  an  attempt  to  characterize  the  growth  of  the  streamer  system,  the  para¬ 
meter  &Q  was  defined  (section  3.5.1),  the  charge  generation  or  injection  rate  and  it's 
growth  frequency,  ,  measured  (section5.2.6,Fig.5.13).  Conceptually  this  is  a  little  diffi¬ 
cult  to  visualize  at  first,  but  consider  that  the  charge  injection  (or  generation)  rate  is 
just  the  charge  generated  in  the  discharge  (allowing  for  attachment)  needed  to  satisfy 
charge  conservation  per  unit  time.  It  is  hence  representative  of,  if  not  actually  equal  to 
the  current  in  the  streamer  active  region.  Thus  the  growth  frequency,  ,  of  this  para¬ 
meter  can  be  thought  of  as  the  growth  frequency  of  the  charge  in  the  streamer  tip 
region.  This  is  what  has  beren  measured. 

It  will  be  recalled  that  F  was  found,  within  a  large  experimental  spread,  to  be 
roughly  linearly  dependent  on  E/p,  and  independent  of  pressure  and  gap  length.  Actually 
a  better  formulation  for  F  is  a  square-root  function,  for  it  must  fall-off  to  zero  at  the 
stability  field.  Undue  significance  should  not  be  attached  to  the  particular  form  given: 
the  spread  of  the  data  is  such  that  many  different  formulations  may  be  compatible  with 
the  measurements. 


6.1.2  b  Spatial  Growth  of  Qgfrnax) 

Plotting  the  spatial  growth  (Fig.5.6)  of  Qgfmax)  fleshes  out  a  problem  already 
implied  by  the  above:  how  can  the  spatial  growth  apparently  not  depend  strongly  on  the 
applied  reduced  field,  when  the  growth  rate  (T1)  increases  broadly  linearly  with  the  same 
parameter?  Especially  considering  that  the  possibility  of  a  charge  or  energy  availability 
mechanism  limiting  Qg&nax)  is  thought  (section  5.2.3)  to  be  incompatible  with  the  pheno¬ 
mena,  and  that  the  cathodic  injection  and  charge  along  the  channel  are  both  thought  to 
be  proportional  to  the  system  size  (in  the  case  of  exponential  growth),  and  so  high  charge 
injection  from  the  cathode  cannot  be  invoked  to  explain  the  phenomena  either  -  a  weakly 
growing  streamer  system  will  have  leas  charge  injection  from  the  cathode  and  less  charge 
residing  along  the  channel:  the  process  cannot  be  separated.  In  any  case,  the  charge  in- 
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jection  is  mall  compared  with  the  free  charge  residing  on  the  channel. 

The  weak  spatial  growth  in  Qg(max),  then,  is  not  an  experimental  construct  but  a 
genuine  property  of  the  discharge.  In  attempting  to  show  how  this  can  be  consistent  with 
an  exponentially  growing  streamer  function,  A(x),  it  is  advantageous  to  collect  together 
the  results  derived  earlier  for  the  free  charge  in  the  gap. 


6.1.2.C  Channel  Model;  Summary 

A  channel  function,  A(x,t)  and  streamer  function,  A(x)  were  introduced,  A(x,t)  being 
the  instantaneous  linear  charge  density  along  the  channel,  A(x)  being  the  charge  density 
added  by  the  streamer  active  region  at  a  particular  positon,  x=v  t.  They  are  related,  in 
a  coordinate  system  where  the  charge  produced  at  X  drifts  to  x  at  time  t,  by  the 
relation:(Fig.5.7) 

w  .v  * x)  V-  /».  x 

(3.17) 


X(*.t)  •  A(>0e',0M- 


The  charge  on  the  channel  as  a  whole  can  be  found  by  integrating  Eq.3.17  for 


specific  forms  of  A(x).  In  the  case  of  an  exponential,  A(x)=_A.e 
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This  is  plotted  in  Fig.3.14  for  typical  values  of  the  parameters,  providing  at  least 
qualitative  agreement  with  the  phenomena.  In  the  case  where  the  current  in  the  channel 
base  has  decayed  to  zero,  Eq.3.18  can  be  integrated  backwards  along  the  channel  to  «, 
yielding  the  free  charge  in  the  gap  in  the  limit  where  the  conduction  current  is  zero: 


Qg  =  I" A(('- Y)y .  (4.2) 
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This  evaluates  to  the  followinq  for  different  forms  of  A(x): 
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6.1. 2.d  Experimental  values  of  Parameters 
In  trying  to  relate  these  formulations  to  the  experimentally  determined  para¬ 
meters,  one  first  notes  that  the  growth  frequency,  P ,  of  <$Q  defined  by 

BePt  (5.5) 

can  be  related  to  A(x),  according  to  the  discussion  above: 

A(*.).{e/„yx/Vs 

(6.2) 

Experimentally,  the  observed  variation  of  the  parameters  can  be  represented  by: 

P=  (MHz), E/p  in  V/cm.torr  (6.5) 


-j£=  <m-t°rr)*1 

ij;  =  a*  +  a7|X-qgj  (km/s),E/P  i 


,  E/p  in  V/cm.torr 


in  V/cm.torr 


)JC  «  (km /s),  E/p  in  V/cm.torr 


(6.6) 


(6.7a) 


(6.7b) 


Where  8^=4  MHz(V/ cm. torr)"  ^  .  B2=4.5V/cm.torr,  aj=0.3(m.torr)"^,  84=2x10'^  jy 
8j=6V/cm.torr,  8g=10(l<rn/s),  8^=0.715  and  in  Eq  6.7  the  values  of  the  coefficients  are  as 


follows: 


p 

torr 

a6 

km/s 

a7 

km/s(V/cm.torr)’^ 

a8 

V/ cm.  torr 

600 

90 

26.3 

4 

500 

95 

30.0 

4 

400 

104 

38.3 

4 

These  are  only  a  crude  nimerical  fit  to  the  data  and  are  shown  in  Fig.6.2  super¬ 
imposed  on  the  appropriate  graphs.  No  particular  claims  are  teingmade  about  the  part¬ 
icular  form  of  Eq.6.5  to  6.7  above,  the  intention  is  to  provide  a  simple  characterization 
of  the  gross  variability  of  the  experimental  parameters. 


6,1.2.e  Predictions  of  Channel  Model 

To  see  if  the  simple  model  presented  above  can  explain  the  observed  behaviour  of 

Q  ,  Eq.  3.18  can  be  rearranged  by  substituting  V=v  v  /(v  +v  )  and  J"  =T/v  ,  yielding 
y  s  0  6  s  s 


-A-  r+714 


(6.8) 


which  is  expressed  entirely  in  terms  of  parameters  measured  experimentally  and  which 
agree  roughly  with  the  range  of  values  found  by  other  experimenter  (except  ve>  which  is 
calculated  from  a  fit  (Gallimberti,  1971a)  to  other  worker's  results).  The  form  of  this 
variation  can  be  seen  to  be  compatible  with  the  stability  field  requiring  T =0  for 
differentiating  with  respect  to  d, 


A  dd 


(6.9) 


and  equating  to  zero  to  find  d _ 

max 


yields: 
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(6.10) 


c^nax  will  be  undefined  for  all  realistic  values  of  vg,*j  and  P  .  furthermore,  substituting 
in  Eq.  6.8  yields  Qg£max)  corresponding  to  this  value  of  dmax  when  the  streamers  are 
below  stability: 


Dependence  on  E/p 

Using  the  values  of  P,  ve,  vg  (and  of  V),  reported  above,  and  plotted  in 
Fig.6.2,  Q  (dVjt  has  been  evaluated  for  various  values  of  E/p  and  d  at  600  torr  and  the 
result  is  plotted  in  F ig.6.3.  The  stability  field  corresponds  to  P =0;  this  is  artificially 
imposed  on  the  nunerical  fit  and  corresponds  to  E/p=4.5V/cm.torr.  It  is  clear  that  the 
interplay  of  the  variation  of  P,  i^and  vg  results  in  a  fairly  flat  growth  in  Qg(max),  except 
near  the  stability  field,  when  the  curves  for  alt  distances  considered,  d,  plunge  downward. 

This  simulates  the  behaviour  of  the  experimentally  determined  plots  of  Qgfrnax),  Fig.5.5 
(reproduced  for  comparison).  The  levelling-off  of  this  plot  has  been  shown  not  to  be  the 
result  of  some  external  parameter  and  is  a  real  property  of  the  discharge. 


Long  Discharges 

This  level  growth  in  Qgfrnax)  is  somewhat  at  odds  with  the  expected  form  at  large 
distances  from  the  anode  (Table  6.1  above).  For  exponential  growth  in  A(x), 

r0- 

and  this  relation  is  plotted  in  Fig.6.4»  It  can  be  seen  that  there  is  indeed  extremely 
strong  growth  in  Qg  as  the  discharge  propagates,  Qg/J[  reaching  ~1051  for  d~2n  and 
E/p  ~8. 5V/ cm  .torr 


Fig.  5.5 
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To  obtain  Qg  requires  an  estimate  of  the  value  of  -A.,  the  charge/unit  length  at 
the  origin.  This  can  be  estimated  from  the  minimum  avalanche  charge  required  for 
streamer  formation,  say  canin’  t09ether  with  the  assumption  that  this  charge  is 
distributed  over  a  channel  length -✓channel  diameter  D  so^-q^^D.  With 
10'^C  (10®)  ions)  and  D-vlO'^m  then.A*-10’^C  m”^.  Using  this  value  Qg  can  be 
estimated  as  a  function  of  reduced  field  and  propagation  distance.  Typical  values  are 
shown  in  Fig.  6.4. 

The  problem  of  extrapolating  a  model  which  works  perfectly  well  in  the  present 
electrode  geometry  to  rather  larger  propagation  distances  is  immediately  apparent  from 
this  data.  A  reduced  field  of  only  5  yields  a  free  charge  of— 10^C  in  only  2m  for 
example  which  is  unrealistically  large.  These  figures  do  however  at  least  emphasise  the 
rapid  growth  to  be  expected  once  the  stability  field  is  exceeded.  It  must  be  remembered 
also  that  this  is  not  a  net  charge  figure  -  there  is  an  equal  and  opposite  positive  charge 
since  the  streamer  channels  contain  a  quasi -neutral  plasma. 

When  assessing  the  growth  of  aircraft  streamers  from  Fig.6.4  it  is  important  to 
include  two  further  factors: 

(1)  The  streamers  are  unlikely  to  propagate  in  fields  much  greater  than  the 
stability  field  for  the  relevant  pressure  since  any  greater  fields  would  most  probably 
already  have  produced  a  breakdown: 

(2)  This  stability  field  is  pressure  and,  therefore,  altitude  dependent  (Fig.4.1). 
Given  that  just  above  the  stability  field  growth  would  be  expected  to  be  significantly  less 
than  that  indicated  by  Fig.6.4  then  the  growth  in  free  charge  would  also  be  much  lower 
than  the  10*®C  calculated  above  in  the  example  given.  More  detailed  measurements 
closer  to  the  stability  field  than  has  been  possible  in  the  present  investigation  are 
required  before  making  valid  growth  predictions. 
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Section  7;  Simmary  and  Conclusion 

The  concept  of  a  positive  corona  streamer  stability  field  such  that  streamer 
growth  at  large  fields  occurs  leads  naturally  to  the  question  of  exactly  how  much  growth 
might  be  expected  for  fields  of  large  extent  and  what  effect  this  growth  might  have  on 
electrical  breakdown.  The  fact  that  large  scale  quasi  -uni  f  orm  fields  (i.e.>10's  of 
metres)  only  occur  in  nature  and  that  the  behaviour  is  pressure  dependent  implies  that 
this  curiously  neglected  corner  Of  discharge  physics  may  be  enormously  relevant  to 
thundercloud  conditions  and  the  aircraft  that  fly  in  their  vicinity. 

This  study  has  been  aimed  at  characterising  the  streamer  growth  in  quasi-uniform 
fields  since  the  large  scale  fields  in  nature  will  necessarily  be  of  this  type  (div  E  =  o / 
and  o,  the  space  charge  density,  is  low).  These  streamers  are  typical  of  those  that  would 
be  generated  at  aircraft  surfaces.  The  3-electrode  arrangement  employed  together  with 
the  analytical  procedure  used,  based  on  the  Schockley-Ramo  theorem,  has  allowed  the 
free  charge  in  the  streamer  system  to  be  determined  as  a  function  of  time  and  for  a 
range  of  fields  and  pressures.  In  addition,  the  rate  of  addition  of  free  charge  -  and  gen¬ 
eration  current  -  has  been  found  for  the  first  time.  The  approach  adopted  also  permits  a 
simultaneous  measurement  of  important  parameters  such  as  attachment  frequency  (co¬ 
efficient)  and  streamer  velocity  together  with  the  growth  frequency  which  charac¬ 
terizes  the  streamer  growth  pattern  in  terms  of  a  single  parameter, P.  Furthermore,  the 
stability  field  for  a  given  pressure  may  be  determined  unambiguously. 

There  is  a  natural  range  of  fields  and  pressures  -  or  more  appropriately  a  range  of 
reduced  fields  -  the  lower  limit  set  by  the  streamer  onset  and  the  upper  limit  by  break¬ 
down.  From  onset  to  the  immediate  prebreakdown  condition,  the  behaviour  is  one  of 
variation  on  a  single  theme.  Although,  from  a  practical  point  of  view  the  breakdown  is 
arguably  the  most  crucial  event  it  is,  of  course,  necessary  to  establish  norm  behaviour  so 
that  departures  from  this  norm  may  be  recognised.  The  data  referred  to  above  (Section 
5)  indeed  quantify  this  norm. 

Prebreakdown  conditions  are  rather  more  difficult  to  quantify  because  of  their 
statistical  nature.  In  all  cases,  however,  the  usual  primary  streamer  event  with  the  asso- 
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dated  decay,  due  to  electrical  attachment,  of  both  conduction  and  displacement  currents 
-  equivalent  to  the  decay  of  the  free  charge  -  is  interrupted  by  the  appearance  of  new 
ionisation.  The  attempt  to  categorise  this  phenomena  is  contained  in  Section  6.  There  is 
nothing  particularly  novel  about  this  observation:  a  similar  pattern  is  observed  in  long 
rod-plane  breakdown  with  the  appearance  of  the  second  corona  in  the  vicinity  of  the 
highly  stressed  electrode.  What  is  unusual  though  is  the  low  reduced  field  at  which  this 
phenomena  is  observed.  A  breakdown  has  been  observed  for  a  reduced  field  of  only 
6V/cm.torr  whilst  examples  of  multiple  corona  are  shown  in  the  text  (e.g.  Fig. 5.14)  for 
similar  values;  this  kind  of  event  has  a  low  probability  but  which  is,  nonetheless,  finite. 

At  a  pressure  of  500  torr  the  breakdown  field  might  than  be  3kv/cm  which  is  within  the 
range  measured  for  a  thundercloud.  An  equally  important  observation  concerns  the 
effect  of  propagation  distance. 

Throughout,  the  effect  of  streamer  propagation  distance  has  been  particularly 
closely  observed  because  of  the  interest  in  extending  any  conclusion  drawn  to  distances 
of  several  metres  so  that  a  valid  comparison  with  aircraft  streamers  might  be  made. 
Generally,  the  effect  of  extending  streamer  length  for  otherwise  similar  conditions  of 
field  and  pressure  makes  little  difference  to  the  measured  parameters.  Total  free  charge 
obviously  increases  because  of  increasing  propagation  but,  surprisingly,  peak  conduction 
currents  (due  to  the  primary  streamer)  actually  decrease  (Section  5).  However,  for  those 
prebreakdown  events  which  include  the  appearance  of  additional  ionisation  there  is  a 
clear  effect  of  distance.  For  the  smallest  gap  studied  (5.3cm)  then  the  second  corona 
event  exhibits  a  peak  current  of  the  same  order  as  the  primary  streamer  current  whilst, 
for  longer  gaps  (up  to  18.3cm)  the  second  corona  current  may  be  7-8  times  larger  than 
the  primary  current  for  the  same  E/p.  Table  5.1/2  contain  data  illustrating  this  point. 
Unfortunately  an  insufficient  number  of  events  were  recorded  to  provide  an  exhaustive 
study  -  because  of  the  relative  rarity  of  the  events  and  the  need  to  make  comparisons  for 
similar  values  of  E/p  at  different  propagation  lengths  -  but  the  trend  is  unmistakable.  If 
confirmed  then  this  could  well  lower  further  the  field  required  for  breakdown  in  the  pre¬ 
sence  of  tong  streamers  or,  alternatively,  increase  the  probability  of  a  breakdown  at  the 
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higher  reduced  fields. 

Phenomenology  apart,  the  mechanism  of  such  a  breakdown  is  not  clear.  The 
model,  developed  in  Section  3  and  extended  in  Section  6,  is  useful  because  is  calculates 
reasonably  well  the  growth  of  free  charge  with  distance  as  a  function  of  reduced  field 
and  so  provides  useful  predictive  information  since  any  breakdown  must  depend  on  the 
generation  of  such  free  charge.  It  is  difficult  to  assess  the  role  of  tha  cathode  defin¬ 
itively  although  two  effects  are  observable,  namely: 

(i)  electon  injection  at  streamer  arrival,  and 

(ii)  the  appearance  of  fast  ionising  waves,  the  generation  of  which  appear  to 
coincide  with  streamer  arrival  at  the  cathode.  There  is  some  evidence  that  the 
waves  -  if  indeed  they  are  waves  -  can  enhance  the  ionisation  near  the  origin  of 
the  streamer  system  (e.g.  Fig.4.7). 

The  magnitude  of  the  electron  injection  is  relatively  small  and  takes  place  at  such 
large  distances  from  the  origin  (i.e.  d/ve  timescale  of  the  breakdown)  so  that  it  is  dif¬ 
ficult  to  see  how  a  contribution  to  the  breakdown  could  be  made  in  the  present  geometry 
and,  if  a  valid  mechanisn  must  be  equally  relevant  to  aircraft  streamers,  then  such  a 
contribution  must  be  ruled  out  since  the  only  cathode-like  structure,  namely  the  preci¬ 
pitation,  would  provide  considerably  less  local  electron  injection.  That  streamer  impact 
with  such  precipitation  might  yield  ionising  waves  must,  however,  remain  a  plausible 
possibility.  An  assesanent  of  this  has  already  been  presented  elsewhere  by  Bicknell  and 
Shelton  (1986). 

The  growth  of  free  charge  represented  by  equation  6.2 

A(x)  =.AeTx.  2T=T,/vs 

suggests  an  alternative  approach  to  the  problems  of  providing  an  explanation  for  long 
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streamer  breakdown.  Therm alisation  of  the  streamer  channel  requires  an  electron  den 

23  -3  21 

sity  of  ^lO  m  compared  with  a  measured  value  for  small  gap  streamers  of  *'10 

m*^,  If  n(x)  is  the  electron  density  at  x  then  assuming  a  similar  constant  cross-section 

for  each  channel 


n(x)  =  n(0)N(0)erx/N(x) 


where  N(x)  is  the  number  of  channels  at  x.  Clearly  the  requirement  that  the  electron 

21  3 

density  remains  at  the  value  n(0)  (*>10  m'  )  implies  that  the  channel  number  growth 
rate  is  similar  to  the  charge  growth  rate  or 


N(x)  =  N(=  )e** 


There  is  little  experimental  data  available  regarding  the  behaviour  of  N(x)  particularly  at 
sub-atmospheric  pressure  but  replotting  the  data  of  Bicknell,  Sadik  &  Tang  (1980)  for  the 
growth  of  N(x)  at  E/p  =  7.92  (atmospheric  pressure)  suggests  a  value  for  JT of  20.5.  At  the 
same  E/p,  results  from  this  present  study  (Fig.6.2)  provide  2T(=P/v  )  of  35.6.  On  this 
evidence  then  the  electron  density  is  growing  with  propagation  distance  as 

n(x)  =  n(#)e  15‘lx 

2  23  21 

or  if  the  required  enhancement  is  10  (10  / 1 0  )  then  the  required  x  is  lnl00/15.1  ~> 

30cm.  For  lower  reduced  fields  this  distance  would  be  larger  because  of  the  smaller  dif¬ 
ference  between  the  two  growth  rates.  An  investigation  conducted  in  gaps  larger  than 
the  present  20cm  arrangement  with  facilities  to  monitor  2 f  both  for  A(x)  and  N(x)may 
well  be  fruitful  and  conclusive. 
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Appendix  A 


proqram  "h“ 

by  Rod  Shelton,  upd at  ed  :  06 . 04  .  97 

this  program  read*  data  f l \ *s  ertattd  by  “corona"  and  "sample", 
the  data  can  be  analysed  in  real  time,  using  cur«or-drutn  routines 
the  input  file  must  be  called  ‘B:D 


PROGRAM  DATAUSEIi 
CONST 

SKEW  -  1 05 i 
TYPE 

MATRIX-ARRAYC0. . 3,0. . 33  OF  REAL* 

VECTOR-ARRAYC0. . 3]  OF  REAL; 

RA-ARRAYC0. .2331  OF  REAL; 

VAR 

I NTC  , INTD1 , INTD2 , INTD7: REAL; 

A: MATRI X; 

C « ERR: VECTOR! 

NA: INTEGER; 

HEADING: ARRAY Z 1. .801 OF  CHAR; 

ICMAX  ,  TF'K  ,  IDFOSMAX  ,  I DNEGMAX  f  T X  ,  QGMAX  ,  ORATE  ,  QERROR ,  DOST AR  T  ,  DQEND :  PEAL ; 
DORATE , DGERROP , DQ I NJECTED , REDUCED  I ELD ; REAL : 

CALL23 ,CALL26,P0LYFITF0UND: BOOLEAN; 

STATUS, PSTATUS: INTEGER; 

YH , VS , F I TA , F I TB , ERARATE : REAL; 

YGi INTEGER; 

GAPLEN , STRLEN , K VOLTS , VE , VS , VV , F I ELD , PRE : REAL ; 

FILENO: INTEGER; 

F ILENUMBER , DATAFILE , DBASE: TEXT ; 

BUTTONPRESSED , F I RSTT I ME , UNDEF I NED , DMP , GRSW , F I TFOUND , CORRECTED s  BOOLEAN : 
ENDP AGE , PF I TFOUND : BOOLEAN ; 

XI ,X2,Yl , Y2 , ZEROLEVEL , ARATE , I NTEGR AL : RE AL ; 

XCUR , YCUR , XCHAN, YCHAN , BUTCHAN , I X 1 ,1X2, IY1 , IY2, IRESP,R8UTt INTEGER; 
IFREEiRA; 

PRESS , CLOCK V ALUE , HUM l D , X  T 1 MB , Y  AMPL , RES 1 ST  OR ; RE  AL ; 

MINUS1 , MINUS2: INTEGER; 

P2EROLE VEL , PARATE , P INTEGRAL : REAL ; 

PI  XI ,PIX2,PIY1 ,PIY2* INTEGER; 

ICONDiRA; 

PPRESS , PCLOCK VALUE , PHUM I D , PXT I MB , P YAMPL , PRES  I STOR : REAL ; 

PNRUN , PN  X  P  TS , PNPO I NTS , PF I RSTSAMPLE  s  l N  T  EGER ; 

PCHANNEL i CHAR ! 

NR UN , NXPTS , NPO I  NTS , F I RSTSAMPLE , OLD XCUR , OLDYCUR » INTEGER ; 

OLD IX 1 , OLD IX 2: INTEGER; 

ACHAR, SELECT  I ON, CHANNEL i CHAR; 


pr  oc  edur  e  wa  1 1 
no  global  usage 

waits  for  a  time  determied  by  the  passed  parameter 


PROCEDURE  WAIT (PERIOD: INTEGER) | 

VAR  COUNTER I  INTEGER; 

BEGIN 

FOR  COUNTER* *0  TO  PERIOD  DO  COUNTER; -COUNTER 
END;  (•  of  wait  *> 

(•**#••«**««#»*••••««#•»•««»•«**•«••«»»•#•«#*««««*•••«•«••#• 
procedure  svccode 
no  global  usage 

inserts  machine  code  patch  to  enaole  sve  data  transfers 


PROCEDURE  SVCCODE I 
8EGIN  ( eSVCCODE*  > 

(*  code  for  output  to  sve  begins  here  *> 

(*  the  first  n  characters  in  memory  starting  at  eF90l  are  read  out  *> 
(•  n  is  the  value  in  #F900  *> 

C*  see  assembler  Listing  for  more  details  •) 

(*  ROD  SHELTON...  DEC  83  •  > 

POKE 16 ( HF910 , S003A) t  POKE  1 6 < #F9 12, S4FF9) * 

POKE  1 6 ( #F9 1 4 , S00 1 1 )  |  POKE 1 6 ( 6F  916,61 3F  9 )  t 
POKE  1 6 ( #F91 8 , 447 1  A) ;  POKE  1 6 ( *F9 1  A , #B2DB ) } 

POKE  1 6  ( *>F 9 1  C  ,  » DA0F  >  j  POKE  1 6  ( »F9 1  E  ,  #F9 1  A  >  ; 
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PQKEl*i#F9.:0,#D“73>  :  PO  El6(#F922«  #0D&1  •  : 

PO  Eld  (#F924  .  #1  7C2)  s  PO  E  1  ©  i  #F  926  ,  #C9F  9  '  • 

•'  *  cod»  far  input  from  ivc  begin*  here  ♦) 

(*  n  character*  are  read  in  from  tne  svc.  starting  at  #F90i  *) 

i*  n  is  the  value  in  #F900  *> 

(*  ROD  SHELTON...  DEC ' 33  *> 

POKE  16  <  #F930  ,  #003A  )  j  POKE  1  6  (  #F932  ,  #4FF9  *  ; 

POKE  16  <#F97.4  ,  #01  11  )  ;  P0KE1  6  «#F936,  #DBF9V  • 

POK  E 1 6 ( #F938 , #0782) j  POKE  16 ( #F93A , #37DA>  : 

POKE  1 6  <  #F93C • #DE*F9 ) j  POKE  16  (  #F93E  ,  #1281 )  ; 

POKE  1 6 ( #F940 , #0D 1 3) s  P0KE16<#F942,#37CC)  : 

P0KE16 <#F944 .#C9P9) j 

(*  CODE  TO  DUMP  SVCSCREEN  FOLLOWS  *> 

POKE (63324  , CHR (62  >>; 

POKE (63823  ,CHR (32  )); 

POKE  (63826  t  CHR  ( 30  >M 
POKE (63327  , CHR (233  >>x 
POK  E ( 63823  . CHR ( 253  >  > s 
POKE (63829  , CHR (30  >  > i 
POKE (63830  , CHR ( 253  ) > ; 

POKE (63831  , CHR (235  >); 

POKE  <  63832  , CHR (14  ) >  j 
POK E (63833  , CHR (3  >  > : 

POKE ( 63834  , CHR ( 30  ) > ; 

POKE (63833  , CHR ( 1 3  >); 

POK  E (63836  , CHR (203  >  >  j 
POKE (63837  , CHR (5  )); 

POKE (63838  , CHR ( 0  >>j 

POKE (63339  , CHR (14  )); 

POKE (63840  , CHR ( 5  >>$ 

POKE ( 6384 1  , CHR (30  > ) } 

POKE (63842  , CHR (10  )); 

POKE (63843  , CHR ( 203  ))j 
POKE (63844  ,CHR<3  >)j 
POKE (63845  tCHR(0  )); 

POKE (63846  , CHR ( 1 4  ))j 
POKE (63347  ,CHR<3  >>* 

POKE ( 63848  , CHR (30  ) ) t 
POKE (63849  , CHR (27  >)| 

POKE (63830  , CHR (203  >); 

POKE <63831  ,  CHR (3  >)| 

POKE <63832  ,CHR(0  )), 

POKE  <  63833  , CHR (14  >  )  j 
POKE <63834  , CHR ( 3  ))j 
POKE (63833  , CHR <30  ) ) j 
POKE ( 63836  , CHR (76  ) ) j 
POKE (63837  f CHR <203  >  > ; 

POKE (63838  ,CHR<3  >>j 
POKE (63839  ,CHR<0  >)j 
POKE ( 63860  , CHR (14  ) ) ; 

POKE ( 6386 1  , CHR (3  ) ) j 
POKE <63862  , CHR (30  ) ) ; 

POKE (63863  « CHR (0  )), 

POKE ( 63864  , CHR ( 205  > ) j 
POKE (63863  ,CHR<3  )); 

POKE ( 63866  , CHR ( 0  ))j 
POKE (63867  t CHR <14  )); 

POKE (63868  , CHR  <  3  >>: 

POKE ( 63869  , CHR  <  30  )  >  x 
POKE (63870  , CHR ( 2  >)j 
POKE ( 6387 1  , CHR  <  203  >  >  ; 

POKE <63872  , CHR ( 3  )){ 

POKE (63873  ,CHR<0  >); 

POKE (63874  , CHR ( 6  >>| 

POKE (63875  , CHR (233  )>| 

POKE (63876  , CHR <197  ))| 

POKE (63877  , CHR (219  >)j 
POKE (63878  , CHR (178  ))* 

POKE <63879  , CHR ( 7  >>| 

POKE <63880  , CHR (218  ))| 

POKE (63881  , CHR ( 1 32  >>j 
POKE <63882  , CHR (249  >>j 


POLE  *0  733  3 

,  CHP  •  719  , 

)  >  i 

POLE  <  6  7884 

,CHP  177 

>  )  : 

PO E ( 63883 

,  CHP '95  > 

>  i 

POLE (63888 

, CHP ' 1 4  1 

>  ; 

POLE  <6  7387 

, CHP  *  5  ) > 

POLE (63888 

,CHF -213 

> > : 

POK E  <63889 

, CHP  <245 

>  > : 

POLE (63890 

,  CHP-  <  1 97 

>  > : 

POLE (63391 

,  CHP < 205 

>  >  ; 

POLE  <63392 

, CHP  <5  >  > 

i 

POLE (63393 

, CHP ‘0  >  > 

\ 

POLE <63394 

, CHR  « 1 93 

>  >  i 

POLE (63095 

, CHP (241 

>  > ; 

FO E  <6  3896 

, CHP ( 209 

>  > : 

POKE  < 63397 

.CHP -  205 

» ) : 

P0»  E  <63898 

,  CHR  5  >> 

ROVE  <  67899 

.CHP‘0  o 

: 

POLE  <63900 

, CHR < 193 

)  >  : 

POKE (63901 

, CHR  <5  ) ) 

j 

POLE (63902 

, CHF <194 

)  >  s 

POKE  <63903 

, CHR  < 1 32 

)  > : 

POKE  <63904 

. CHR ‘249 

>  > : 

POLE (63905 

, CHP  <219 

)  >  ; 

POLE ( 67906 

, CHR <173 

)  )  : 

POLE (63907 

, CHR  1 7  >  > 

; 

POKE  <63908 

, CHR *218 

)  >  ; 

POLE <63909 

, CHP ' 161 

)  )  ; 

POLE  <63910 

, CHR <249 

)  ) : 

POLE  <6391 1 

.CHP <219 

>  >  i 

POKE  <63912 

< CHR <177 

>  > ; 

POKE  <63913 

, CHR (95  )>; 

POKE  <63914 

.CHR <14  ) 

)  ; 

POKE  <63915 

. CHP  <  5  )) 

; 

POKE  <63916 

,CHR(213 

)  > : 

POKE  <63917 

, CHR (245 

>  >  • 

POKE ( 639 1 8 

, CHR < 197 

>  )  : 

POKE <63919 

, CHR  <  203 

>  > ; 

POKE  <63920 

,  CHR (5  >  > 

; 

POKE  <63921 

,  CHR ( 0  )) 

» 

POKE  <  63922 

, CHP ( 193 

) )  i 

POKE <63923 

, CHP (24 1 

)  >  i 

POKE (63924 

, CHR  <  209 

)  >  i 

POKE <63925 

, CHR (205 

>  )  • 

POKE (63926 

. CHP (3  ) > 

5 

POKE <63927 

. CHP  <0  >> 

} 

POKE (63920 

, CHR (30  >>| 

POKE (63929 

, CHR (233 

>  >  1 

POKE (63930 

, CHR (253 

)  >  J 

POKE (63931 

, CHR (61  ) 

)  1 

POKE (63932 

. CHP ( 194 

>  > i 

POKE (63933 

, CHR (85  > ) | 

POKE (63934 

, CHR (249 

)  >  l 

POKE (63935 

, CHR <201 

)  )  ; 

<♦  ROUTINE  FOR  DOUBLE  DENSI^v  •> 

END;  (*  OF  SVCCOOE  •> 

(*  N.B. 

th«  jpjcf  from  4F900  tc  •F^43  i«  thu*  used  by  these  routines 
the  is  from  IIF902  to  4F90F  and  can  hold  IS  character*  only 

for  the  first  character  1*  clearly  a  control  constant.  *) 

function  entier 

return*  the  greatest  integer  x 

FUNCTION  ENT  IER  (XsREAL)  :  INTEGER? 

BEGIN 

IF  X*0  THEN  ENTIERi -TRUNC < «  • -1  ELSE  ENT I£Ri -TRUNC < X ) x 
END| 

f unct i on  f rac 

return*  the  value  x-entitr  (x) 


FUNCTION  FRAC (Xi REAL ): REAL: 


) 


IV 


rr 


d£  j  I  N 

FSAC:  -x-ENTIER.  X '  ; 
END: 


*  uOCt  1  On  yt:;) 


FUNCTION  v t X: INTEGER) : INTEGER; 

BEGIN 

Y :  »YG-*-ENTIER<  IFREECX3/YH)  ; 

END: 

PROCEDURE  3HOWMATRI X (VAR  A: MATR I  X ; VAR  NA:  INTEGER); 

VAR 

ROW, COL: INTEGER; 

BEGIN 

FOR  ROM: =0  TO  NA  DO 
BEGIN 

WRITEt  t  >: 

FOR  COL : *0  TO  NA  DO  MR I TE i AC  ROW , COL  3 ,  >; 

WR I TELN  <  ' 3 ' ) s 
END; 

END; 

PROCEDURE  REDUCE ( VAR  A , B: MATRI x j vAR  NA . NB , I , J : INTEGER ) ; 
V  AR  , 

ROW, COL: INTEGER; 

BEGIN 

FOR  ROW: =  0  TO  NA  DO 
FOR  COL : =0  TO  NA  DO 
BEGIN 


IF 

'ROW 

*  I) 

AND 

(COL  ‘ J) 

then 

BC  ROW 

, COL  3 : =AC  ROW , COL  J 

IF 

(ROW 

•'  I  > 

AND 

(COL  J) 

then 

B  C  ROW 

t COL- 1 ] : =AC ROW , COL  3 

IF 

(ROW 

'  I  > 

AND 

(COL  vj) 

then 

BCROW-1 

,COL  ] ; = AC  ROW , COL  3 

IF 

(ROW 

>1) 

AND 

(COL  vj) 

then 

B  C  ROW- 1 

, COL -13: =A [ ROW , COL  3 

END; 

NR  j  =  NA- 1  ; 

ShCWMATR  I  x  <  B ,  NB )  ; 

END: 

FUNCTION  DETERMINANT (VAR  A: MATRIX; VAR  NA : I NTEGER ): REAL ; 

VAR 

DETjREAL; 

BlMATRI X  j 

NB, I ,F: INTEGER; 

COL: INTEGER; 

BEGIN 

IF  NA-0  THEN  DET:»AC0,0] 

ELSE 
BEGIN 
DET: -0; 

FOR  I : *0  TO  NA  DO 
BEGIN 
COL: -0; 

REDUCE  (A,B,NA,NB,COL,  I  )  ;  (*  ELIMIATES  COL  ( I )  *<  ROW(0>  .  .RESULTj*B  *> 

IF  ODD  (  I  )  THEN  F  i  ■—  1  Ei.SE  F:-l; 

DET;  »DETfF*AC0, I ] *DETERM I NANT ( B , NB )  ; 

END; 

END: 

DETERMINANT: -DET j 
ENDl 

FUNCTION  COFACTOR (VAR  At  MATRI X ; VAR  NA ,  I  , J ;  INTEGER >: REAL ; 

VAR 

BiMATRIX ; 

ROW , COL , NB , F : INTEGER; 

BEGIN 

REDUCE ( A , B , NA , NB , I , J ) ; 

IF  ODD  <  I  ♦  J  )  THEN  Ft  — l  ELSE  Fi-1j 
COFACTOR: -F *DETERM 1 NANT ( B , NB ) ; 

END; 

FUNCTION  POLY (VAR  X  *  I NTEGER >: REAL ; 

VAR 

A i REAL; 

N: INTEGER; 

BEGIN 
A;  -0; 


FOR  N:-a  ro  DO  4:  *A-CCNl*FQwER  \  X  ,N)  ; 

POLY: -Ai 
END; 

PROCEDURE  POLYFIT: 

VAR 

J,KfN&,ROW,COL, I: INTEGER; 

RESIDUAL . DETA ,HHH: REAL; 

B: MATRIX; 

•/:  VECTOR; 

BEGIN 

FOR  J:«0  TO  v  DO 
BEGIN 
YCJ]j»0: 

CCJ3: -0: 

FOR  K : -0  TO  3  DO  ACJ,k3j*0j 
END; 

FOR  TO  1x2  DO 

BEGIN 

FOR  K:»0  TO  3  DO 
BEGIN 

YCK3:*YCK3^IFREECJ]*POWER(J  ,»•  >  t 
AC0»KJ:*AC0,t-  1  ♦POWER (J,K>  ; 

END; 

FOR  ki-1  TO  3  DO  ACK ,33:-ACh , 33+POWER < J , <K*3> > ; 

END; 

FOR  ROW: *1  TO  3  DO 
FOR  COL: *0  TO  2  DO 
A  CROW  ,  C0L3  s  -ACROW— 1  ,COL*U; 

NA  x  *  3 ; 

WRI TELN (  MATRIX  COMPUTED > ; 

SHOWMATRIX ( A, NA)  ; 

DETA: -DETERMINANT (A, NA) ; 

WRITELN(  DETERMINANT  FOUND'); 

FOR  I : *0  TO  3  DO 
FOR  J : -0  TO  3  DO 

CC  U  t »C C  U ♦COFACTOR <  A , NA , J , I >  * Y  C  J  D  ; 

WR I TELN <' COFACTOR  MATRIX  S.  INVERSION  COMPLETE ')  ; 

WR I TELN ( ' DET  A : ■  '  , DETA) ; 

FOR  I:-0  TO  3  DO  CC 13 i -CC II/DETA; 

RESIDUAL: -0; 

FOR  Ii-IXl  TO  1X2  DO 

RESIDUAL: -RES I DUAL^SQR ( I FREE C 1 1 -POLY ( I ) ) ; 

HHH* -RESIDUAL/ < DETA* ( I X2- I  X 1 ♦ l -4 > ) j 
WR I TELN ( 'RESIDUALS  COMPUTED  > ; 

FOR  I:-0  TO  3  DO 
BEGIN 

REDUCE ( A , 8 , NA , N8 , 1 ,  I)  ; 

ERRC I  3 i -SORT (HHH*DETERM INANT (B ,NB) > ; 

END; 

WR I TELN <  ABSOLUTE  ERRORS  IN  COEFFICIENTS  CALCULATED '> ; 
POLYP ITFOUND i -TRUE ; 

WRITE < ' > ' )  |  READLN; 

END; 


Auction  yfit; 

finds  th*  (real)  value  of  the  exponential  fit  done  by  expfit 


FUNCT ION  YF I T ( X : INTEGER ) i I NTEGER ; 

VAR  J : REAL ; 

BEGIN 

IF  F ITFOUND  THEN 
BEGIN 

Ji-EXP(FITB) ; 

IF  (IFREECIXl J<0>  THEN 
J : -J*EXP (F 1 TA*X ) /YH; 

<*  YFIT; -YG+ENTIER < J > j *> 

END 

ELSE  IF  POLYFITFOUND  THEN  J : -POL V ( X ) / YH; 
YFIT:-YG*ENTIER(J)  ; 

END; 


procedure  normal  screen 


vi 


no  global  usaqe 

r«stor»s  80  column  format  to  svc  display 


PROCEDURE  NORMALSCREEN; 

VAR  C: INTEGER; 

BEGIN 

Ci -CRM (6, 27 > ; 

C; »CPM16,49>  j 

END;  <*  of  nor*Alscr«#o  *) 


procedure  clearscreen 
no  global  usage 

clears  the  screen  exiting  with  svc  in  graphics  -format 


PROCEDURE  CLEARSCREEN; 

VAR  Cs INTEGER; 

BEGIN 

NORMALSCREEN; 

Ci*CPM(b,26>  ;  (*  CURSOR  HOME  $<  CLEAR  SCREEN  *> 

C: *CPM<6,27> ; 

C S *CPIi <  6 , S2 )  ;  <*  BACK  TO  GRAPHICS  MODE  *> 

END; 

i  *♦*♦****■**************-*♦**************♦*****♦*♦*** 

procedure  point 

global  usage:  .none 
requires  graphics  format 
draws  a  point  at  XX, YY 


PROCEDURE  PO I NT (XX, YY : CHAR ) ; 

VAR 

SKEWX; CHAR; 

BEGIN  (»  point  *> 

IF  ( YY <  *CHR ( 235 ) )  AND  ( YY>*CNR (0) >  THEN 
BEGIN 

SKEWX : »CHR  <  ORD (XX) +SKEW ) j 
POKE(*F900,CHR(A) >  ; 

POKE <#F901 , CHR (27)  )  j 
POKE ( 4F902 , CHR (83) >  ; 

POKE ( #F903 , SKEWX )  | 

POKE ( #F904 , CHR ( 00 ) > ; 

POKE ( #F905 , YY ) | 

POKE ( *F906 , CHR ( 00 ) ) ; 

USER (4F910) j 
END; 

END;  <*  of  point  *> 

(»««••«••»*»«**#«*••»*•«*««««•**•**««»•• 
procedure  unpoi nt 

global  usage i  none 
requires  graphics  format 
rubs  out  anything  at  XX, YY 


PROCEDURE  UNPOINT < XX, YYiCHAR)  ; 

VAR 

SKEWX i CHAR; 

BEGIN 

IF  <YY>-CHR<0> >  AND  < YY< -CHR (235) )  THEN 
BEGIN 

SKEWX  t *CHR ( ORD (XX) +S* EW ) ; 

POKE ( #F900 , CHR ( 6 ) ) ; 

POKE (4F901 ,CHR  <27) ) ; 

POKE < *F902 , CHR ( 02  >  > } 

POKE ( #F903 , SKEWX ) ; 

POKE (#F904 , CHR (00) ) j 
POKE  <  #F903 , YY ) ; 

POKE ( #F906 , CHR ( 00 ) ) ; 

USER  <#F910) | 

END; 

END; 

(a************************************************************************* 
procedure  cur  sorrow  (:<) 

♦••♦a************************************************************#********) 


PROCEDURE  CURSORROW ( X s INTEGER) i 


BEGIN 

P0> E  * #F900 , CHR ( 4 )  > ; 
POfc E < #F90l , CHR (27) )  i 
POk  E ( #F902 , CHR  < #3D > >  l 
POk  E ( #F903 , CHR ( X )  )  i 
POk  E ( #F904 , CHR ( 5)  >  l 
USER (#F910) ; 

END; 


proctdur*  move 


PROCEDURE  HOVE < XX ,YY« CHAR) ; 
VAR 

SKEWX ; CHAR ; 

BEGIN 

3V  EWX  J  =*CHR  ( ORD  (XX)  ♦St  EW )  ; 
POKE t #F900  ,CHR (6) ) ; 

POKE (#F901 ,CHR  <27) > ; 

POKE ( #F902 , CHR (109)); 

POKE ( #F 903 , SKEWX  > ; 

POKE  <  #F904 , CHR ( 0  >  > } 

POKE ( #F905 , YY  > ; 

^Ok  E  (  #F  906  y  CHR  (  0  )  )  ; 

USER  (#F910> j 
END; 

'•*«•*»*«»*»***•***»»*«**«• 
procedure  complement 


PROCEDURE  COMPLEMENT (  XX  , YY i CHAR ) } 

VAR 

RESULT , SKEWX ; CHAR; 

BEGIN 

SKEWX ; *CHR (ORD  ( XX ) ♦SKEW) ; 

POKE  <  #F900 , CHR ( 6 ) ) ; 

POKE ( #F90 1 , CHR ( 27 ) ) ; 

POKE ( #F902 , CHR  <  #54  > ) ; 

POKE ( #F903 , SKEWX ) ; 

POKE  <  #F904 , CHR <  0  >  )  • 

POKE  (  #F905 , Y Y ) ; 

POKE <#F 906, CHR (0) ) , 

USER  <#F9t0> ; 

POKE ( #F900 , CHR  < 1) ) ; 

USER (#F930) J 

RESULT; -PEEK (#F901 f CHAR) | 

CASE  RESULT  OF  CHR (0) *  PO INT ( XX , YY ) ; 

CHR  < 1 ) i  UNPOINT < XX ,YY) 
END; 

END; 


procedure  line 

PROCEDURE  LINE <XX, YYiCHAR) ; 
VAR 

SKEWX: CHAR; 

BEGIN 

SKEWX; -CHR (ORD (XX) +SKEW) ; 
POKE  <  #F900 , CHR ( 7 ) ) ; 

POKE ( #F901 , CHR (27) ) | 

POKE ( #F902 , CHR ( 1 08) ) : 

POKE ( #F903 , SKEWX ) ; 

POKE ( #F904 , CHR  <  0 )  )  ; 

POKE (#F905, YY) ; 

POKE ( #F906 , CHR (0) ) ; 

POKE < »F9®7 , CHR ( 1 )> ; 

USER (#F910> | 

END  | 


procedure  dmprq 

global  usage;  none 
requests  sve  to  dump  screen 


PROCEDURE  DMPRQ j 


) 


vii1 


VAR 

C: INTEGER; 

BEGIN 

Ci -CPM(6,#1B> ; 

C:-CPM<6,#64) ; 

C*-CPM<6, #50) ; 

END? 

(••*•••••*«••*•••»•***«** 
proctdurt  graphic* 
global  usagct  none 
enters  graphics  mode 


PROCEDURE  GRAPHICS; 

VAR  C: INTEGER; 

BEGIN  (*  graphics  *> 

C?  -CPM (6,27) ; 

C: -CPU (6,52) ; 

END;  (*  of  graphics  *) 


procedure  initialise 

global  usage: UNDEFINED, FIRSTTIME, XI ,X2,V1 ,Y2 

0LDX1 ,0LDX2 , RBUT , XCHAN , YCHAN , 6UTCHAN , XCUR,  YCUR  ,  I  X  1  , 1X2, IY1 , IY2 
OLDXCUR , OLD YCUR , PRESS , F I RST SAMPLE , CLOCK  VALUE , HUM  I D , X  T I  MB , YAMPL 
RES  I ST OR , RUN ,NXPTS , NPQ I  NTS 
sets  initial  values 


PROCEDURE  INITIALISE; 

VAR  A: INTEGER: 

BEGIN 
ICMAX : =0; 

I DNEGMAX : =0; 

TPK: =0; 

TX : =0; 

IDPOSMAX: -0; 

QGMAXi«0; 

□ERROR* -0; 

DOST ART ;=0; 

DOENDi *0; 

DQRATE: -O; 

DOERROR; -0; 

PFITFOUND: -FALSE;  FITFOUNDi -FALSE; 

WRITELN<  THIS  VERSION  READS  DATAFILES  FROM  DRIVE  “A- .  ); 

WRITELN; 

WRITELN(  V20.05.07,  corrected  for  cable  impedance 
WRITELN; 

WAIT ( 1000) ; 

VSi-0; 

I X 1 1 -30; 

1X21-90; 

ENDPAGE t -FALSE ; 

GRSW; -FALSE; 

DMP: -FALSE; 

MINUS 1 i-l ; 

MINUS2: -~1 ; 

UNDEFINED: -TRUE? 

FIRSTTIME* -TRUE; 

XI i-0;  Y 1 : — 0;  X2: -0:  Y2:-0; 

OLD I  X 1  *  —  I  X 1 ;  OLD I X2: -1X2; 

RBUT  *  - 100; 

XCHAN* -11: 

YCHAN* -12; 

BUTCHAN* -8; 

XCUR I  —  I  X 1 ; 

YCUR: -200; 

IYli -100; 

IY21-150; 

OLDXCUR* -XCUR; 

OLD YCUR* -YCUR; 

PRESS* -0; 

F 1RSTSAMPLE* -0; 

CLOCK VALUE *  —0 ; 

HUMID* -0; 


ix 


XTIMB?  -0? 

YAMFL? -0? 

RESISTOR? «0? 

NRUNi -0? 

NXPTS: -0; 

NPO I NTS*  -0? 

END?  (*  of  initial***  *> 


proc*dur*  beginning 

global  usage: DATAFILE 

read*  th*  heading  in  the  input  file 


PROCEDURE  BEGINNING? 

VAR  A* INTEGER? 

BEGIN 

RESET (DATAFILE,  As D«  .  '>? 

RESET (FILENUMBER,  A: LOOKING  .AT  '? 

REWRITE (DBASE, ’ A: DATABASE . NEW ' > ; 

A:«l: 

REPEAT  R£AD ( DATAF ILE , ACHAR) ? 

HEAD  I NG  C  A ] : * ACHAR  ? 

WRITE (ACHAR) ? 

A:-A+l? 

UNTIL  ACHAR-  '**? 

WRITELN; 

WRITELN? 

WRITECENTER  filenumer  »»>');  READ  (F ILENUMBER ,F ILENO)  ;  WRITELN  (FILENO)  ? 
WRITELN? 

WRITELN? 

WRITELN (' ENTER  D  (cm)  needle  height  (mm)  V  <KV)  and  p  (mmHg)  ')? 

WRITE < ' >>> 1 ) ? 

READ ( GAPLEN  ,  STRLEN  ,  KVOLTS , PRE ) ; 

WRITE (' E/p-  ’ , <kVOLTS/GAPLEN/PRE*1000) , '  V/cm.to rr  ?  enter  approx,  value  ' ) i 
READ (REDUCEDFIELD) ;  WRITELN? 

STRLEN  s -GAPLEN- ( STRLEN/ 1 0  >  ? 

FIELD* -KVOLTS/GAPLEN*1.0E+5?  <*  VOLTS  PER  METRE  ♦> 

VE*-1.0E+6*POWER< (FIELD* 1 . 0E-3/PRE* l . 0E*3> , 0. 713) *  1 . 0E-2; 

(*  from  gallimberti  ref.  data*  V*  (cm/s) -10~6* (e/p)  '0. 713 

where  */p  is  in  v/cm.to rr  ...  I  have  converted  to  V*  (m/*)  *> 

END?  <*  of  beginning  *> 


procedure  zero 

global  usage: RAW, FIRSTSAMPLE , 

finds  th*  mean  of  th*  ten  RAW  datapomt*,  the  zerolevel . 

if  any  on*  differs  by  more  than  ♦  /-  6  from  th*  first,  then  th* 
GARBAGE  flag  is  set 


PROCEDURE  ZERO  (VAR  A:  READ? 

VAR 

B* INTEGER? 

BEGIN 
At  -0? 

FOR  Bi -FIRSTSAMPLE  TO  FIRSTSAMPLE  -*-4  DO  A* -A*  I  FREE  I B3  ? 

A* -A/3 

END?  (*  of  zero  *> 


procedure  newevent 

gl obal  usage*  CHANNEL , NXPTS , NPO I NTS , PRESS , CLOCKVALUE , F I RSTSAMFLE 
HUMI D, XT  I MB, YAMPL, RESISTOR, RAW 
reads  th*  parametric  datafield  for  the  next  record 
reads  th*  integer  data  into  th*  array  RAW 


PROCEDURE  NEWEVENT? 

VAR  A, Bt INTEGER? 

ACHAR  *  CHAR  ? 

BEGIN 

REPEAT  RE AD (DATAF ILE, ACHAR)  UNTIL  ACHAR- ' \ ' ? 

READ (DATAFILE, CHANNEL , NXPTS , NPO I NTS , PRESS , CLOCKVALUE , F I RSTSAMPLE ) ? 
READ ( DATAF I LE , HUMID , XT  I MB , YAMPL , RES  I STOR  > | 

IF  CHANNEL-  A'  THEN  YAMPL * -YAMPL *RES I STOR/ 33. 3?  (*  correction  •  *) 

(*  Z0-70  OHMS,  R— 64 . 3  OHMS  *  FINAL  CORRECT  I ON > 

IF  CHANNEL-  B  THEN  YAMPL * -YAMPL *RES I STOR/ 33, 3| 


X 


l*  20-74  OHMS,  P-60.6  OHMS  *< 

FOR  A:  -F  IFSTShMPLE  TO  NP0INT3  DC 
BEGIN 

READ ( DATAF I LE , B  >  s 
I FREE  C A  3 : -B*v AMPL/ 1  28: 

END; 

ZERO  <  2EROLEVEL ) J 

FOR  Aj  -FIRSTSAMPLE  TO  NPOINTS  >0  I  —  EE C A3  i  * IFREEC A 3 -ZEROLEVEL ; 
END  j 


proctdur*  display 

g  1  oDal  usag*  i  NRUN  ,  NPO I  NTS ,  N  t  -  ~  .  -  3  RSTSAMPLE ,  CLOCKVALUE  ,  RES  I  STOP 
PRESS,  HUMID,  XT  1-3  .  »^MPU 

tabulate*  the  parametric  data-ield  for  this  record 

♦»»#»»«»»«♦•»#♦»»»»»«♦«»»»»»»»#♦•»<  ♦  *  »♦»♦*«♦•♦♦♦»**♦*»♦♦♦♦»♦»♦»**♦♦♦♦**»»♦«») 

PROCEDURE  DISPLAYS 
VAR  A i INTEGER; 

BEGIN 
A  ;  =  1  ; 

REPEAT  WRITE  (HEADINGCA3  )  ;  As-A-M  HEADINGC  A3*  '  S  '  ; 

WRITELNj 

WRITE  <  * '  , F I LENO: 4  ,  '  .  '  )  ; 

WRITE<NXPTSs2, '  >; 

CASE  STATUS  OF 

Is  WRITEt'Xc  *  )  ;  <*,YH*2.5E^.  ss  3,  mA/cm  •  ,  XT  I  MB*  1 . 0E+6:  8  s  3 ,  ‘  ui/cn  )  j  ♦  > 

2  s  WR I TE  <  '  Id  ');(*,  YH*2 .  5E-4  A  ;  2:1,  mA/cm  ’  ,  <T  I  MB*  1 . 0E+6 :8:3, 'us/cm' );*) 

3;  WRITE !'Ie  >;(*, YH*2.  ZE-*z  E:  3,  mA/cm  '  , XTIMBM . 0E+6: 0s 3 ,  ‘ ut/cm ' ) ; a) 
4;  WRITE  (  ‘  Qg  '  )  ;  <  *  ,  YH*2 . 5E+ 12:  5s  3 ,  '  nC/cm  '  ,  XT IMB*1 . 0E+6: 8 a  3 ,  ‘ us/cm ' ) j *> 

5s  WR I TE  ! ' dQ  '  > j  !*, YH*2. 5E>i2i S: 3,  nC/cm  '  ,  XT  I MB* l , 0E+6: 9:3,  '  us/cm ' ) ;•) 

6:  BEGIN 

WRITE (  Ax  ‘  >  ;  <  *  ,  YH*2. 5E-*  1 0:  8 :  3  •  ’nC/c/n  *  ,  XTIMB*VE*VS/  (VE+VS)  *10Cs  8s  3)  ;  *> 
WRITE  (  'cm*  '  ) 

END 

END*  <♦  OF  CASE  *> 

WRJTELNi 

WRITE  (GAPLENs  4j  1  ,  cm  '  ,  PRE:  3:0,  '  n.mHg  *  )  ; 

WRITE<FIELD*1.0E-5:5*2,  *  KWcm  '  ,  t- lELDM . 0E-5/PRE*l .  0E*3>  i  5:  2)  j 
WRITELN (  V/cm.torr ’)} 

WRITE! 'V*-  '  , (VS/1000) sSs 1  ,  '  Km/s' >; 

WRITELNj 

WRITE! ICMAX, '  j  ' , TPK ,  j  '.IDPOSMAx,  j)j 
WR I TELN ! I DNEGMAX ) j 

WRITELN!TXf'  ;  ,VS,'  j  '  , QGMAX ,  *  j  , QERRQR) j 

WRITE !DQSTART ,  j  ' , DQEND ,  '  j  '  ,  DGRATE )  ; 

WRITELN! ‘  j  , DQERROR) } 

WRITE tDQINJECTED,  '  |  ' ,CALL23,  '  )  >  ; 

WR I TELN  <  CALL26 ) j 
ENDj  ! *  of  display  •) 


procedure  anal  menu 

global  usages  XI  ,  X2,  YI  ,  Y2,  INTEGRAC,NU,UNDEFINED,NXPTS,CHANNEL 
displays  the  cursor  data  and  the  calculated  quantities 
then  displays  a  user  prompt  for-  response  choice 

PROCEDURE  ANALMENUs 
VAR 

SFsREAL; 

BEGIN 

IF  NOT  DMP  THEN  NORMALSCREENj 
DISPLAY; 

CASE  STATUS  OF 

1,2,3s  SF  t -1 . 0E+-3| 

4,5s  SFs-1.0E+9 
END \  ! *  of  case  *> 

(*  WRITE!'uS')| 

CASE  STATUS  OF 

1 ,2,3i  WRITE! ' ,mA  ) j  _ 

4,5s  WRITE ! ' ,nC ' ) 

END;  *) !•  of  case  •) (• 

WRITE!  ' s  ' , X 1*1 .0E+6I 8s 3,  ' ,  ' ,Yl*SFt8i 3,  ’  j  , X2* 1 . 0E*6a 8s 3) j 
WRITELN!  , ‘ , Y2*SFs  8s  Z) | *> 

WRITE  !  *  dT»  '  , (X2-X1 >*1 . 0E+6 a  5s 2 ,  uS  Ve-  ,  !  VE/ 1 000) s 4 s 2 ,  Km/s' >| 


V 

t 


X 


WR I TEuNj 

IF  STATUS  4  then 
IF  NOT  UNDEFINED  THEN 

WRITE*  sint.  ,  INTEGRAL*  1 . 0E*9*  12:  4  ,  nC  >s 
IF  FITFOUND  THEN 

WRITE!  satt.  fr.  (k l )  ,ARATE*i.0E-o*.  12:4,  KHz  ) i 
WRITELN: 

IF  NOT  DMP  THEN 
BEGIN 
WRITELNs 


WR I  TELN  i  '  enter  sel  tction  to. ...  > 

WRI TELN j 

WRITE!  'i-:  add  (X1,Y1> 

IF  CORRECTED  THEN  WR I TELN !  '  ■  I  d  ♦  I  c . 
WRITELN!  '2;«  add  <  X2 , Y2> 

WRI  TELN  <  ’  •-  3  find  attachment  rate 
WRITELN!  '*.4  :  find  integral 

WRITELN!  -.*  :  go  to  next  record 

WRITELN!  6-s  quit  program 

WRITELN!  7  :  snow  data 

WRITELN!  -3  *  swap  cursor  1$<2 

WR I TELN !  ' •  9  ;  s  swap  old  4  new 
WRITELN!  ’ -.21  Ic  pi-  ‘  22--  -rid  pi  23 
WRITELN!'  25*  dQ  growth,  start  end 
END; 


11'$  f i nd  I e  > ; 

ELSE  WRITELN!  *ld^Ic  ); 

•■I  2i  find  Qg  >  i 
'lZ*i  screendump  )  ; 

14.;  gr at i cul e  > ; 

16  :  ti no  dC  ' ; 

find  st r earner f unc t l on  A*:) 
10  :  find  streamer  velocity 
19":  cubic  polynomial  fit) 

•  20 : :  le*  I d* Ic  > ; 
pi.  24  Og  pf  V  decay  )  ; 

6  -dQ  pfc  27  output  to  disc.  >; 


!d-n ) /d  ) 
at  xi 

1  7 

dataf i el d 

T m  Sr  -Id 


END;  ( *  of  analmenu  *> 


procedure  adc 

global  usagesnone 

returns  the  voltage  for  a  particular  channel  of  the  adc  card 


PROCEDURE  ADC ( VAR  ANS , CHANNEL : I NTEGER > j 
VAR 

A , B: INTEGER; 

INPUT: CHAR; 

RESULTiREAL; 

BEGIN 

out ( 3 2 ♦CHANNEL , CHR ( 0 )  )  ; 

FOR  At-0  TO  3  DO  A: -A; 

INPUT: -INP  < 32) ; 

ANS » -ORD! INPUT) ; 

END;  (*  of  adc  ♦) 


procedure  cursor 

global  usages  QLDXCUR , OLD YCUR , XCUR , YCUR , RBUT 

hnd»  new  cursor  posi  1 1  on ,  stores  old  one  so  it  can  be  erased, 
and  assigns  an  appropriate  value  to  BUTTONPRESSED 


PROCEDURE  CURSOR; 

BEGIN 

OLDxCURi -XCUR; 

OLDyCURi -YCUR; 

ADC (XCUR, XCHAN) j 
ADC (YCUR, YCHAN) ; 

ADC ! RBUT , BUTCHAN  > ; 

IF  RBUT; 125  THEN  BUTTONPRESSED * -TRUE  ELSE  BUTTONPRESSED i -FAL SE ; 
IF  XCUR«.  FIRST  SAMPLE  THEN  XCUR*  -F  IRSTSAriPLE; 

IF  XCUR ; NPO INTS  THEN  XCURs -NPOINTSs 
END:  !*  of  cursor  *) 


procedure  dr awgr a t i cu 1 e 
global  usage: none 

draws  scope  graticule  on  svc  screen 
requires  graphics  mode 


PROCEDURE  DRAWGRA  T I CULE ; 
VAR 

A, 6, C» INTEGER* 

BEGIN 

<•  SQUARE  GRID  •) 

(*  FOR  Bi-0  TO  230  DO 
FOR  A*-«  TO  10  DO 


T" 


begin 

POINT  (CHR  A*25>  ,  CHR  (B)  m 
POINT  (CHR  (B>  ,CHR  (A*25)  >  : 

END ;  *> 

\*  AXES  ONLY  *> 

FOR  Aj-0  TO  250  DO 
PEG  IN 

POINT (CHR (0) ,CHR ( A) ) j 
POINT <CHR( l 30> , CHft ( A ) ) ; 

ENDj 

FOR  A:»0  TO  l50  DO 
BEGIN 

POINT (CHR (A)  , CHR (0)  >  ; 

POINT <CHR (A) , CHR (250) > ; 

END; 

FOR  Bx-1  TO  9  DO  FOR  C:-l  TO  3  DO 
BEGIN 

POINT (CHR (C) ,CHR (25*B> » ; 

POINT  t CHR ( 150-C) , CHR (25*6) >  : 

END: 

FOR  B:»l  TO  5  DO  FOR  C:-l  TO  3  DO 
BEGIN 

POINT (CHR (25*6> ,CHR (C) ) : 

POINT .CHR (25*6) , CHR ( 250-0 ) ; 

END; 

C:  «CPn  (6,#i  6)  ;  C :  *CPt1  ( 6  ,  *0C )  ; 

END; 


procidur#  4utorang»(sc«l0:r*«l  .ranges  integer > 
global  usage: I FREE 

Find*  the  parameter*  needed  by  drawtrace  to  plot  the  data  in  IFREECal 


PROCEDURE  AUTQRANGE (VAR  C; REAL ; VAR  Ax  INTEGER) ; 

VAR 

8 : INTEGER; 

G, I NT .MAX .MINI  REAL; 

BEGIN 

MAXx-0;  MINx-0; 

FOR  Bi -FIRSTSAMPLE  TO  NPOINTS  DO 
BEGIN 

IF  IFREE C  B  3  >MAX  THEN  MAX : - I FREE t B J ; 

IF  IFREEC83<MIN  THEN  M I Ni - 1  FREE C 8 3 ; 

END* 

INTi -MAX-MINI 
A i -TRUNC <-l*MlN/ INT*10) ; 

IF  A;0  THEN 

G i -LN<-1 *MIN/A) /LN( 10. 0) 

ELSE 
BEGIN 
A:  •  t  ; 

G:-LN(MAX/9> /LN (10.0) ; 

END; 

Ci - (  (POWER ( 10,FRAC (G)  ) *10) *1 ) / 10aPOWER  <10,ENTIER(G> )  ; 
Ci"C/25|  A i •  A*25 | 

END; 


procedure  curdriwixi j 


PROCEDURE  CURDRAW ( X i INTEGER) ; 

VAR 

A i INTEGER; 

BEGIN 

(•FOR  Ax  »-2  TO  2  DO 
BEGIN 

COMPLEMENT  (CHR  <  X*A)  , CHR ( Y ( X )  *A)  )  ; 
COMPLEMEN i (CHR ( X*A) , CHR < Y (X) -A) ) | 
END; 

COMPLEMENT  « CHR ( X  >  . CHR ( V  <  X )  )  > ; * ) 

FOR  A* --2  TO  l  DO 
BEGIN 

COMPLEMENT (CHR ( X*A> , CHR ( Y ( X ) -2) ) ; 
COMPLEMENT (CHR (*♦ 2)  ,CHR  < V ( X ) -A)  )  ; 


L 


UmPPl£n6n T  CHF  *  A  —A )  ,  CAR  <  i  ■.  K  i  * 2  >  i  i 
COMPLEMENT  :CHR  X -2>  ,  CHR  •.  v  '  X  *  -A>  »  ; 
EN  Ls 
END: 


procedurt  drawtrac* 

gl obal  u*ag»: FIRST SAMPLE .NPOINTS, RAW 

r»qtiir»»  graphic*  mod* 

draw*  crudt  data  itraiqht  trom  RAW 


PROCEDURE  ORAWTRACE; 

VAR 

A  ,  B  ,  C  ,  D : INTEGER; 

BEGIN 

CLEARSCREEN; 

AUTORANGE (VH.YG)  ; 

IP  GPSW  THEN 
BEGIN 

DR AWGRAT I CULE: 

FOR  A  j  -F  I RST SAMPLE  TO  NPOINTS  DO 
IP  (A  S)  OR  (  Y  l  A )  •  S )  OR  ( Y ( A )  245)  THEN 
BEGIN 

D; «vG*ENTIER ( IFREEC A] /YH) ; 

FOR  B: --3  TO  3  DO 
FOR  C : --3  TO  Z  DO 
UNPQINT iCHR ( A*B)  ,CHR  <  0*0  >  ; 

END: 

END; 

MOVE (CHR (F IRSTSAMFLE)  , CHR  < Y  < F I RSTSAMPLE >  )  >  ; 

FOR  Ai *F I RSTSAMPLE  TO  NPOINTS  DO 
L I NE  < CHR ( A  >  , CHR ( Y ( A )  )  ) ; 

IF  FITFOUND  OR  POLYFITFOUD  THEN 
BEGIN 

A: »FIRST5AMPLE; 

REPEAT 
A:  -A-M  ; 

UNTIL  (YFITtA)  ®>  AND  t  YF  IT  ( A)  *.  2SB)  ; 

A;  -AM  i 

WHILE  <YFIT(A)i0)  AND  <  YF  I  T  <  A )  •.  2S0  >  AND  (A-  -NPOINTS)  DO 
BEGIN 

HOVE ( CHR ( A- 1 ) , CHR ( YF I T ( A- 1 > ) > ; 

L I NE  <  CHR (A)  , CHR ( YF l T ( A )  )  >  ; 

IF  ( A  * • I  X 1 )  AND  (Av*lX2)  THEN  L I NE ( CHR ( A )  , CHR ( Y ( A >  >  >  ; 
A; -AM | 

END; 

END; 

Ci -CPM<6,#lB>  ;  C: -CPM <6,#0C> ; 

FOR  C  j ■ 1  TO  4  DO  WRITELN; 

WRI TE  <  '  * i g  d  ,  F ILENOi 3,  ' .  *  ,NXPTSj  1 ) ; 

CASE  STATUS  OF 
1 :  BEGIN 
WRI TELN ( • a ' ) | 

WRITELN; 

WR I TELN ( ' conducti on ' ) ; 

WR I TELN ( ' currtnt ,  Ic  >; 

END; 

2:  BEGIN 

WRITELN (  b  >  ; 

WRITELN; 

WR  I  TELN<  displac  fflitnt  '  )  ; 

WRI TELN  <  currtnt ,  Id  '  )  ; 

END; 

3s  BEGIN 
WRI TELN <  c  > 1 
WRITELN; 

WR I TELN  <  ' * i mp 1 •  * > ; 

WR I TELN ( ' currant ,  I*  ) ; 

END; 

4  i 

BEGIN 

WRI TELN (  d  >  j 
WRITELN; 


XIV 


wRITELN*  total  *r*«  ■; 

WRITELN'  chai  ge,  Oq  >i 
£N0; 

S:  BEGIN 
WR  [TELN  v  e  •; 
wPI TELN: 

WRITELN <  gtn*ritioo  o+  >; 

WRITELN<  charge  dO/sec  )  ; 

END; 

BEGIN 

WR  I  TELN  t  *  )  } 

WRITELN: 
wR  I TELN t  A>: 

WRITELN: 

END 

END  l 

WRITELN; 

WRITELNJ  a 
WRITELN*  p 
wR  1  TELN  *  E 
WRITELN* 

WfcITELNt  E/p 
WR I TELN  < 

WRITELN*  vs 
wR I TELN  <  Ve 
IF  FITFOUND 
BEGIN 

IF  iARATE-  0)  THEN  WR I TELN (  '  dec  ay  <rf q.  *  *  ELSE 
WRITELN*  *ARATE*1 . 0E-6)  s  ^2)  ; 

MR  I  TELN  i  ♦  , <EF ARAT£*1 . 0E-6) : 6: 2 ,  MHz  >  ; 

END: 


,3APLEN: 4:1,  cms  >  ; 

,PRE:  7:0,  mrm-fg  )  ; 

,  -FIELD*! . 0£-5> : 4:  1 )  ; 

KV/cm  >j 

,  (FIELD* 1 . 0E-S/PRE*1 . 0E*3> : 4; 2> : 
‘J  /cm.  torr  '  >  ; 

, <VS*1 .06-3) «5i  1 ,  '  ►  <*/*>; 

,  i VE*  1 . 06-3*  :5i  2,  '*  >  ; 

THEN 


wR I TELN  t 


growth 


♦req.  * ' > 


IF  NOT  UNDEFINED  THEN 
BEGIN 
WRITELN i 

WRITELN!  integral  ■  > : 

WR I TELN ( ( INTEGRAL*! . 0E*9) :*:2,  nC  ) ; 

END: 

IF  UNDEFINED  and  vFITFQUND*FAL3£>  Then  BEGIN  WRITELN:  WRITELN;  WRITELN  END: 
WR I TELN; 

WR  I  TELN*.  ‘  v*r  /hor  sc»l»'>| 

WR  I  T  E  <  y  I  )  ; 

CASE  STATUS  OF 

4:  WRITELNJ <YH*2.3E*10) :6*2, *  nC  > t 
1 f 2 ,3*  WR I TELN (  I YH*2 . 5E*4 ) i 6 1  2 «  *A  >  ; 

3,©i  wR I TELN (  ( YH*2. 3E* 1 > : 6: 3 «  C/e  > 

END: 

WR I TELN (  h:  ' ,  ( X  T I  MB* 1 . 0E*6  > i 7 1 2 ,  ‘  uS ' > J 
IF  UNDEFINED  THEN 
BEGIN 

CURDRAW (1X1)1 
CURDRAW ( I X2>  j 
END 
ELSE 
begin 

FOR  A: ■ I  X  1  TO  1X2  DO 
BEGIN 

MOVE  c  CHR  t  A )  , CHR ( YG ) ) : 

L I N£ ( CHR ( A )  , CHR ( Y  < A)  >  >  j 
END) 

ENDi 

IF  < NOT  UNDEFINED)  OR  <STATUS«3>  THEN 
BEGIN 

MOVE (CHR (FIRST5AMPLE) , CHR ( YG) ) ; 

LINE  iCHR (NPOINTS) ?  CHR ( YG) > ; 

END  | 

END; 

PROCEDURE  DRAWCUR: 

BEGIN 

CUr<DRAW  (OLD  I  X  1  )  j  CUR  DRAW  (OLD  1X2)  ; 

CURORAW ( I  X 1 ) i  CURDRAW < I X2>  | 

ENDI 

(  ****** **************** •  4**************************************************1 


XV 


procedure  rout  1 

gl  oo*l  usaqe:  Kl  ,  Ixi  ,  aCuR ,  x  T  I  MB  ,  / 1  ,  1  v  l  ,  RAW ,  ZEROLEVEL  ,  RBut  ,  » AMF  l 
ouDrxi • old x cur 
r«quir»x  qr aphicft  mod* 

Finds  cursor  posi 1 1 on , ma>  **  this  new  cursor  1 .stores  oi d  v«lue 
repeat*  until  Dutton  is  pressed 


PROCEDURE  ROUT  1 ; 

BEGIN 

UNDEF I  NED: -TRUE: 

DRAWTRACE; 

OLDIX l : -OLDxCUR; 

REPEAT 

DRAWCURj 

CURSOR; 

xi;»  <CUR*X T IM&/2S: 
n : -I FREE C <CUft  3 ; 

I x 1 : *  X  C UR  s  OLD! < I : -OLD*CUF; 
UNTIL  RBUT  1 2Z: 

END; 


procedure  r out 2 

g 1 oD * 1  usages XCUR,X2, I x2 , OLD  1X2, OLDXCUR , XT  I  MB 
V2,RAW, ZEROLEVEL , /AMPl ,RBUT 

same  as  rout  1  but  gets  Cursor  position  and  maps  it  onto  CursorZ 
calculating  the  reel  values  X2,Y2  oF  the  cursor  positon. 
again  the  old  ones  are  saved  For  erasure 


PROCEDURE  R0UT2; 

BEGIN 

UNDEFINED: -TRUE; 

DRAWTRACE; 

OLD I  X 2; -OLDXCUR; 

REPEAT 

DRAWCUR: 

CURSOR; 

X2: »xCUR*XTIMB/25; 

V2: -IFREECXCURI; 

I X2: -XCUR ;  OLD  I X2: -OLDXCUR; 
UNTIL  RBUT :  125; 

END; 


procedure  expFit; 

ptrform*  a  least  squares  Fit  to  the  data  between  the  cursors 
to  the  Function  y«a  exptbx) 


PROCEDURE  EXPFIT; 

VAR 

YS,SX ,SY ,SXX  f SYY ,SXY ,SD»DD1EQ: REAL ; 

S ,  N; INTEGER; 

FUNCTION  Q(X; INTEGER) iREAL; 

VAR 

J  ,  K : REAL; 

Cx INTEGER; 

BEGIN 

J  x  - 1 FREE  C  X  3 ; 

IF  J-  0  THEN  J; -J*  < - 1 >  ; 

IF  J  >0  THEN 
Ki *LN< J) 

ELSE 
K x  » -40  j 

IF  K--40  THEN  Ci-CPM<6,7>;  (*  BEEP  *> 

Qx  -K; 

END; 

FUNCTION  D  <  VAR  Xx  INTEGER) x REAL; 

BEGIN 

Dx  -F I  TA#X+F  I  TB-0  <  X  >  I 
END; 

BEGIN 

FITFOUNDi -FALSE; 

SXx-01  SY:-0;  SXXx-01  SYYx-0;  SXYj«0;  3D: -0; 
FOR  S : - 1  X 1  TO  1X2  DO 


XVI 


BEGIN 

3X:»SX*S; 

SXX:  -3XX«-SQR<S>  ; 
vs i -Q  <S  J ; 

3Y: -SY+VSi 

SYY; -SYY+SQR  < YS)  J 

SXYj -SXY^S*YS{ 

END  j 

Nj -I xr-I XI M  ; 

DD  : -N*SX  X -30ft (SX> | 

FITAj - (N*SXY-SX*SY> /DD; 

FITS: - ($Y*SXX-SX*SXY> /DD; 

ARATEj -FITA*25/XTItt&; 

FOR  Sj«IX1  TO  1X2  DO  SD: -5D+SQR <D<S) > ; 

EQi -SD/ <N-2) ; 

ERARATE* -SORT  <N*EQ/DD> *25/XTIMB; 

WRI TELN '  C  x  ]n  :=  ,SX,  CXXln  ,SxX,  C/]n  ,Sv,  C'Yln  :=  ,  3  Y  Y )  ; 

WRITELN'  CXYln  ,3xv,  C  DD  In  ,SD,  N: *  ,N>; 

WRIT£LN<  nCXXl.N-  EXjn.CXJn  ;*  ,  DD  ,  CDDln/'.N-C)  ' ,EG> ; 

WRITER  a  ,FITA,  b;«  *,FITB,'  -  att.+rtq  :«  , ARATE ,  */- 

WR I TELN  i  ERARA  TE >  ; 

READLNt 

END; 


procedure  rout4 

global  useage :  INTEGRAL, 1X2, 1X1 , RAW , ZERQLEVEL , / AMPL , X T I MB , UNDEF I  NED 
requires  graphics  mode 

displays  data,  calculates  integral  between  X 1  -V  X  2 ,  displays  result, 
on  return  to  program,  negates  undefined  flag 


PROCEDURE  R0UT4; 

VAR 

A: INTEGERS 
BEGIN 

FITFOUND: -FALSE; 

IF  (STATUS  .  4)  THEN 
BEGIN 
DRAWTRACE; 

INTEGRAL! -0; 

DRAWCUR: 

IF  1X2; IX 1  THEN 

FOR  Ai-lXl  TO  1X2  DO  INTEGRAL: *  INTEGRAL* I FREE C A3 
ELSE 

FOR  A:-IX2  TO  I  X 1  DO  INTEGRAL: -  I NTEGRAL+ 1 FREE C A 3 ; 

INTEGRAL: - I NTEGRAL*XT IMB/2S; 

WRITE ( INTEGRAL* 1 . 0E+9: 12: 4 , '  nC ' > ; 

UNDEFINED: -FALSE: 

READLN; 

END 

ELSE 

BEGIN 

WRITELN('  meaningless  *>; 

UNDEFINED: -TRUE } 

READLNi 

END;  (*  of  if  status  <  4  *) 

END; 

(••«««««*«•••**••«*•*•••«•*«••«•«•*«*•••*•«•»*«*»«*«««»*««*««««••***«»«••• •< 

procedure  rout* 

gloal  usage: F IRSTT IME, PZEROLEVEL ,ZEROLEVEL,PARATE,ARATE 
P INTEGRAL , INTEGRAL ,PRAW, RAW, PFRESS, PRESS 
PCLOCKVALUE , CLOCKVALUE , PHUM I D , HUM  I D , F  X  T I MB , X  T \ MB 
P Y AMPL , Y AMPL , PRES  I STQR , RESISTOR , PNRUN , NR UN , PNX PTS , NXPTS 
PNPO I  NTS  ,  NPO  I  NTS  ,  PF  I RST  SAMPLE  ,  F  I  RSTSAPIPLE  ,  PCHANNEL  ,  CHANNEL 
XI  ,  I  X 1 ,Y1 ,X2, 1X2, Y2 

stores  the  parameters  and  gets  more  data  from  file 

i.e.  < f i el d2> i • < f i el d 1 ) i  then  new  (fieldl).  <field2>  is  lost. 


PROCEDURE  ROUTS; 
VAR 

A; INTEGER; 
RTEMP: REAL ; 
BEGIN 


) 


XVI  1 


POL  .-F  i  T FOUND:  *FAl5£: 

CALL.IT :  -FALSE; 

CALL C £> :  -FALSE  ; 

FIRSTTIME: -FALSE: 

<*  SAVE  OLD  PARAMETERS  AND  DATA  *> 

P ZEROLEVEL: -ZEROLEVEL; 

PARATE: -ARATE: 

P INTEGRAL: -INTEGRAL: 

FOR  Ai-0  TO  255  DO  ICQNDC A 3 : - 1 FREE C  A  ]  ; 
PPRESS: -PRESS: 

PCLOCK VALUE: -CLOCKVALUE; 

F MUM ID: -HUMID: 

PXTIMR: -XTIMB; 

PYAMFLj -VAMPL; 

PRESISTOR: -RESISTOR; 

F'NRUN :  -NR UN ; 

PNXPTS; -NXPTS; 

F NF  0 1 NTS : -NPO I NTS : 

PF I RSTSAMFL  E : -F I FST3AMF LE : 

F CHANNEL : -CHANNEL ; 

PF I TFOUND : -F I TFOUND ; 

(*  END  OF  SAVING  PARAMETERS  *> 

FI TFOUND: -FALSE; 

NEWEVENT; 

x i J - 1 x l *  x  T I  MB / 25 : 

Y 1 : *  I FREE C I  X 1 1; 

X2:-IX2*XT I MB/ 25: 

Y2: - IFREEC 1X23; 

UNDEFINED: -TRUE: 

END; 

<»•««•••«•*«**••••«•*«*•*•«»»••••»*•*•••»« 

procedure  rout  7 

global  usage:none 

calls  other  routines  to  draw  (fieldl) 


PROCEDURE  R0UT7; 

VAR 

A,£r:  INTEGER; 

BEGIN 

DRAwTRACE; 

IF  NOT  DMP  THEN  READLN; 
END; 


procedure  r out 3 

global  usage:  I  XI , XCUR, IYl ,RAW,0LDIX1 ,OLDXCUR,Xl , XT  I  MB 
Y 1  ,  ZEROLEVEL , YAMPL , ARATE , &UTTQNFRESSED 
requires  graphics  mode 

gets  a  new  cursorl  value,  calculates  the  attachment  freq. 
here,  and  displays  the  value, 
stops  when  button  is  pressed. 


PROCEDURE  ROUT3; 

BEGIN 

EXPFITj 

F I TFOUND: -TRUE; 

ROUTT; 

END; 

<•••••*••»****••*«••*•««•••*«•**«•*«*«**«*****«»»»»»«»•«««»««•«••««•«»« 
procedure  routS 

global  usaqe: IX1,IX2,1Y1, I Y2 , XCUR , YCUR , X 1 ,X2, Y1 ,Y2 

interchanges  the  cursor  fields;  cur sor l : -cur sor 2 ,  cur sor 2 : -cur sor l 


PROCEDURE  ROUTS; 

VAR 

ITEMPt INTEGER; 

RTEMP : REAL : 

BEGIN 

I  TEMP j -I  XI ;  IX  1 : -I X2:  I X2: -I  TEMP; 
I  TEMP ; •! Y1 |  IYl J -IY2:  IY2;-ITEMP; 
XCUR » - r xi ; 

YCUR :  »•  »  Y  1  j 

RTEMP : -X 1 ;  X 1 ; - X2 ;  X2:=RTEMP; 


<2-.  -RTEMP; 


R  "EMP:  <Mi 
END; 


procedure  rout9 

global  usage: FIRSTTImE  ,  ZEROLEVEL  , P ZEROLEVEL  ,  Aft  ATE  ,  PAR ATE 

INTEGRAL, F INTEGRAL , RAW ,PRAW  ,  PRESS, PPRESS , HUM  ID ,PHUMID 
CLOCK VALUE , PCLOCK  VALUE  ,  XT  I  MB  ,  PXT IMP , VAMPL , PYAMPL 
RESISTOR , PRES I  ST OR , NRUN , PNRUN , NXPTS , PNXPTS 

NPO I NTS , PNPO I N TS , F I RSTSAMPLE  ,  PF I RSTSAMPlE , CHANNEL , PCHANNEL 

XI , I X 1 ,Y1 , X2 , 1 X2,  Y2 

interchanges  tne  dat  a*  i  el  ds;  ( *  a  el  d  1  )  ;  3  i  el  d2  >  ,  <  t  i  el  d2>  :  *  <  f  i  el  d  1  >  . 


PROCEDURE  ROUT 9; 

VAR 

A, I  TEMP: INTEGER; 
ft TEMP: REAL; 

RTEMP: BOOLEAN; 

CTEMP* CHAR; 

BEGIN 

UNDEFINED; -TRUE; 

IF  FIRSTTIME  THEN  ROUTS; 

ft  TEMP:  -ZEROLEVEL;  ZEROLEVEL:  -PZEROLEVEL;  PZEROLEVEL:  =RTEMP; 

RTEMP: -ARATE;  ARATE : -PARATE;  PARATE : -RTEMP ; 

ftTEMP: -INTEGRAL;  INTEGRAL: -P INTEGRAL;  P INTEGRAL :  -RTEMP; 

RTEMP: «FI TFOUND;  F I TFOUND : -PF I TFOUND;  PF I TFOUNDi -BTEMP; 

FOR  Aj-0  TO  255  DO 
BEGIN 

RTEMP : - I FREE  C  A3  ; 

IFREECA3: -ICONDCA3; 

ICONDC  A3 i -RTEMP ; 

END; 

lTEMPs -STATUS;  STATUS : -PST ATUS ;  PSTATUS: - I  TEMP; 

RTEMP; -PRESS;  PRESS; -PPRESS;  PPRESS: -RTEMP; 

RTEMP; -CLOCKVALUE;  CLOCKVALUE ; -PCLDCKVALUE ;  PCLOCKVALUE ; -RTEMP ; 

R .EMP: -HUMID;  HUMID: -PHUM ID;  PHUMID: -RTEMP; 

RTEMP: -XTIMB,  XTImBj -PXTIMB;  PXTIMB* -RTEMP; 

RTEMP: -VAMPL}  VAMPL : -PYAMPL ;  PYAMPL : -RTEMP J 

RTEMP: -RESISTOR;  RES  I  STOP; -PRES  I  STOP;  PRES  I  STOP: -RTEMP; 

I  TEMP: -NRUN;  NRUN; -PNRUN;  PNRUN :- I  TEMP ; 

ITEMP: -NXPTS;  NXPTS: -PNXPTS ;  PNXPTS; - I  TEMP; 

I  TEMP; -NPO I NTS;  NPO INTS; -PNPOINTS ;  PNPO I  NTS; - 1  TEMP; 

I  TEMP: -F I RSTSAMPLE :  F IRSTSAMPLE : -PF  3  RSTSAMPLE;  PF IRSTSAMPLE; - I  TEMP ; 
CTEMP : -CHANNEL ;  CHANNEL : -PCHANNEL ;  PCHANNEL : -CTEMP ; 

XI : -IXleXTlMB/25; 

V 1 1 —  t FREE  C  I X  1  3  ; 

X2; -IX2*XTIMB/2S; 

Y2: -IFREEt 1X23; 

END; 

PROCEDURE  ROUT  10; 

VAR 

At  INTEGER; 

BEGIN 

FOR  At -F IRSTSAMPLE  TO  NPOINTS  DO  IFREE C A 3 : -LN ( SORT ( SQR ( I FREE Z A 3 )  >  >  ; 
END; 


pr oc edur e  rout  U ; 

g  1  ot>*l  usage :  RAW  t  ZEROLEVEL  ,  PRAW  ,  PZEROLEVEL  »  YAMPL  ,  PYAMPL 
performs  the  mapping  described  in  routlB. 
note:  data  in  RAW  is  lost  ' 


PROCEDURE  ROUT  11; 

VAR 

At  INTEGER; 

FACTOR*  t  REAL ; 

BEGIN 

FI TFOUND: -FALSE; 

IF  CORRECTED  THEN  FACTORX ; -STRLEN/ GAPLEN  ELSE  FACTORXj-l; 

UNDEFINED: -TRUE: 

IF  <STATUS-2)  AND ( PSTATUS- 1 )  THEN 
BEGIN 

FOR  A:«FIRSTSAMPLE  TO  NPOINTS  DO  IFREE  C  A  3  :  -  I  FREE  C  A  3 ♦FACTORX* ICONDC  A3 ; 
STATUS* -3; 


XIX 


- V 


END 

ELSE 

IF  tSTATUS=l>  AND  (PSTATUS-»2>  THEN 
BEGIN 

WR I  TELN  i  *  9  "*  )  ; 

READLN 

END 

ELSE 

BEGIN 

W RI  TELN  <  '  <-  5  ‘  )  : 

READLN 

END: 

END) 


procedure  rout  12 

global  usagt:  none 

calls  routS, routS, rout  I  1 ,rout4  ;  a  convenience 


PROCEDURE  ROUT 12* 

VAR 

Ai INTEGER; 

B : REAL; 

BEGIN 

UNDEFINED: -TRUE; 

FITFOUNDi -FALSE; 

IF  <STATUS-3>  AND  (PSTATUS-t)  THEN 
BEGIN 

B: -VE/GAPLEN*100;  <*  I/D  in  metres  •) 

FOR  A: -F IRSTSAMPLE  TO  NPQINTS  DO  I FREE C A ] : - I  FREE C A ] / B ; 
STATUS; -4; 

END 

ELSE 

IF  (STATUS- 1)  AND  'PSTATUS-3)  THEN 
BEGIN 

WR I  TELN  (  \9,r')j 
READLN 
END 
ELSE 
BEGIN 

WR I TELN  <  *<11>  ?T>* ) ; 

READLN 

END;  (*  OF  IF  *> 

END; 

(#*«*•••**#»••*#*««**•*»«•#«»•*»»»«**#«««««#«»««**•*«•«•« 
procedure  rout  13 
global  usage; none 

draws  the  current  field,  outputs  it  to  the  printer 
then  lists  the  parameters  associated  with  it 


PROCEDURE  ROUT  13; 

VAR 

A , B , C: INTEGER; 

BEGIN 

GRSWi -TRUE; 

DMP: -TRUE; 

C: -CPM<5,27> ; 

C ; -CPM (S , 31 >  ; 

C* -CPM  <S ,24) )  <*  SET  EPSON  LINESPACING  FOR  GRAPHICS  *) 

GRAPH I Co ; 

DRAWTRACE; 

DMPRQ ; 

USER ( HF950) )  <*  GRAPHICS  DUMP  ♦) 

C* -CPM  <S , 27 ) ; 

C; -CPM  <5 , 30) ;  (*  RESTORE  DEFAULT  LINESPACING  TO  EPSON  *> 

DMP; -FALSE; 

<*  IF  ENDPAGE  THEN 
BEGIN 

Cs -CPM<3, 12) ; *>  (*  FORM  FEED  *> 

(*  ENDPAGE; -FALSE; 

END 

ELSE  ENDPAGE; -TRUE; *) 

WRJTELN* 


XX 


PQh  Elo'4ElS°,*EBQ6<  :  '  *  RESTORE  D-  P  JUMP  In  BIOS  ’j  CRT:  *> 

GRSW:  -FALSE ; 

END; 


proc«dur*  rout  14 

global  usage: GRSW 

turns  on/ot  f  the  graticule  drawing  routine 


PROCEDURE  ROUT  1 4 ; 

BEGIN 

IF  GR’SW  THEN  GRSW: -FALSE  ELSE  GRSW: -TRUE 
END; 


procedure  rout  IS 
global  usagemone 

reads  in  a  new  record  (ic  and  id>  and  displays  id 


PROCEDURE  ROUT  IS; 

BEGIN 

(*  SET  RETURNED  VALUES  TO  0  * > 
ICMAX : -0; 

IDPOSMAX  t -0; 

IDNEGMAX : *0; 

TX:«0* 

VS: *0; 

DGMAX : -0; 

DQ I NJECTED: =0; 

□RATE: *0; 

DQRATE: -0; 

□ERROR: -0; 

DQERROR: =0; 

TPk:»0: 

DOST  ART : =0 ; 

DQEND: >0; 

INTCi -0: 

INTD2: -0: 

INTD3i -0: 

ROUTS;  ROUTS; 

STATUS: -2; 

PSTATUS: *1 ; 

R0UT7; 

END; 


procedure  rout  16 

global  usageinone 

finds  gap  charge  generation,  dO 


PROCEDURE  ROUT  16; 

VAR 

I CTERH , OMYGODWHATAMESS , RTEMP , RR  TEMP i REAL ; 

A i INTEGER; 

BEGIN 

UNDEFINED: -TRUE; 

IF  (STATUS-4)  AND  (PSTATUS-1)  AND  FITFOUND  THEN 
BEGIN 

OMYGODWHATAMESS : - I CONDCF I RSTSAMPLE  3 ; 

FOR  A! «F I RSTSAMPLE* 1  TO  NPOlNTS  DO 
BEGIN 

\  ICTERMi  -  ( I COND C  A  3  +-OMYGODWHAT AMESS )  * X T I  MB / S0 ; 

OMYGODWHATAMESS i - I COND  C  A  3 ; 

I COND C  A ] : *  I FREE  C A3  - < IFREEC A- 13- ICTERM) *EXP ( ARATE-XT IMB/2S) : 
ICONDCA3 : - IC0NDCA3/XT IMB*2S;  (*  dQ/*ec  ♦  > 

END; 

PSTATUS* -5; 

END 

ELSE  IF  (STATUS- 1>  AND  (PSTATUS-4) 

THEN  BEGIN  WR I TELN <'< 9  >■??')  ;  READLN  END 
ELSE  BEGIN  WRI  TELN  <  '  15>">  ‘  )  ;  READLN  END; 

IF  NOT  FITFOUND  THEN  WRI  TELN  <  '  <Z>'*r* '  )  ; 

PF I TFOUND : -FALSE ; 

I COND  C  F I RSTSAMPLE  3 i -0 1 
F I  RSTSAMPLE :  -F I  RSTSAMPLE-*- 1 ; 


XXI 


Poor?: 

END: 


proctdur*  fOutlT 

form*  the  ttrtamtrf  unct  i  on  A(?:>  in  IFREE 


(•PROCEDURE  ROUT 17; 

VAR 

A,B» INTEGER; 

CC , vELO « D: REAL ; 

BEGIN 

PF I TFOUND: -FALSE: 

IF  (PSTATUS*1)  and  FITFOUND  AND  (VS-  0}  THEN 
BEGIN 
STATUS: -e,; 

VELO: -VE*VS/ < VE+VS) » 

CC: - * SORT ( SOR  <  ARATE > > /VE ) *VELO*XT INB/25: 

HP I TELN*  CC: -  .CC): 

MR  I TELN (  ARATE: -  , ARATE) ; 
mR I TELN  t  VELO: -  . vELO) : 

FOR  A: -  1X1  TO  I  X2  DO  I  FREE C A-I X 1 3 :  «  < I CONDC A3 /VELO) *EXP (  ( A- I  X 1 ) *CC> : 
FOR  A*  -  <  IX2-M-IX1  )  TO  NROINTS  DO  I  FREE  C  A]  :  -0; 

END 

else  wrjteln<  •:  is  .**  v3  ^  •  12, 

READLN; 

FITFOUND: -FALSE; 

F-FI  TFOUND:  -  TRUE; 

END:  *> 


procedure  routlS; 

finds  the  streamer  velocity 

t*«««««***»«*«*«***««***«**«t*«»»«<»#*«#**««««*«*«*»4****««4t«44*»#*4«««4«*) 

PROCEDURE  ROUT 18» 

BEGIN 

VS: -STRLEN/ ( ABS ( X2-X 1 ) *1 00) \ 

END: 

PROCEDURE  ROUT  19: 

BEGIN 

POLYFIT* 

R0UT7; 

END: 


procedure  rout20 

global  usage:  CORRECTED; 

changes  ROUTil  such  that  it  finds  le-Id+Ic;  and  not  the  default 
which  is  Ie-Id* (d-n) /d. Ic 


PROCEDURE  RCUT20: 

BEGIN 

IF  CORRECTED  THEN  CORRECTED: -FALSE  ELSE  CORRECTED: -TRUE : 
IF  CORRECTED  THEN  WRI TELN ( * f i nds  Ie-Xd+Ic. (d~n > /d ' ) 

ELSE  WRITELN<  finds  Ie-Id+Ic ' > ; 

READLN: 

END; 

PROCEDURE  ROUT21: 

VAR 

A: INTEGER: 

BEGIN 
I CMAX : -Y2j 
TPK : - <  X2-X l )  | 

INTC: -0: 

FOR  A: -F1RSTSAMPLE  TO  NPOINTS  DO  INTC: - INTC* IFREEC AJ ; 
INTC:«INTC*XTIMB/23| 

END; 

PROCEDURE  ROUT22: 

BEGIN 

IDPOSMAX : -Y2; 

END: 

PROCEDURE  ROUT23; 

VAR 

A: INTEGER; 

BEGIN 


XXI 1 


CALU2'!  s- TRUE.: 

IDnEGMAXs  «v2; 

TXi - ( XZ-t 1 > ; 

ROUT l 3; 

INTD2* -05 
INTD3* -0; 

POP  A i  —  ( I  X 1 )  TO  1X2  DO  INTD2: - INTD2* IFREEt A3 ; 

FOR  A:  -  (  I  X  1  )  TO  NFOINTS  DO  I NTD3 :  -  I  NTD3+ 1  FREE  C  A  ]  ; 

INTD2: - INTD2*XT ln&/ 25: 

INTD3: * INT03*XT I MB/ 25: 

END* 

PROCEDURE  ROUT 24 ; 

BEGIN 
QGMAX  j - Y 1 1 
QRATE: -ABATE; 

QERROR: -ERARATE; 

END; 

PROCEDURE  ROUT 25; 

BEGIN 

DQSTART ; - V 1 ; 

DQEND: -Y2; 

DORATEx -ARATE; 

DOERROR : -ERARATE ; 

END; 

PROCEDURE  ROUT 26; 

BEGIN 

CALL26: -TRUE; 

DQ INJECTED; -V2; 

END; 

PROCEDURE  R0UT27; 

VAR  A: INTEGER; 

BEGIN  • 

WRITELN (DBASE, ’  >; 

WRITE  (DBASE,  '  \  ,FIl_ENO  ,  ',NXPTS,' 

WRITE (DBASE, FI  ELD ,  •  ‘  ,  (FIELD/PRE)  ,  '  ',GAPLEN,  ',PR£,'  ')j 

WR I TE ( DBASE , I CMAX ,  0  ’  ,  I  DP OSMAX ,  *  0  '  , I DNEGMAX ,  '  0  >; 

WRITE ( DBASE , QGMAX  ,  0  ’ , DO I NJECTED ,  0  .ORATE,'  '  , QERROR,  '  >; 

WRITE (DBASE, DOB ATE, '  ' , DQERROR , '  ',VS,'  0  ,  VE,'  0  ); 

WRITE (DBASE, INTO, '  0  ' , INTD2, '  0  ' , INTD3, *  0  >; 

WRITE (DBASE, VE, '  0  ' , TX , '  0  ,STRLEN,  0  ',TPK,  0  ', DQSTART,'  0  '>; 
WRITE (DBASE, DQEND, '  0  '); 

IF  CALL23  THEN  WRITE  ( DBASE  ,  '  1  ')  ELSE  WR  I  TE  ( DBASE  ,  '  0 

IF  CALL26  THEN  WR  I  TE  ( DBASE  ,  '  1  ')  ELSE  WR  I  TE  ( DBASE  ,  '  0  '); 

END; 


proc«dur»  analysis 
global  u**qt: IRESP 

drives  the  sub-menu  for  the  cuh  sor -or  x  anted  analysis  package 

PROCEDURE  ANALYSIS; 

BEGIN 

REPEAT 

NORMALSCREEN; 

(*<*  WRITELN< ' YOURE  SO  SMART  !');*) 

WAIT (2) | 

NORMALSCREEN; 

ANALMENU; 

READ ( IRESP) ; 

CASE  IRESP  OF 
1  *  ROUT  1 ; 

2 l ROUT 2; 

3i RQUT3; 

4 i R0UT4 ; 

5*  ROUT  15; 

<*  6; R0UT6; * ) 

7iROUT7; 

8 i ROUTS; 

9*  R0UT9; 

10i ROUT  10; 

1 IjROUTI 1 | 

1 2j  ROUT 1 2 1 
1 3; ROUT 1 3 i 


1 4 : RGu  r 1 4  j 
1  6:  ROUT lb; 

(*  l 7s ROUT l 7*  *- 

lasftouTia: 

1 9 : ROUT 1 9 ; 

20s  ROUT20; 

2i : ROUT 21 ; 

22  s  ROUT 22  s 
27: ROUT27; 

24 ; R0UT24 ; 

25*  R0UT25? 

26  s  ROUT 26* 

27  s  R0UT27 
END; 

UNTIL  1RESF-6; 

END;  i *  of  analysis  • 


program  begins  “e^e 


&EGIN  t  *  program  *» 
SVCCOOEi 
INITIALISE; 
BEGINNING: 
ANALYSIS; 

END.  <*  of  pregram 


Appendix  B 


A  complete  data  list  is  available  from  the  author. 

The  data  has  not  been  included  in  order  to  reduce  the  bulk  of 


the  report 


END 

DATE 

FILMED 


